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ABSTRACT 


Nakkash,  Tammam  Zaki.   Ph.D.,  Purdue  University,  June, 
1969.   ACTIVITY-ACCESSIBILITY  MODELS  OF  TRIP  GENERATION. 
Major  Professor:   William  L.  Grecco. 


This  research  examined  the  trip  generation  process 
with  the  specific  purpose  of  evaluating  the  effect  of 
activity-accessibility  variables  on  trip  generation.   An- 
other aspect  of  accessibility  was  studied  by  stratifying 
the  zones  of  the  study  area  by  location.   In  contrast  to 
the  traditional  trip  generation  procedures  where  the  trips 
generated  by  a  zone  are  considered  to  be  a  function  of  the 
characteristics  of  only  the  zone  itself,  the  models  pro- 
posed by  this  research  take  into  account  also  the  character- 
istics of  all  the  other  zones  and  the  transportation  net- 
workconnectingthein. 

Operationally,  the  accessibility  of  a  zone  to  each 
activity  was  defined  as  the  sum  of  the  products  of  the  size 
of  each  activity  in  each  zone  of  the  study  area  and  the 
friction  factor  corresponding  to  the  travel  tine  between 
each  zone  and  the  zone  under  consideration.   The  friction 
factors  were  those  derived  from  calibrating  a  gravity  model 
of  trip  interchange  for  various  trip  purposes.   The  access- 
ibility of  each  zone  was  then  expressed  as  a  percentage  of 
total  accessibilities  of  all  zones  of  the  area  and  referred 


XI  V 


to  as  relative  accessibility.   Criteria  v;ere  also  set  to 
differentiate  the  zones  of  the  study  area  into  two  strata: 
central  and  non  central  areas. 

Utilizing  data  from  the  1964  surveys  for  a  comprehen- 
sive transportation  study  in  Indianapolis,  Indiana,  four 
sets  of  trip  generation  regression  equations  were  developed 
for  each  of  thirteen  trip  purposes.   One  of  the  developed 
sets  was  a  recall  brat  ion  of  the  Indianapolis  Regional  Trans- 
portation and  Development  Study  (IRTADS)  trip  generation 
equations.   This  set  was  developed  by  the  traditional  pro- 
cedures, that  is,  the  independent  variables  were  limited  to 
demographic,  socio-economic,  and  land  use  variables.   A 
second  set  was  developed  in  which  relative  accessibility 
variables  were  included  among  the  independent  variables. 
Two  more  sets  were  developed,  each  corresponding  to  one  of 
the  former  two  sets  but  calibrated  with  data  stratified 
according  to  the  zone's  location  in  the  central  or  non 
centra  1  area . 

The  developed  models  were  compared  as  to  their  statisti 
cal  strength.   Considering  the  factor  of  location  always 
improved  the  statistical  strength  of  the  trip  generation 
models.   Including  relative  accessibility  variables  im- 
proved the  statistical  strength  of  trip  attraction  models 
more  than  that  of  trip  production  models.   Stratification 
alone  improved  the  models  more  than  including  relative 
accessibility  variables  only. 
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1985  forecasts  of  the  demographic,  socio-economic,  and 
land  use  variables  together  with  estimates  of  travel  time 
on  the  proposed  future  network  were  inputs  to  the  four  sets 
of  developed  models  to  forecast  trip  productions  and  attrac- 
tions.  Comparison  of  these  forecasts  indicated  that  strati- 
fied models  forecasted,  on  the  average,  more  trips  for 
zones  of  the  non  central  area  and  less  for  zones  of  the 
central  area  compared  to  forecasts  made  by  unstratified 
models.   Forecasts  by  models  with  accessibility  variables 
and  stratification  were  significantly  different,  on  the 
average,  than  forecasts  by  basic  IRTADS  models.   There  were 
some  trip  purposes  for  which  IRTADS  models  indicated,  on 
the  average,  larger  trip  forecasts  than  the  accessibility 
models  with  stratification.   Further  analysis  on  those  trip 
purposes  showed  that  the  zones  where  the  IRTADS  models  did 
not  overforecast  trips  were  in  the  vicinity  and  along  cor- 
ridors defined  by  the  major  thoroughfares  of  the  study  area 

The  results  of  this  research  recommend  that  future 
relative  accessibility  of  zones  of  the  study  area  be  con- 
sidered in  estimating  future  trips.   The  process  would  be 
iterative  and  would  be  terminated  when  equilibrium  is 
reached  between  the  forecasted  demand  for  transportation 
(future  trips)  and  planned  for  supply  of  transportation 
(  future  network ) . 


CHAPTER  I.   INTRODUCTION 

In  recent  years  one  has  witnessed  a  great  increase  in 
urban  transportation  planning  activities  in  many  American 
cities.   This  increased  activity  was  in  response  to  the 
challenging  urban  transportation  problem,  one  of  the  major 
urban  problems  of  contemporary  cities.   The  safe  and  effi- 
cient intracity  movement  of  goods  and  people  is  very 
essential  for  the  social,  cultural,  and  economic  health  of 
an  urban  area . 


Urban  Transportation  Planning 
The  urban  transportation  problem  is  the  product  of 
many  interacting  factors.   The  enormous  population  growth 
in  urban  areas  and  their  expanding  areal  extent  as  a  result 
of  the  redistribution  of  population,  the  imoroved  standard 
of  living  due  to  increased  affluency,  and  the  subsequent 
greater  reliance  on  private  automobiles  are  only  some  of 
those  factors.   Together  with  those  size-related  features 
of  the  problem,  the  temooral  aspects  induce  periodic  high 
demands  for  transportation.   This,  of  course,  is  due  to  the 
interdependency  of  human  activities  which  simultaneously 
occur,  essentially,  during  the  eight-hour  work-day;  starting 
and  ending  at  rather  definite  times. 


A  recognition  of  the  immense  complexity  and  the  vast 
dimensionality  of  the  urban  transportation  problem  is  a 
prerequisite  to  any  attempt  at  solving  the  problem. 

The  symptoms  of  the  urban  transportation  problem  were 
first  most  apparent  in  the  congestion  of  the  streets;  thus 
early  attempts  at  reducing  the  severity  of  the  situation 
VI ere    traffic  volume  oriented.   The  failure  of  projections 
based  on  traffic  volumes'  factoring  to  cope  with  the  need 
for  accurate  forecasts,  to  be  used  in  highway  planning  and 
design,  has  directed  the  effort  of  transportation  planners 
since  the  1940's  toward  more  sophisticated  approaches.   The 
development  of  Urban  Transportation  Planning  (UTP)  was  a 
step  in  thi  s  di  recti  on . 

UTP  views  the  traffic  condition  of  the  city  streets 
as  an  element  within  the  framework  of  the  larger  urban 
environment.   Projections  for  transnortatinn  planning  in- 
clude in  addition  to  the  expected  amount  of  travel  its 
temporal  and  locational  pattern;  and  consider  all  transpor- 
tation modes.   UTP  is  a  continuing  process,  projections 
and  plans  are  periodically  adjusted  to  conform  to  changes 
in  the  complex  picture  of  the  urban  environment-. 

UTP  has  developed  into  an  intricate  process  which 
relies  heavily  on  surveys.   An  enormous  body  of  data  is 
collected,  usually,  on  a  sampling  basis.   The  data  include 
information  on  the  demographic,  social,  and  economic 
characteristics  of  the  area.   An  exhaustive  inventory  of 


the  existing  transportation  facilities  and  their  character- 
istics is  also  conducted.   The  land  use  survey  supplies 
the  information  on  the  quantity,  intensity,  and  location 
of  land  classified  by  use.   Detailed  floor  area  information 
are  collected  for  commercial,  services,  and  educational 
uses  . 

The  backbone  of  the  urban  transportation  planning 
surveys  is  the  Origin-Destination  (0-D)  Study.   This  study 
consists  of  a  home- interview,  a  truck-taxi  survey,  and  a 
cordon  count.   The  results  of  this  study  establish  the 
travel  patterns  in  the  base  year  and  provide  the  data  in- 
puts for  model  formulation. 

The  socio-economic,  demographic,  and  land  use  informa- 
tion are,  usually,  summarized  by  zone.   The  zone  is, 
usually,  the  smallest  areal  unit  considered  in  all  further 
analyses  and  projections. 

The  transportation  network  is  coded  into  a  topological 
network  of  links  and  nodes.   Access  to,  and  egress  out  of 
the  network  is  accomplished  via  zone  centroids.   The  topo- 
logical network  carries  all  the  information  of  the  physical 
transportation  network  relevant  to  the  planning  process. 
This  abstract  numerical  representation  of  the  road  and 
transit  network  is  amenable  to  mathematical  manipulation. 

The  travel  information  obtained  from  the  home- i ntervi ew , 
and  truck-taxi  survey  are  utilized  to  explain  trip  genera- 
tion of  the  zones  of  the  study  area.   The  socio-economic 


and  land  use  characteristics  of  the  zone  are  considered  to 
be  the  explaining  variables. 

A  model  of  zonal  trip  interchange  is  calibrated  uti- 
lizing the  characteristics  of  the  transportation  network 
and  the  present  travel  pattern.   The  parameters  of  the 
calibrated  model  reflect  the  propensity  of  each  zone  to 
interchange  trips  with  the  other  zones,  and  the  rate  of  de- 
cay of  this  trip  interchange  with  increased  zone-to-zone 
separation,  measured  in  travel  cost,  travel  time,  or  travel 
distance.   These  relationships  are  assumed  either  to  stay 
the  same  for  the  future,  or  to  change  in  a  predictable 
manner. 

The  distribution  of  trips  is  then  assigned  to  the 
transportation  network  yielding  link  volumes.   A  successful 
assignment  v/hich  reasonably  matches  ground  counts  in  the 
base  year  indicates  that  a  similar  assignment  of  forecast 
trip  interchanges  should  be  valuable  in  indicating  defi- 
ciencies . 

Socio-economic,  demographic,  and  land  use  data  are  then, 
separately,  forecasted  for  the  future  design  year.   This 
provides  the  basis  for  a  future  trip  generation  forecast. 
Future  trip  generation  is  an  input  to  estimate  future  zonal 
interchanges  using  the  calibrated  distribution  model.   Fu- 
ture trip  assignment  is  the  basis  for  identifying  possible 
future  deficiencies,  evaluating  proposed  networks,  and 
proposing  alternative  plans.   Figure  1  shows  the  major  steps 
in  UTP. 
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FIGURE    I    -     URBAN     TRANSPORTATION      PLANNING 


Trip  Generation  in  the  Urban  Transportation  Planning  Process 

Apart  from  the  socio-economic,  demographic,  and  land 
use  forecasts  trip  generation  constitutes  the  first  step 
toward  establishing  the  future  travel  pattern.   The  accuracy 
of  the  future  trip  distribution  in  forecasting  design  year 
trip  interchange  cannot  be  any  better  than  the  accuracy  of 
the  trip  generation  forecasts,  except  due  to  chance. 

The  ultimate  purpose  of  the  trip  generation  analysis 
is  to  arrive  at  an  estimate  of  the  trip  ends  generated  at 
each  analysis  unit  of  the  study  area.   Trip  generation  tech- 
niques try  to  establish  a  relationship  between  the  demo- 
graphic and  socio-economic  characteristics  of  the  population 
of  an  analysis  unit  and  its  trip  generation.   Similarly, 
the  intensity,  character,  and  location  of  different  land 
uses  are  related  to  trip  making  of  the  analysis  units.   These 
procedures  are  based  on  the  hypothesis  of  a  causal  relation- 
ship between  population  characteristics,  land  use,  and  the 
trip  making  behavior  of  people. 

Traditionally,  trip  generation  forecasts  are  estab- 
lished independently  of  any  direct  consideration  to  the 
transportation  network.   This,  of  course,  assumes  that  trips 
produced  at,  or  attracted  to  a  zone  are  a  function  only  of 
the  attributes  of  the  zone  itself;  and  are  not  directly  a 
function  of  the  transportation  network  which  provides  the 
roadway  for  the  trips.   Apart  from  some  consideration  of 
the  future  network  in  forecasting  land  uses,  this  ignores 
the  feedback  and  interaction  between  some  significant 


components  of  the  complex  system  of  the  urban  environment. 
A  flow  chart  of  the  traditional  trip  generation  process  is 
presented  in  Figure  2. 

Purpose  and  Scope 

The  purpose  of  this  research  was  to  study  the  trip 
generation  process,  with  the  specific  aim  of  investigating 
the  feedback  from  the  transportation  system  on  the  rate  of 
trip  making.   Conceptually,  there  is  no  strong  basis  for 
assuming  that  trip  making  is  independent  of  the  transporta- 
tion system.   On  the  contrary,  it  seems  that  trips  produced 
by,  or  attracted  to  a  zone  should  be  a  function  of  the 
relative  accessibility  of  the  zone  to  different  land  uses, 
in  addition  to  the  characteristics  of  the  zone  itself. 

Trio  making  is  a  product  of  the  desire  for  human 
interaction.   Within  present  day  technology,  way-of-life, 
and  the  requirement  of  compatible  land  uses,  different  daily 
activities  have  to  be  performed  at  different  locales.   The 
possibilities  and  merits  of  substituting  travel  by  improved 
communications  are  fully  recognized  [1].    However,  it  is 
not  anticipated  that  such  substitution  will  drastically 
modify  travel  needs  and  habits,  in  our  cities,  within  the 
next  few  decades  . 

Basically,  the  rate  of  trip  making  is  a  function  of 
two  categories  of  variables.   One  category  represents 
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variables  which  tend  to  increase  the  potential  of  trip 
making;  the  other  constitutes  the  restrictive  forces.   Vari- 
ables such  as  the  availability  of  the  vehicle  to  the  resi- 
dents of  a  zone,  and  the  percent  participation  of  the  resi- 
dents in  the  labor  force;  or  the  number  employed  in  a  zone, 
and  the  amount  of  floor  area  of  different  land  uses  in  it 
are  examples  of  the  first  category  of  variables.   They 
measure  respectively  the  potential  of  trip  production  or 
trip  attraction.   The  penalties   incurred  by  travel  measured 
in  cost,  travel  time,  or  travel  distance  represent  the 
variable  which  belong  to  the  second  category,  the  restrictive 
forces  . 

This  study  utilized  data  obtained  from  the  surveys  for 
the  Indianapolis  Regional  Transportation  and  Development 
Study  (IRTADS).   Multiple  linear  regression  predictive  models 
of  person-trip  productions  and  attractions,  by  purpose,  were 
developed.   The  developed  models  differ  from  the  traditional 
trip  generation  models.   The  independent  variables  were  not 
restricted  to  socio-economic  and  land  use  measures  of  the 
zones,  but  included  also  measures  of  the  relative  accessi- 
bility of  the  zone  to  different  activities  and  land  uses. 

The  locational  aspects  that  affect  trip  generation  were 
also  investigated.   It  was  hypothesized  that  central  loca- 
tions in  the  study  area,  generally,  afford  greater  accessi- 
bility; and  the  convergence  of  the  street  network  on  the 
city  center  favors  the  core  location.   The  zones  of  the 
study  area  were  stratified  into  two  groups:   a  central  and 
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a  non  central.   This  stratification  was   entered  as  an 
independent  dummy- vari abl e  in  the  trip  generation  analysis 

A  comparison  of  the  forecast  trip  generation  by  the 
suggested  approach  with  the  forecast  by  the  traditional  ap- 
proach was  conducted. 


1 1 


CHAPTER  II.   REVIEW  OF  LITERATURE 


General 

The  ultimate  purpose  of  the  trip  generation  analysis  is 
to  establish  a  reasonable  estimate  of  the  trips  produced  at, 
or  attracted  to  some  subunit  of  the  urban  area,  for  the 
design  year.   This  problem  has  been  approached  differently 
by  various  transportation  studies.   Due  to  the  difficulty 
of  projecting  trips  directly,  it  has  been  general  practice 
to  establish  some  relationship  betv/een  trips  and  other 
attributes  of  the  urban  area  and  to  use  this  relationship 
to  provide  a  forecast. 

The  three  main  approaches  that  have  been  used  d.re: 
(1)  trip  rate  analysis,  (2)  cross-classification  analysis, 
and  (3)  multiple  regression.   These  three  approaches  and 
their  variations  will  be  discussed  later  separately. 

It  is  significant  to  observe  that,  as  yet,  "The  Model" 
of  trip  generation  is  not  available.   That  is,  there  is  no 
simple  model  or  set  of  models  which  have  well  defined 
structure  and  mathematical  formulation  that  is  applicable 
to  any  study  area   when  calibrated  with  data  obtained  from 
it.   Apparently,  a  valid  conceptual  basis  for  trip  genera- 
tion has  not  yet  been  reached. 
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Trip  generation  analysis  is  usually  conducted  v/ith 
reference  to  some  subunit  of  the  study  area.   In  this  case, 
and  since  the  0-D  study  furnishes  information  on  the  dwell- 
ing unit  level,  some  aggregation  of  data  has  to  be  done.   A 
valid  model  of  trip  generation  should  predict  the  same 
total  number  of  generated  trips,  regardless  of  the  way  the 
study  area  is  subdivided,  or  the  level  of  data  aggregation 
[2].   The  subdivision  of  the  study  area  into  subunits,  such 
as  traffic  zones,  could  be  done  in  an  infinite  number  of 
ways.   Urban  geographers  have  recognized  this  fact 

...any  set  of  area!  boundaries  is  essentially 
arbitrary...  When  one  computes  the  density  of 
an  area!  unit,  he  obtains  a  figure  which,  though 
well  defined,  is  necessarily  relative  to  areal 
boundaries  and  hence  fictitious  in  a  sense 
[3,  p.  37]. 

When  dividing  the  study  area  into  zones,  the  interest 

will  be  in  the  among-zones  relationship,  rather  than  the 

variation  within  the  zone.   The  assumption  of  homogeneity 

within  the  zone  is  of  questionable  validity.   Supposedly, 

the  smaller  the  zone,  the  better  the  chances  that  this 

assumption  is  valid;  but  rather  the  location  of  the  boundary, 

not  only  the  size  of  the  zone,  will  affect  its  homogeneity 

[3]. 


Trip  Generation  Analysis  Procedures 
The  procedures  used  in  trip  generation  analysis  attempt 
to  establish  a  relationship  between  the  amount  of  trips  gen- 
erated per  zone,  and  the  various  demographic,  economic  and 
land  use  characteristics  of  the  zone.   They  differ  only,  in 
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the  manner  in  which  the  relationship  is  established. 


Land  Area  Trip  Rate  Analysis 

The  method  utilizes  the  numerical  relationship  between 
land  use  and  trips  which  it  generates  [4].  The  rate  of 
trips  generated  by  land  areas  of  different  character,  loca- 
tion, and  intensity  are  established.  Character  of  the  land 
area  denotes  its  functional  use,  location  refers  to  the 
spatial  arrangement  within  the  urban  area,  and  intensity  of 
use  is  usually  measured  by  density. 

Trip  ends  data,  obtained  from  the  0-D  study,  are  related 
to  land  use  information  to  establish  the  trip  rate  per  acre 
for  each  different  land  use.   The  procedure  was  first  used 
and  documented  by  the  Chicago  Area  Transportation  Study 
(CATS),  after  it  was  originally  introduced  by  the  San  Juan 
Study  [5].   The  experience  of  CATS  with  this  procedure  indi- 
cated that  a  forecast  of  the  trips  by  applying  the  trip 
rates  of  the  base  year  to  the  land  use  forecast  of  the 
future  year  resulted  in  an  underestimate.   This  discrepancy 
resulted  because  the  trip  rate,  itself,  is  a  function  of 
other  variables  which  were  not  considered.   The  low  estimate 
had  to  be  adjusted  by  a  control  total  established  by  a 
second  projection  based  on  car  ownership  and  residential 
density  [6].   An  example  of  trip  rates  by  land  area  for  dif- 
ferent land  uses  is  presented  in  Table  1. 
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Cross-Classification  Analysis 
Cross-classification  is  a  technique  applicable  when 
large  numbers  of  observations  are  available.   The  technique 
requires  data  to  be  summarized  in  cross -tabu  1 ati ons  ,  so 
that  the  frequencies  or  the  average  values  of  a  dependent 
variable  can  be  obtained  for  different  combinations  of 
values,  by  classes  or  ranges,  of  two  or  more  independent 
variables  [7].   This  type  of  analysis  does  not  rely  heavily 
on  any  assumptions,  about  the  distribution  of  the  variables. 
Thus,  it  could  be  thought  of  as  a  nonparametri c  technique 
[8]. 

Cross-classification  has  been  used  as  a  research  tool 
in  the  analysis  of  trip  generation.   Oi  and  Shuldiner  re- 
ported cross-tabulations,  such  as  trips  per  household  at 
various  levels  of  car-ownership  and  persons  per  dwelling 
unit  [9].   Recently,  cross-classification  is  gaining  recogni- 
tion by  operational  transportation  studies  [10,11,12]. 

A  modification  of  the  cross-classification  analysis 
was  developed  at  the  Puget  Sound  Regional  Transportation 
Study  (PSRTS).   Later  it  was  used  by  the  Albuquerque 
Transportation  Study  to  compare  results  with  a  regression 
procedure  forecast  of  person-trips  [12].   For  home-based 
person-trip  productions,  trips  per  household  were  calculated 
Analysis  zones  of  similar  household  and  enviromental  char- 
acteristics were  grouped  together,  and  average  trip  rates 
per  household  were  calculated  for  each  group  of  zones.   In 
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selecting  the  independent  variables,  results  of  previous 
research,  and  variability  of  the  data  were  the  criteria. 

The  zones  were  ranked  according  to  the  dependent  vari- 
able from  high  to  low,  similarly,  they  were  ranked  according 
to  each  of  the  independent  variables.   A  Spearman's  rank 
correlation  coefficient  rho,  of  the  ranking  based  on  the 
dependent  variable  as  compared  to  each  of  the  rankings  ac- 
cording to  each  of  the  independent  variables  was  calculated. 

Average  person-trips  per  household  were  cross-tabulated 
at  three  levels  of  intensity  of  environmental  characteris- 
tics and  six  levels  of  household  characteristics.   The  18- 
cell  table  was  referred  to  as  the  rank-classification  matrix, 
and  is  shown  in  Figure  3 . 

Projections  of  the  independent  variables  for  the  fore- 
cast year  were  developed  at  the  zone  level.   In  order  to 
establish  a  forecast  of  the  trips  produced  at  any  zone,  one 
must  use  the  rate  of  trips  per  household  in  the  cell  of  the 
rank-classification  matrix  corresponding  to  the  forecasted 
environmental  and  household  characteristics  of  the  zone. 
An  example  is  shown  in  Figure  4. 

A  similar  approach  was  used  to  project  trip  attractions. 
It  is  significant  that  some  zones  had  a  small  employment 
base,  therefore,  in  order  to  insure  statistical  stability 
of  employment  trip  generation  rates  "...zones  for  which  the 
trip  information  indicated  less  than  a  selected  level  were 
omitted  from  the  analysis  for  determining  the  groupings  of 
zones ."  [12,  p.  95]. 
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Cross-classification  is  simple,  does  not  rely  heavily 
on  assumptions  underlying  the  distribution  of  the  variables, 
keeps  the  study  staff  closer  to  the  data;  but  it  has,  how- 
ever, serious  disadvantages.   Large  number  of  observations 
are  needed  for  conclusive  results.   The  number  of  observa- 
tions needed  increases  with  the  increase  in  the  dimensional- 
ity of  the  classification  matrix.   Major  shortcomings  of 
this  technique  were  stated  as  follows  by  Ezekiel  and  Fox 

...it  provides  no  measure  of  how  important  the  rela- 
tion shown  is  as  a  cause  of  variation  in  the  factor 
being  studied,  or  of  how  closely  that  factor  may  be 
estimated  from  the  others  on  the  basis  of  the  rela- 
tions shown... [cross-classification]... does  not 
determine  the  relationships  where  many  variables  are 
involved  so  satisfactorily  as  does  multiple  regres- 
sion [7,  p.  394]. 

Multiple  Regression  Analysis 

General    Multiple  regression  is  a  statistical  technique 
to  fit  a  response  surface  satisfying  the  least  squares  cri- 
terion.  The  increased  use  of  multiple  regression  could 
probably  be  attributed  to  the  recent  availability  of  computer 
codes  that  can  handle  the  estimation  of  regression  parameters 
efficiently. 

Generally  the  functional  relationship  between  the 
dependent  and  independent  variables  could  be  represented 
mathematically  by 


Y  -  f(X-|,  ...,  XlB-j,  . 


) 
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where 

Y   =  the  dependent  (or  response)  variable 

X.  =  the  ith  independent  variable  (i=l,  ...,  p) 

B.  =  the  jth  parameter  in  the  function  (j=l,  ...,  q) 

and  f  stands  for  the  assumed  form  of  the  function 

[13]. 

In  order  to  make  probability  statements  on  the  parameter 
estimates  three  basic  assumptions  have  to  be  satisfied: 

1.  the  independent  variables  are  fixed  and  measured 
wi  thout  error ; 

2 .  for  each  given  value  of  the  independent  variables, 
the  values  of  the  dependent  variable  are  normally 
and  independently  distributed;  and 

3.  the  variances  of  the  distributions  of  the  dependent 
variable  are  all  the  same  [13]. 

In  a  trip  generation  analysis  based  on  multiple  regres- 
sion the  previous  assumptions  dire    likely  to  be  violated. 
Statistical  tests  are  available  to  check  if  the  second  or 
third  assumption  is  violated.   Some  transformations  of  the 
data  could  be  used  to  correct  for  a  violated  assumption. 
The  topic  of  regression  is  well  documented  in  the  literature 
on  statistics  [7,13,14]. 


Linear  Versus  Nonlinear   Regression .   The  distinction 
between  linear  and  nonlinear   regression  is  in  the  functional 
form  of  the  model.   Linear  regression,  as  the  name  indicates, 
refers  to  regression  models  which  are  linear  in  the 
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parameters.   Some  models  with  exponential  parameters  are 
intrinsically  linear  since  they  can  be  put  in  linear  form 
by  a  transformation  [14].   The  estimation  of  the  parameters 
of  a  linear  model  is  reached   by  solving  a  set  of  simultane- 
ous linear  equations,  referred  to  as  the  Normal  Equations. 
The  estimation  of  the  parameters  of  a  nonlinear  model  var- 
ies in  complexity  with  the  form  of  the  model.   Search  tech- 
niques are  usually  utilized  for  this  purpose. 

Linear  regression  has  been  more  popular  than  nonlinear 
regression  in  trip  generation  analysis.   Linear  regression 
is  simpler,  and  does  not  require  ai  priori  knowledge  of  the 
exact  functional  form  of  the  model.   One  of  the  earliest 
documentations  of  the  use  of  linear  regression  was  by  Wynn 
[15].   Nonlinear  regression,  however,  has  the  advantage  of 
forcing  the  researcher  to  look  at  the  relationship  between 
the  variables.   Nonlinear  regression  was  used  to  estimate 
recreational  trip  rates  in  research  reported  by  ^'atthias 
[16]. 


Constrained  Regression  Estimation .   The  assumptions  of 
regression  state  that  the  independent  variables  are  fixed 
and  uncorrelated.   The  independent  or  explanatory  variables 
in  a  regression  analysis  of  trip  generation,  in  general, 
violate  these  assumptions.   In  the  extreme  case,  where  a 
very  high  correlation  exists  between  some  of  the  independent 
variables,  the  computational  procedure  fails,  since  this 
implies  matrix  singularity  and  the  Normal  Equations  do  not 
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have  a  unique  solution.   If  this  occurs,  enough  warning 
vn'll  be  received.   A  1  ess-than-perfect  correlation  between 
some  of  the  independent  variables  will  result  in  dividing 
their  effect  on  the  dependent  variable  rather  arbitrarily. 
Many  times  this  yields  conceptually  invalid  parameter  esti- 
mates (regression  coefficients).   If  theory  can  provide 
some  bounds  on  the  values  of  some  of  the  coefficients  of 
the  independent  variables,  then  it  seems  appropriate  to  use 
this  information  in  developing  the  regression  models;  thus 
avoiding  structural  mis-specification.   Finding  the 
parameters'  estimates  of  a  constrained  regression  model  is  a 
quadratic  programming  problem.   Meyer  and  Glauber  provide 
the  basic  mathematical  formulation  for  this  problem  [17, 
pp.  184-196].   Recently  reported  research  utilized  this 
technique  in  estiriating  urban  passenger  travel  behavior  [2]. 


The  Uses  of  Dummy-Va ri abl es  in  Regression.   Some  of  the 
independent  variables  that  are  relevant  to  trip  generation 
are    impossible  to  measure  on  a  continuous  scale.   Cross- 
classification  is  one  of  the  useful  techniques  in  such  a 
case,  however,  it  has  the  disadvantages  stated  earlier. 
The  dummy- vari abl es  technique  handles  this  case  in  the  frame 
work  of  multiple  regression.   In  the  context  of  trip  genera- 
ation  household  characteristics  could  be  stratified  with 
respect  to  various  factors,  with  different  levels  for  each 
factor.   Each  household  characteristic  will  have  a  dummy- 
variable  associated  with  it,  except  one  class  to  avoid 
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having  more  coefficients  to  estimate  than  independent  Normal 
Equations.   The  dummy-vari abl e  associated  with  each  class 
will  have  a  value  of  zero  except  for  the  one  class  to  which 
the  observation  belongs  will  have  a  value  of  one.   The 
earliest  reference  in  the  literature  of  statistics  on  the 
use  of  dummy- vari abl es  is,  probably,  the  article  by  Suits 
[18]. 

The  use  of  dummy- vari abl es  in  trip  generation  studies 
is  gaining  acceptance  and  is  being  suggested  by  many  of  the 
recent  publications  [8,19,20].   Its  use  was  reported  in  a 
research  paper  by  Fleet,  Stowers,  and  Swerdloff  [21]. 


Site  Analysis 

Data  obtained  from  the  0-D  survey  is  of  questionable 
applicability  in  predicting  nonresidential  trip  generation, 
especially  if  the  purpose  of  the  trip  estimation  is  system 
design  rather  than  area- wide  planning  [22].   While  steady 
growth  could  be  accounted  for  on  the  zone  level,  the  ernergence 
of  a  major  trip  generator  would  probably  best  be  treated 
separately.   Data  for  trip  generation  analysis  of  major 
generators  should  probably  be  collected  at  the  site  (genera- 
tor), rather  than  at  the  dwelling  units. 

In  the  base  year  analysis,  individual  site  investiga- 
tion can  supplement  the  zone-level  analysis.   In  the  con- 
tinuing phase,  periodic  checks  at  the  site  can  provide 
continuous  updating  of  the  trip  generation.   Site  analysis 
has  not  yet  been  used  extensively  in  operational  transportation 


24 


studies.   Significant  v/ork  has  been  reported  by  the  Cali- 
fornia Division  of  Highways  [23,24].   Special  traffic 
counts  are  done  at  selected  sites  for  a  period  ranging  from 
one  day  to  seven  days,  depending  on  the  situation.   For  the 
seven-day  or  five-day  (Monday  through  Friday)  period,  trip 
averages  per  24  hours  are  related  to  some  of  the  character- 
istics of  the  site  such  as  floor  area,  number  of  employees, 
number  of  patients,  or  student  enrollment  for  generators 
such  as  shopping  centers,  offices,  hospitals,  or  schools 
respectively  [25]. 

Both,  interview  techniques  and  volume  counts,  were  re- 
ported to  have  been  used  by  the  Chicago  Area  Transportation 
Study  to  analyze  trip  generation  at  a  high  rise  apartment 
building  and  at  O'Hare  International  Airport  [26]. 


Factors  in  Trip  Generation 
Urbanization  is  an  expression  of  the  i nterdependency 
of  human  activities.   The  concentration  of  people  in  urban 
areas  results  from  the  pursuit  of  opportunities  for  inter- 
action.  Face-to-face  interaction  and  person-to-place  con- 
tact presume  spatial  proximity.   Urban  travel  is  done,  in 
general,  for  these  purposes.   Excluding  any  possible  utility 
of  the  mere  act  of  travel,  the  demand  for  transportation  is 
a  derived  demand.   It  is  not  sought  for  its  own  sake,  but 
rather  to  satisfy  some  desire  for  interaction. 
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The  Trip  Generation  Mechanism 

For  the  proper  understanding  of  the  generation  of 
trips,  and  for  predicting  future  trip  generation,  travel 
should  be  observed  within  the  framework  of  the  total  urban 
picture.   Stated  otherwise,  trip  generation  is  better 
analyzed  and  described  by  the  same  scales  which  measure  the 
complicated  activities  of  the  urban  area. 

Travel  is  best  and  most  completely  described,  as  sug- 
gested by  Mitchell  and  Rapkin,  by  the  following  character- 
istics:  (1)  volume,  (2)  composition,  (3)  distance,  (4)  time 
rhythms,  (5)  location,  (6)  density,  (7)  kind  of  trio,  (8) 
kind  of  establishments  at  origin  and  destination,  and  (9) 
trip  maker  characteristics  [27].   Among  these  characteris- 
tics, the  last  three  are,  probably,  the  i^iost  relevant  to 
trip  generation. 

Trips  could  be  classified  according  to  purpose  such  as: 
work,  shopping,  school,  recreation,  or  personal  business. 
For  the  purpose  of  the  distribution  models  trips  are  usually 
classified  also  into  production  trips  and  attraction  trips, 
and  also  as  home-based  or  non  home-based  trips.   Person- 
trips  could  also  be  classified  according  to  the  mode  of 
travel  such  as  walking-trips,  auto-trips,  or  transit-trips. 
Most  0-D  studies  ignore  walking-trips  which,  of  course, 
introduces  a  bias  in  the  calibration  of  any  person-trips 
generation  model. 

In  establishing  a  method  for  identifying  independent 
variables,  three  methods  were  suggested  by  Mitchell  and 
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Rapki  n : 

1.  systems  based  on  kind  of  trip,  that  is,  trip 
purpose ; 

2.  systems  related  to  kind  of  establishment  or  area, 
that  is,  land  use,  either  by  parcel  or  by  area; 
and 

3.  systems  based  on  process  of  action  or  roles  [27]. 
A  comparable  set  of  methods  reported  by  Shuldiner  include: 

1  .   land  use  ; 

2.  activity -purpose;  and 

3.  others  [22]. 

Mitchell  and  Rapkin  suggest  that  travel  is  composed  of 
dispersive,  assembling,  and  random  components  [27].   The 
trips  generated  at  the  home  are  considered  dispersive,  they 
might  be  predicted  by  household  characteristics.   The 
destination  end  of  these  same  trips  could  be  predicted  by 
the  land  use  at  each  of  the  destinations.   In  the  prediction 
of  trips,  two  types  of  variables  are  available,  causative 
and  symptomatic  [28].   For  example,  the  school  enrollment 
in  a  zone  is  a  causative  variable  of  the  school  trip  attrac- 
tions.  Also,  the  number  of  single  family  dwellings  in  the 
same  zone  could  be  correlated  with  school-trip  atractions. 
The  number  of  single  family  dwellings  do  not  cause  the  school 
trip  attractions,  rather  it  is  a  symptomatic  variable  [28]. 
Since  most  of  the  characteristics  of  the  urban  area  are  re- 
lated, relationships  based  on  causative  variables  would  be 
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more  stable  with  time  than  those  based  on  symptomatic  vari- 
ables; thus  providing  a  better  predictive  model. 

Factors  Used  to  Estimate  Residential  Trip  Generation 

Residential  trips  are  of  the  dispersive  type.   They  are 
best  estimated  by  relating  then  to  the  characteristics  of 
the  househol d . 

Household  Size .   Oi  and  Shuldiner  stress  the  effect  of 
family  size  on  trip  generation  [9,  p.  85].   Stowers  and 
Kanwit  reported  a  higher  standardized  regression  coefficient 
(Beta  coefficient)  in  correlating  trips  to  family  size  than 
to  automobiles  owned  [20,  p.  50].   Williams  and  Robertson 
also  found  a  high  effect  of  family  size  on  trips  [29]. 

Car  Ownership .   Oi  and  Shuldiner  state  that  car  owner- 
ship is  the  one  variable  which  exhibits  the  closest  associa- 
tion with  reported  trip  generation  rates  [9,  p.  86].   They 
further  state  that  car  ownership  depends  on  other  variables 
such  as  real  income,  price,  and  financing  terms,  as  sup- 
ported by  studies  of  the  automobile  market  [9,  p.  87].   The 
impact  of  car  ownership  on  trip  rates  is  twofold.   It  first 
increases  the  proportion  of  the  residents  v;ho  are  trip  makers 
Second,  increased  car  ownership  enhances  the  intensity  of 
travel  [9]. 

Janes  reported  that  the  average  number  of  trips  per 
vehicle  for  two-vehicle  households  is  almost  one-third  less 
per  vehicle,  than  for  one-vehicle  households  [30,  p.  21]. 
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That  is  the  increase  in  vehicles  per  household  does  not 
proportionately  increase  the  trips  generated  by  these  house^ 
holds. 


Income .   Oi  and  Shuldiner  suggest  deleting  income  from 
the  analysis,  because  it  is  highly  correlated  to  many 
factors  that  are  more  highly  correlated  to  trips  [9,  p.  104], 
Stowers  and  Kanwit,  however,  rank  income  third  after  family 
size  and  car  ownership  in  the  order  of  significance  rela- 
tive to  trip  making  [20,  p.  50]. 

In  his  study  of  the  effects  of  car  ownership,  Wynn 
divided  the  population  into  four  income  classes  and  compared 
1948  to  1955  data  and  found  the  following: 

1.  the  number  of  cars  for  the  highest  three  income 
brackets  was  almost  constant,  and 

2.  between  1948  and  1955  the  cars  per  resident  was 
approaching  a  constant  for  all  income  classes. 

He  thus  concluded  that  "car  ownership  is  much  less  signifi- 
cant than  level  of  income."  [15,  p.  26]. 

Occupation  of  Head  of  Household,  and  Social  Status.   In 
a  study  by  Walker,  rate  of  trips  generated  by  a  household 
were  found  to  be  a  function  of  the  occupation  of  the  head 
of  household,  if  the  population  per  household  and  car  owner- 
ship were  controlled  [31].   Households  were  occupati onal  ly 
stratified  by  "white  collar"  and  "blue  collar"  employment, 
or  as  high,  medium,  and  low  occupational  groupings. 
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Earlier,  Oi  and  Shuldiner  tested  the  hypothesis  of  the 
significance  of  occupation  of  head  of  household  on  the  rate 
of  trips  generated  per  dwelling  unit  by  a  covariance  model. 
Although  the  test  was  significant,  they  questioned  the 
validity  of  the  linearity  assumptions  used  [9,  p.  116]. 

In  their  attempt  to  study  the  effect  of  social  status, 
Oi  and  Shuldiner  investigated  the  significance  of  the 
social  attributes  of  an  area  described  by  other  indices: 
(1)  social  rank,  (2)  degree  of  urbanization,  and  (3)  extent 
of  segregation.   These  social  area  indices  are  referred  to 
as  the  Shevky-Bell  Typology  [32]. 

Another  social  factor  proposed  by  Stowers  and  Kanwit 
is  the  "stage  in  the  family  life  cycle"  [20,  p.  50].   It 
ranked  fourth  after  family  size,  automobiles  owned,  and  in- 
come in  explaining  the  variability  in  household  travel 
habits. 


The  Distance  from  the  Central  Business  District  (CBD). 
It  has  been  suggested  by  Wynn  that  the  distance  from  the 
CBD  is  an  important  factor  in  trip  generation  rates  [15]. 
Some  researchers  have  argued,  however,  that  the  distance 
from  the  CBD  could  be  a  proxy  variable.   The  distance  of  a 
dwelling  unit  from  the  CBD  is,  in  fact,  very   much  related  to 
many  household  and  environmental  characteristics,  such  as 
family  size,  automobiles  owned,  income,  stage  in  the  family 
life  cycle,  and  residential  density  of  the  neighborhood. 
This  idea  was  illustrated  by  Mertz  et  al.  [33].   Oi  and 


30 


Shuldiner  have  shown  that  the  importance  of  this  factor  is 
insignificant  when  the  correlation  of  trip  generation  rates 
with  car  ownership  and  household  size  is  accounted  for  [9]. 
A  similar  conclusion  was  reached  by  Stowers  and  Kanwit  [20]. 

The  Type  of  Dwelling  Unit.   The  type  of  a  dwelling 
unit  reflects  the  style  of  life  of  its  occupants;  and  it 
is,  probably,  related  to  its  distance  from  the  CBD.   In 
studying  the  relationship  between  dwelling  unit  type  and 
rate  of  trip  generation,  Oi  and  Shuldiner  conclude  that  the 
effect  of  the  relationship  is  not  nearly  so  pronounced  as 
that  of  income  or  occupation  [9,  p.  121]. 

Additional  Remarks .   The  previous  discussion  on  the 
factors  affecting  residential  trip  generation  was  oriented 
toward  a  dwelling  unit  analysis.   When  the  analysis  is  at 
the  zone  level,  appropriate  equivalent  variables  are  usually 
used.   Population,  labor  force,  and  school-age  population 
are  some  examples.   Such  demographic  variables  are    corre- 
lated with  residential  trip  generation  at  the  zone  level. 
Population  was  noted  to  be  highly  correlated  with  residential 
trip  generation  by  Sharpe  et  al.  [34].   It  should  also  be 
noted  that  measures  of  the  intensity  of  the  use  of  land 
are  only  possible  to  calculate  in  areal  units  of  analysis, 
such  as  the  zone  or  district. 
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Factors  Used  to  Estimate  Nonresidential  Trip  Generation 
In  developing  models  for  estimating  nonresidential 
trip  generation  two  basic  approaches  are  in  use:   (1)  Land- 
Use  Based  Models,  and  (2)  Activity -Purpose  Models  [22,  p.  74] 
The  basis  for  the  above  classifications  is  to  facilitate 
the  review  of  each  model  and  does  not  necessarily  reflect 
basic  differences  between  them. 


Factors  Related  to  Land  Use .   A  basic  concept  in  trans- 
portation planning  is  the  recognition  of  the  functional 
relationship  between  travel  patterns  and  land  use,  which 
was  extensively  analyzed  by  Mitchell  and  Rapkin  [27]. 
Land-use  based  models  establish  relationships  between  the 
amount,  location,  and  intensity  of  use  of  different  classi- 
fications of  the  use  of  land  and  the  amount  of  trips  gener- 
ated by  them. 

In  providing  the  basis  for  such  a  procedure,  a  general- 
ized land  use  classification  system  is  needed.   The  same 
system  is  not  used  by  all  urban  area  transportation  studies 
Even  within  the  same  urban  area,  variations  in  trip  rates 
within  generalized  land  use  classes  have  been  far  too  large 
[22,  p.  74].   Another  source  of  difference  between  urban 
transportation  studies  is  introduced  by  the  methods  by 
which  trip  rates  are  calculated.   Trips  per  acre,  and  trips 
per  thousand  square  feet  of  floor  area  are  not  comparable 
because  of  the  range  in  establishment  sizes,  building 
heights,  and  land  coverage.   In  addition,  gross  acre  is 
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sometimes  used,  while  on  other  occasions  net  acre  is  used. 

The  use  of  the  floor  area,  as  a  basis  for  the  analysis, 
instead  of  the  acreage,  has  the  advantage  that  floor  area 
is  a  better  measure  of  the  intensity  of  use.   It  reflects 
the  effect  of  the  mult i -story  structures,  and  thus  the  char- 
acter of  the  area  also.   The  limited  use  of  floor  area 
measures  can  be  attributed  to  the  difficulty  and  expense 
of  obtaining  those  measures. 

Black  compared  land  area  versus  floor  area  as  a  pre- 
dictor of  nonresidential  trip  generation.   He  concluded: 
"Contrary  to  expectation,  floor  area  does  not  seem  consis- 
tently better  than  land  area.   Further  more,  floor  area 
trip  rates  are  not  uniform  throughout  a  metropolitan  area, 
but  increase  as  the  density  decreases."  [35,  p.  1]. 


Factors  Related  to  Land-Use  Activity .   This  approach 
recognizes  that  the  activity  which  takes  place  on  land  is, 
probably,  more  related  to  trip  generation  than  the  land  area 
in  different  uses.   Variables  such  as  employment  and  its 
various  subclass ificat ions,  retail  sales,  and  school  enroll- 
ment are  examples  of  activity  measures.   Land-use  activity 
measures  are  good  predictors  of  trip  generation  by  purpose. 

Black  compared  three  parameters  of  nonresidential  trip 
generation:   land  area,  floor  area,  and  employment.   He  used 
data  collected  by  CATS.   His  findings  indicate  that  no  one 
particular  measure  is  best  for  predicting  trip  generation 
for  all  land  use  categories.   He  concluded:   "Floor  area 
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seems  best  for  commercial,  employment  for  manufacturing, 
and  land  area  for  public  buildings,  public  open  space,  and 
transportation."  [35,  p.  1]. 

Accessibility  Considerations  in  Trip  Generation 
in  Previous  Investigations 

Little  has  been  done,  as  yet,  on  accounting  for  the 
effects  of  the  changes  in  the  transportation  netv;ork  on  the 
trip  generation  in  an  operational  transportation  study. 
Almost  in  all  the  cases,  the  network  consideration  were 
within  the  framework  of  the  competitive  models  of  trip  dis- 
tribution. 


Fort  Worth  Study 
The  Fort  Worth  Study  attempted  to  allocate  the  total 
number  of  trips  of  a  certain  purpose  which  were  derived 
from  estimates  of  household  trip  productions  from  all  zones 
This  was  used  in  lieu  of  estimating  the  trips  by  purpose 
attracted  to  each  zone.   Shuldiner  comments  that  this  ap- 
proach is  most  suited  for  estimating  highly  competitive 
trips,  such  as  shopping  trips  [22].   The  following  presenta- 
tion of  the  method  used  by  Fort  Worth  Study  relys  heavily 
on  Shuldiner's  article  [22].   A  set  of  attraction  indices, 
whi chare  weights  that  assign  to  each  land  use  factor  its 
relative  importance  in  the  trip  attraction  model  was 
developed.   Considering  the  home-based  other  trip  purpose 
category,  for  example: 
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Let 

Y.   =  the  basic  attractiveness  of  zone  i  for  home- 
based  other  trips  relative  to  all  other  zones. 
Basic  attractiveness  is  used  here  to  designate 
the  attractiveness  of  zone  i  without  regard  to 
its  location  or  accessibility  relative  to  all 
other  zones; 
P.   =  population  in  zone  i, 
C-   =  commercial  employment  in  zone  i, 
I.   =  industrial  employment  in  zone  i,  and 
0.  -    other  employment  in  zone  i. 
llith  the  variables  P.,  C.,  I.,  and  0.  is  associated  a  set 
of  attraction  indices:   Alp,  Alp,  AI,  and  AI„  respectively, 
which  are  weights  that  assign  to  each  variable  its  relative 
importance  in  the  trip  attraction  model.   Then,  Y.  = 
AIp(P.)  +  AI^(C.)  +  AIj(I.)  +  AIq(0.). 

"The  Y.'s  are  then  combined  with  appropriate  friction 
factors  to  distribute  the  total  number  of  ho me -based  other 
trips  among  the  various  zones  by  means  of  the  gravity  model. 
[22,  pp.  81-82].   This  is,  actually,  a  variation  of  the 
classical  gravity  model  of  distribution;  but,  instead  of 
using  trip  ends  as  the  measure  of  attractiveness  of  the 
zone  the  Y.'s,  as  defined  above,  were  used. 

Southeast  Connecticut  Area  Transportation  Study  (SEATS) 
SEATS  made  extensive  use  of  the  relative  attractiveness 
concept  [22].   Employment  was  used  as  the  single  measure  of 
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activity.   Nonresidential  land  was  classified  into  one  of 
the  eight  following  categories:   (1)  industrial,  (2)  personal 
service,  (3)  business  service,  (4)  institutional,  (5)  recre- 
ational, (6)  commercial  amusements,  (7)  retail,  and  (8) 
other.   Accordingly  industrial  employment  in  a  zone  is  the 
measure  of  industrial  activity,  institutional  employment  in 
the  same  zone  is  the  measure  of  institutional  activity,  and 
retail  employment  is  a  measure  of  the  retail  activity  in 
it,  etc.   This  measure  of  activity  was  called  the  destina- 
tion zone  factor.   The  trips  to  the  various  land  uses  were 
then  grouped  into  3  classes:   home-based  long,  home-based 
short,  and  non  home-based.   Trip  type  factors  for  each  class 
of  trips  to  each  land  use  activity  by  the  total  employment 
in  that  activity.   The  last  step  was  to  develop  an  attrac- 
tion index  for  each  zone  for  each  class  of  trip  as  a  func- 
tion of  the  destination  zone  factors  and  trip  type  factors. 

For  a  single  class  of  trips,  home-based  long,  the  pro- 
cedure in  notational  form  is  as  follows: 
Let 

k    =  a  subscript  denoting  a  particular  land-use 
activity  (k  =  1,  ...,  m);  and 

i     ■=    a  subscript  denoting  a  particular  zone 
( i  -    1  ,  .  .  .  ,  n  )  . 
Further  1 et 


L, 


Lk 


n 

Z   A 
i=l 


,  where 


ki 


36 


Lk 


'ki 


And  let 
I 


Li 


=  trip  type  factor  for  land  use  k, 

=  total  number  of  home-based  long  trips  to  land 

use  k  in  all  zones ,  and 
=  destination  zone  factor  for  land-use  activity 

k ,  in  zone  i  . 


^l^    hk   \i'  '^^^'^ 


Li 


attraction  index  for  home-based  long  trips  to 

zone  i  . 

The  attraction  index,  computed  above,  is  used  to  dis- 
tribute home-based  long  trips  among  the  various  zones.   The 
basic  assumption  of  the  model  is  that  the  trip  type  factor 
for  kth  land  use  F. .  is  the  same  for  all  zones. 

Baltimore  Metropolitan  Area  Transportation  Study  (BMATS) 
The  relative  attractiveness  of  each  traffic  zone  was 
related  directly  to  several  types  of  system  variables.   The 
study  report  stated  that  network  variables  which  indicate 
the  level  of  transportation  service  available  to  residents 
were  difficult  to  determine,  "since  in  past  studies  no 
standard  parameters  have  evolved."  [36,  p.  96].   The  highway 
accessibility  index,  as  a  measure  of  highway  service,  was 
defined  as  the  inverse  of  the  sum  of  the  travel  time  from 
each  zone  to  every  other  zone.   Time  was  measured  from  the 
minimum  time  path  trees.   The  highway  accessibility  index 
reflects  the  standard  of  service  as  measured  by  average  speed 
on  the  street  system  [36].   The  transit  service  index 
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represents  the  frequency  of  transit  service  available  to 
each  zone.   It  was  calculated  by  totaling  the  number  of  trans- 
it trips  accumulated  for  each  route.   These  two  indices  were 
used  with  other  variables  to  build  multiple  regression 
models  of  trip  generation. 

Indianapolis  Regional  Transportation 
and  Development  Study  (IRTADS) 

IRTADS  experimented  with  measures  of  accessibility  that 
were  adaptations  of  the  accessibility  indices  produced  by 
the  gravity  model.   Instead  of  using  trip  ends  as  a  measure 
of  attractiveness  for  zones,  as  normally  is  the  case  in  the 
gravity  model,  IRTADS  used  total  population  or  employment 
in  each  zone.   Using  travel  times  from  both  the  IRTADS 
highway  and  transit  networks,  the  gravity  model  program  was 
run  to  produce  both  highway  and  transit  accessibilities  to 
population  or  to  employment  [37].   Only  three  of  the  multiple 
regression  equations  of  trip  generation  included  the  accessi- 
bility variables,  and  each  with  a  negative  coefficient. 
Eventually  the  accessibility  indices  were  dropped  from  the 
trip  generation  analysis  [38]. 


Other  Transportation  Studies 
In  an  unpublished  paper  by  Corradino  on  the  trip 
generation  of  Sharon-Farrel  Study  Area,  the  set  of  indepen- 
dent variables  included  some  measures  of  the  highway  network 
variables  [39].   These  variables  included: 
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1.  reciprocal  of  time  to  center  of  -^jravity  of 
population; 

2.  reciprocal  of  time  to  center  of  gravity  of 
empl oyment ; 

3.  Log-,^  of  time  to  center  of  gravity  of  population; 
and 

4.  '-OSin  °^  time  to  center  of  gravity  of  employment. 
In  this  study,  Corradino  states  that  the  travel  time  vari- 
ables were  used  as  measures  of  decentralization  [39,  p.  7]; 
rather  than  as  measures  of  accessibility. 

In  a  tabular  summary  of  trip  generation  equations 
developed  by  several  transportation  studies  and  reported  in 
a  paper  by  Bellomo  and  Schultz  [40],  the  Boston  Study  equa- 
tions contained  four  accessibility  variables: 

1.  accessibility  by  highway  to  all  jobs; 

2.  accessibility  by  transit  to  all  jobs; 

3.  accessibility  by  highway  for  all  population;  and 

4.  accessibility  by  transit  for  all  population. 
Unfortunately  no  further  details  were  available. 
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CHAPTER  III.   DEVELOPMENT  OF  THE 
ACTIVITY-ACCESSIBILITY  CONCEPT 


The  Interactance  Hypothesis  and  Gravity  Concepts 
of  Human  Interaction 

In  very  general  terms,  the  gravity  concept  of  human 
interaction  hypothesizes  that  the  force  of  interaction 
between  two  human  activity  areas  is  directly  proportional 
to  some  function  of  their  population,  and  inversely  propor- 
tional to  a  measure  of  the  spatial  separation  between  them. 
The  hypothesis  relates  the  sizes  of  the  two  populations  to 
the  impedence  or  friction  induced  by  the  intervening  space. 
This  could  be  stated  formally 


I 


T  J 


f(P-  ,  P  •) 

g{D.  .) 


where 


I  .  . 
1  J 

P.  ,P  . 
1   J 

D.  . 
1  J 


interaction  between  area  i  and  area  j; 
the  population  of  areas  i  and  j,  repsectively; 
a  measure  of  the  separation  between  area  i  and 
area  j;  and 

f  and  g  stand  for  the  assumed  form  of  the  func- 
tional relationship  of  population  and  separation 
respectively. 
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Early  Formulations  of  the  Gravity  Concept 
The  earliest  known  explicit  formulation  of  the  gravity 
concept  of  human  interaction  was,  probably,  by  Carey  [41]. 
He  argued  that  the  same  fundamental  laws  govern  both  physi- 
cal objects  and  human  behavior.   Carey  stated: 

Man,  the  molecule  of  society,  is  the  subject  of 
Social  Science ... .The  great  Law  of  Molecule 
Gravitation  [is]  the  indispensable  condition  of 
the  existence  of  the  being  known  as  man.... The 
greater  the  number  collected  in  a  given  space, 
the  greater  is  the  attractive  force  that  is 
there  exerted  .... Gravi tati on  is  here,  as  every- 
where, in  the  direct  ratio  of  the  mass,  and  the 
inverse  one  of  distance  [42,  pp.  42-43]. 

This  concept  was  later  used  by  Ravenstein  and  Young  to 
explain  migration  phenomena.   Ravenstein  presented  empirical 
data  as  supporting  evidence  that  migratory  movement  tends 
to  be  toward  the  larger  cities,  and  that  the  volume  of 
movement  decreases  with  the  increase  of  distance  between  the 
origin  and  destination  [43,44].   Young  hypothesized  that  the 
number  of  migrants  drawn  from  several  sources,  to  a  single 
destination,  varies  directly  as  the  "force  of  attraction" 
of  the  destination  and  inversely  with  the  square  of  distance 
between  the  source  and  the  destination  [45]. 

Reilly's  "Law  of  Retail  Gravitation"  uses  the  gravity 
concept  to  explain  the  drawing  power  of  a  retail  center. 
In  his  formulation,  he  postulated,  that  a  retail  center 
will  attract  trade  from  its  surrounding  area  in  direct  pro- 
portion to  the  population  size  of  the  center  and  inversely 
to  the  square  of  the  distance  from  the  center  to  the 
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individuals  attracted  [46].   The  same  formulation  provides 
a  basis  for  determining  the  equilibrium  point  between  two 
retail  centers  competing  for  trade.   All  the  early  formula- 
tions, as  could  be  recognized,  were  partial;  in  the  sense 
that  the  size  of  only  one  of  the  interacting  groups  was 
consi  dered  . 


Structural  Formulation  of  the  Gravity  Concept 
In  the  1940's,  the  concept  of  human  interaction  was 
structured  by  Stewart  and  by  Zipf  following  the  gravita- 
tional analogy  [47,48].   The  formulation  by  Stewart  was  in 
terms  of  Newtonian  physics,  namely  that  the  "force"  of 
interaction  between  two  groups  is  directly  proportional  to 
the  product  of  the  population,  and  inversely  proportional  to 
the  square  of  the  distance  on  the  line  joining  them. 
Stewart  defined  the  "force"  mathematically  as  follows: 


T  J 


p.p. 
_] 1 

1  J 


where, 

F..   =  the  force  of  interaction  between  concentrations 

i  and  j  ; 

all  other  terms  as  previously  defined. 

Following  the  analogy  from  physics,  Stewart  defines 

the  "energy"  of  interaction  between  two  points  i  and  j  as: 

P.  P  . 
T J 


'U"    D, 
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where  k  =  constant  of  proportionality,  analogous  to  the 

gravitational  constant  of  physics. 
Stewart  extended  his  analogy  from  physics  to  include 
the  concept  of  population  "potential,"  which  measures  the 
intensity  of  the  possibility  of  interaction.   The  potential 
at  point  i,  or  the  intensity  of  the  possibility  of  interac- 
tion of  an  individual  at  i,  with  respect  to  a  group  at  point 
j  is 

P  . 


ij 


D.  . 


where 

v..  =  potential  at  i  of  the  population  at  area  j. 

The  total  potential  at  i,  that  is,  the  intensity  of  the 
possibilities  of  interaction  between  an  individual  at  i  and 
the  population  at  all  the  other  points  would  be: 
n   P  . 

J  =  l  1 J 

In  the  above  formulation,  it  is  important  to  include  the 
potential  of  the  population  of  a  certain  location  with 
respect  to  an  individual  at  the  same  location.   That  is,  it 
includes  the  term  P./D..,  where  D..  is  usually  taken  as  the 
average  distance  from  the  center  of  location  i  to  its 
boundaries  [41]. 


Variations  in  the  Distance  Function 
The  function  of  distance  in  the  potential  and  gravity 
formulations  has  been  a  source  of  controversy.   Experimental 
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results  reported  by  Price,  Ikle,  Carroll,  and  others  support 
the  argument  that  the  function  of  distance  in  the  basic 
formulations  of  population  potential  is  a  simple  inverse 
relationship  of  distance  raised  to  some  power  other  than 
one  [49,50,51].   Based  upon  empirical  studies,  the  distance 
exponent  ranged  from  one-half  to  over  three.   Anderson  has 
suggested  that  the  exponent,  itself,  is  a  variable,  inversely 
proportional  to  the  population  size  [52].   Consider,  for 
example,  two  centers  of  population  P.  and  P.  ,  were  P.  >    P,  , 

J         K  J      1^ 

both  centers  are  equidistant  from  a  third  center  P.,  that 
is  D..  =  D.,;  then  the  potential  of  the  smaller  center  P, 
will  be  more  reduced  by  the  distance  than  the  potential  of 


That  is  in  general 


P  . 


1  J 


D 


TJ 


where  , 

a 

Carrothers  argues,  however,  that  the  exponent  is  more 
likely  inversely  proportional  to  distance,  rather  than  to 
population  [41].   "That  is  an  extra  unit  of  distance  added 
to  a  long  movement  is  of  less  importance  than  an  extra 
unit  added  to  a  short  movement."  [41,  p.  97].   Thus  the  ex- 
pression of  the  exponent  of  distance  is: 
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Price  considered  in  addition  to  the  distance  the 
direction  of  the  migration  [49].   He  considered  migration 
movements  as  vectors  having  length  and  direction,  the 
length  corresponding  to  the  distance  of  migration,  and  the 
direction  depending  on  points  of  origin  and  destination. 


Variations  in  the  Population  Function 
In  studying  the  interaction  phenomenon  in  different 
regions,  Stewart  observed  that  the  effect  of  population  is 
not  always  the  same.   He  concluded  that  there  are  population 
characteristics  that  could  make  one  population  group  responsi 
ble  for  a  higher  interaction  potential  than  another.   To 
account  for  this  differential,  he  suggested  using  relative 
weights  to  account  for  the  population  characteristics.   Going 
back  to  the  analogy  from  physics,  this  is  in  fact  equivalent 
to  accounting  for  the  difference  in  the  molecular  weights  of 
two  particles.   Stewart,  thus,  suggests  the  following  formu- 
lation: 

k  m  .  P  .  •  m  .  P  . 
F    =     11     J   J 
1  J 


D.  . 
1J 


where 


m.  =  molecular  weight  of  an  individual  in  i;  and 
ni.  =  molecular  weight  of  an  individual  in  j. 
Dodd  in  his  formulation  of  the  interactance  hypothesis 
introduces  to  the  equation,  factors  other  than  population 
and  distance  [53].   These  factors  referred  to  as  Specific 
Indices  of  Influence  Level  should  account  for  different 
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characteristics  of  the  population,  such  as  sex, educati on ,  or 
income.   The  basic  energy  equation  thus  becomes: 


E.  .  =  k- 
1  J 


P.    (Is  .)    P  . 
1      J     J 


D.  . 


where 


Tr-  =  specific  indices  of  influence  level  for  popula- 
tion P  .  ,  and 
IS.  =  specific  indices  of  influence  level  for  popula- 

J 

t  i  0  n  P  . . 
Anderson  suggested  raising  the  numerator  of  the  basic 
equation  to  some  power  [52].   Carrothers  suggested  raising 
the  individual  population  elements  in  the  basic  equation 
to  a  power,  which  need  not  be  identical  for  both  populations 
[41].   When  all  the  variations  to  the  basic  formulation  of 
the  population  and  distance  functions  are  considered, 
Carrothers  suggests  the  following  general  expression: 


E.  .  =  k 
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a,  b,  and  c  are  empirically  derived  exponents. 


Variations  of  the  Basic  Formulation 
Empirical  investigations  have  suggested  different 
interpretations  in  drawing  analogies  to  man  and  distance. 
So  far  in  this  discussion,  population  has  been  used  as  the 
measure  of  mass.   In  studying  migration,  however,  variables 
such  as  employment,  or  the  income  of  the  regions  or  sub- 
areas  tends  to  be  a  more  significant  measure  of  mass  than 
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population  [54].   In  studying  market  potentials,  factors 
such  as  dollar  volume  of  sales  or  income  could  be  more 
significant  than  population.   The  choice  of  the  measure  of 
mass  is  controlled  by  the  nature  of  the  problem  at  hand, 
and  by  available  data.   The  wide  range  of  variables  which 
could  be  applicable  includes  employment,  school  enrollment, 
income,  telephone  calls,  newspaper  distribution,  etc. 

Similarly,  the  appropriate  measure  of  distance  need 
not  be  restricted  to  straight  line  distance.   Other  measures 
of  impedence  such  as  travel  time  and/or  travel  cost  has 
been  used  [55].   Some  problems  night  make  it  appropriate  to 
consider  road  distance,  fuel  consumption,  delay  time,  or 
even  intervening  opportunities  as  a  measure  of  "social 
distance"  [56]. 


The  Gravity  Concept  and  Urban  Trip  Distribution 
Since  the  mid  1950's,  there  has  been  an  increased 
interest  in  applying  the  gravity  concept  to  intra- 
metropolitan  movements.   Pioneer  work  by  Voorhees  advanced 
the  gravity  model  of  trip  distribution;  it  provided  a 
systematic  procedure  capable  of  synthesizing  zone-to-zone 
movements,  by  purpose,  of  alternative  configuration  of  land 
use  and  transportation  facilities  [55].   The  earliest  use 
of  the  gravity  model  in  an  operational  study  was  reported 
by  studies  at  Hartford,  Connecticut;  and  Baltimore,  Maryland 
[57].   Earlier  the  interactance  hypothesis  was  tested  with 
0-D  data  from  Charlotte,  North  Carolina;  and  St.  Louis, 
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Missouri;  by  Wynn  and  Linder  [58]. 

The  general  formulation  of  the  gravity  model  in  the 

framework  of  i ntra-metropol i tan  interchange  is: 
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where 

number  of  trips  produced  at  zone  i  and  attracted 
to  zone  j  ; 

total  trips  produced  by  zone  i; 
total  trips  attracted  by  zone  j; 
measurement  of  the  separation  between  zones  i 
and  j,  normally  expressed  in  terms  of  travel 
time;  and 
b    =  an  exponent,  to  be  determined  empirically. 
In  this  formulation  trip  productions  and  attractions 
are  the  measure  of  mass;  travel  time  is  the  measure  of 
distance,  and  it  is  raised  to  a  constant  exponent. 

Operationally,  and  for  reasons  v/hich  will  be  stated 
later,  the  following  formulation  is  used: 
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where 


T..,  P.,  A.   =  same  as  defined  above;  and 


F.  . 


a  set  of  "friction  factors,"  which 
measure  the  rate  of  trip  decay  as  a 
function  of  interzonal  travel  time. 
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The  replacement  of  the  inverse  function  of  distance  by  a 
set  of  friction  factors  introduces  computational  flexibility. 
It  does  not  presuppose  a  constant  exponent  because  the 
friction  factors  are  a  function  of  travel  time  and  trip 
purpose . 


The  Activity -Accessibility  Concept  in  Trip  Generation 
As  stated  earlier  in  CHAPTER  II,  the  rationale  of  trip 
generation  has  been  based  on  the  characteristics  of  the 
zone:   demographic,  economic,  and  land  use,  with  no  con- 
sideration for  the  status  of  the  transportation  system. 
This  research  is  based  on  the  concept  that  trips  generated 
by  a  zone  of  the  study  area  are  also  a  function  of  the  status 
of  the  transportation  subsystem  that  serves  a  zone  and  con- 
nects it  to  the  other  zones  of  the  study  area.   The  effect 
of  the  transportation  system  on  trip  generation  will  be 
investigated  in  the  light  of  the  relative  accessibility 
of  each  zone  to  various  urban  activities,  and  the  spatial 
relationship  of  the  different  zones  to  each  other.   The 
development  of  measures  of  relative  accessibility  will  be 
discussed  fully  later.   The  concern  at  this  stage  is  to 
point  out  instances  where  researchers  have  realized  the 
existence  of  a  feedback  from  the  transportation  system  on 
the  trip  generation  phenomenon.   To  date,  little  has  been 
done  to  measure  those  effects. 

Mitchell  and  Rapkin  in  their  discussion  of  the  systems 
of  urban  activities  state: 
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Martin,  Memmott,  and  Bone  in  their  discussion  of  the 

future  trip  generation  characteristics  give  the  following 

statements : 

...it  v/as  shown  that  the  transportation  studies 
assume  trip  generation  to  be  independent  of  the 
quality  of  transportation  facilities  and  the 
resulting  accessibility.   This  assumption  is 
obviously  incorrect,  as  induced  or  generated 
traffic  has  been  observed  to  be  an  important 
segment  of  the  total  traffic  on  new  facilities 
[59,  pp.  201-202]. 

Shuldiner,  in  investigating  nonresidential  trip  genera- 
tion, points  out: 

...as  far  as  prediction  of  future  trips  is 
concerned,  there  is  the  inevitability  of  change. 
All  present  [trip  generation]  rates  are  not 
fixed  and  immutable.   Rather,  they  are  based 
on  a  certain  structure  of  land  uses  and  land-use 
competition  which  is  constantly  changing.   Ag- 
glomeration, competition,  and  new  transportation 
all  affect  the  generation  of  a  given  attractor 
[22,  p.  75]. 

Domencich  et  al .  pointed  out  that  urban  transportation 

planning  has  developed  into  the,  almost  separate,  processes 

of  trip  production,  attraction,  distribution,  assignment, 

and  modal  split.   This  oversimplification  of  the  analysis  of 

a  complex  system  results  in  the  implicit  assumption:   "that 
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the  number  of  trips  generated  is  independent  of  the  per- 
formance of  the  transportation  system."  [2,  p.  66].   They 
add  that  this  assumption  implies  that  the  policies  imple- 
mented by  transportation  planners  have  no  effect  on  the 
total  number  of  trips  made.   Then  they,  conservatively, 
concl ude : 

While  it  is  possible  that  changes  in  the 
transportation  system  will  not  affect  trip 
generation,  there  is  no  good  reason  for  making 
this  assumption  a^  priori  [  2 ,  p .  66]. 

Walker,  in  an  attempt  to  develop  a  theory  of  trip 
attraction,  examined  how  land  was  used  and  for  what  purpose. 
He  pointed  out  that  the  activity  occurring  on  a  piece  of 
land  suggests  to  the  people  which  of  their  needs  could  be 
satisfied  by  interacting  with  the  particular  piece  of  land. 
He  also  observed  that  the  size  of  land  puts  a  limit  on  the 
volume  of  the  activity  that  could  be  carried  on  it.   Com- 
paring two  pieces  of  land  with  identical  size  and  function 
characteristics.  Walker  concludes:   "the  attraction  of... 
land  will  be  directly  related  to  the  land's  accessibility 
to  the  total  number  of  possible  interactions  perceived  by 
the  user."  [12,  p.  98]. 

Schneider  in  studying  the  relationship  of  trip  genera- 
tion to  access  started  by  a  modification  of  the  classical 
formulation  of  the  gravity  model  of  trip  distribution.   He 
then  imposed  the  conservation  law  of  equality  of  productions 
to  attractions,  under  the  simplifying  assumption  of  complete 
symmetry  between  origins  and  destinations.   This  could  be 
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illustrated  in  notational  form  as  follows 


R  .  F.  . 
T    .  p   2 U. 


where 


T..,  p.,  and  F..  are  the  same  as  defined  earlier  in  the 


"i 


formulation  of  the  gravity  model; 
=  an  undefined  quantity  measuring  that  which 

attracts  people  to  a  place,  rather  than  A. 

as  in  the  standard  formulation. 
Satisfying  the  conservation  law,  and  assuming  complete 
symmetry  between  origins  and  destinations,  we  get: 


P.  F.  .  R  . 
Z  F.  .  R  . 


P  .  F  .  .  R. 
J   Ji   1 
E  F  .  .  R. 
i  ^'       ' 


If  we  further  assume  the  friction  factors  are  symmetrical 
that  is,  F..  =  F..,  the  relation  reduces  to: 


P. 


P  . 


R.ZF..R.      R.ZF..R 


1  j   IJ   J 


J  ^-   Ji   1 


This  can  hold  true  for  each  pair  of  origin  and  destination 
zones  only  if  the  two  sides  of  the  equation  are  separately 
equal  to  the  same  constant  c,  which  yields  in  general: 


cZFR 


Stated  in  words  this  relationship  says  that  the  trip  density, 
at  a  point,  in  terms  of  origins  (or  destinations)  per  unit 
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of  R  is  equal  to  a  quantity  which  Schneider  called  the 

"access  integral  around  the  point."  [60,  p.  165]. 

Schneider  concludes: 

...the  number  of  trips  at  a  point  is  proportional 
to  the  accessibility  of  the  point  and  to  its 
attractiveness  to  people;  trip  ends  appear  at  a 
place  because  people  can  and  want  to  get  there. 
This  a  small  shift  from  the  customary  point  of 
view  in  which  trips  are  thought  of  as  occurring 
to  satisfy  the  craving  for  fulfillment  of  trip 
ends  [60,  p.  166]. 

The  cited  references  support  the  concept  of  this  re- 
search:  that  the  trip  generation  of  a  zone  is  not  only  a 
function  of  the  characteristics  of  the  zone  itself,  but 
also  is  affected  by  its  location  in  the  urban  area  relative 
to  other  zones  and  its  relative  accessibility  to  different 
land  use  activities.   It  is  necessary  now  to  make  the  trans- 
ition from  the  conceptual  to  the  operational  level. 

The  Conceptual  and  Operational  Definitions 
of  Access  i  bi 1 i  ty 

The  term  "accessibility,"  frequently  referred  to  in 
the  literature,  is  not  very  well  defined  neither  conceptually 
nor  operationally.   Different  authors  have  different  under- 
standing and  interpretation  of  it. 

Haig  understands  accessibility  as  "contact  with  rela- 
tively little  friction."  [61,  p.  38].   Mitchell  and  Rapkin 
comment  that  overcoming  space  friction  involves  two  types 
of  costs:   transportation  cost  and  site  rental  [27].   Con- 
necting the  two  costs  of  rental  and  of  transportation,  Haig 
maintains: 


53 


Rent  appears  as  a  charge  which  appears  that  the 
owner  of  a  relatively  accessible  site  can  impose 
because  of  the  saving  in  transportation  costs 
which  the  use  of  the  site  make  it  possible  [61, 
p.  38]. 

Hansen  defines  accessibility  as: 

...a  measurement  of  the  spatial  distribution 
of  activities  about  a  point,  adjusted  for  the 
ability  and  desire  of  people  or  firms  to  over- 
come spatial  separation  [62,  p.  73]. 

In  this  definition,  accessibility  measures  "...the  intensity 

of  the  possibility  of  interaction..."  [62,  p.  73].   As  such, 

it  is  a  generalization  of  the  "population  potential"  concept 

developed  by  Stewart  [47]. 

Martin  et  al.  define  the  accessibility  of  a  site  as: 

...the  weighted  summed  travel  resistance  be- 
tween the  particular  site  and  all  other  sites 
on  which  there  are  mutually  interacting  ac- 
tivities [59,  p.  76]. 

The  operational  definitions  of  accessibility  are  even 
more  diversified  than  the  conceptual  definitions.   The 
simplest  operational  definition  is  to  measure  accessibility 
of  a  site  to  another  by  the  inverse  of  some  function  of  the 
distance  or  travel  time  between  the  two  sites.   This  defini- 
tion is  completely  devoid  of  any  measurement  of  the  char- 
acteristics of  the  two  sites  of  concern.   It  reflects,  only, 
the  degree  of  spatial  separation  between  them. 

Shindler  and  Ferrer  1  calculated  the  "highway  accessi- 
bility" and  the  "transit  accessibility"  of  a  residential 
zone  to  a  nonresidential  zone  by  multiplying  the  reciprocal 
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of  the  travel  time  squared,  between  the  two  zones,  by  the 
respective  total  trips  attracted  to  the  nonresidential  zone 
by  highway  or  transit.   The  reciprocal  of  the  travel  time 
squared  was  referred  to  as  the  "friction  factor."    The 
"friction  factor,"  according  to  Shindler  and  Ferrer i, 
attempts  to  measure  the  restraining  influence  which  time 
has  on  travel  between  zones  [63].   The  implied  assumption  is 
that  the  restraining  influence  on  travel  varies  as  the 
square  of  the  travel  time. 

Savigear  attempted  to  measure  accessibility  quantita- 
tively as  a  function  of  travel  time  and  availability  of 
parking.   Computationally,  accessibility  is  derived  as  fol- 
1  ows  : 


Let 


measure  of  demand  on  trips  from  zone  i  to  the 

zone  under  consideration,  and 

travel  time  from  zone  i  to  the  zone  under 

consideration. 


Further  let 


t  - 


Zn.t. 


t  is  a  weighted  mean  average  travel  time  of  trips  to  the 
zone  under  consideration  from  all  other  zones.   An  initial 
estimate  of  the  accessibility.  A,  of  the  zone  under  con- 
sideration is:   A  =  —  [64].   Savigear  adds  that  any  measure 

t 
of  highway  accessibility  should  take  parking  into  account. 
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Let 

f   =  probability  of  being  able  to  park  at  the  zone  of 
destination,  for  as  long  as  desired. 
Then , 

A',  a  better  indicator  of  accessibility  is  A'  =  f.A 
[64]. 
Savi gear's  definition,  as  it  could  be  inferred,  is  only  in 
terms  of  travel  time  weighted  by  the  demand  for  travel  on 
the  system,  and  the  availability  of  parking  at  the  destina- 
tion.  No  considerations  are  given  to  the  activities  at 
di  f f erent  zones  . 

A  similar  operational  definition  of  accessibility  was 
used  by  the  Baltimore  Metropolitan  Area  Transportation  Study 
(BMATS).  As  a  measure  of  the  highway  service,  the  inverse 
of  the  sum  of  the  travel  times  from  each  zone  to  every  other 
zone  was  used  and  termed  the  "accessibility  index"  [36,  p. 
96].  The  travel  times  were  measured  from  the  minimum  path 
trees  . 

A  more  comprehensive  operational  definition  of  accessi- 
bility was  suggested  by  Hansen.   It  is  the  ratio  of  "popu- 
lation over  distance"  [62,  pp.  73-74]: 
S 


A.  . 


D^ 


a 
1  J 


where 


A..   ==  relative  accessibility  of  zone  i  to  an  activity 
1  ocated  i  n  zone  j ; 
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S.    =  size  of  the  activity  in  zone  j  (a  function  of 
population,  employment,  land  use,  etc.); 

D..   =  travel  time  or  distance  between  zone  i  and  zone 
j  ;  and 

a    =  time  or  distance  exponent. 
Considering  more  than  two  zones: 


A. 
1 


n   S 
j  =  l   D 


a 
i  J 


where 


n    =  total  number  of  zones; 

A.^   =  relative  accessibility  of  zone  i  to  an  activity 
in  all  zones. 

Swerdloff  and  Stowers  suggested  that  instead  of  using 
the  travel  time  or  distance  raised  to  an  exponent,  to  use 
the  alternate  formulation: 


A. 
1 


J-1 


S  .   F.  . 
J    1J 


where 


F..   =  the  friction  factor,  corresponding  to  a  D.. 

separation  as  derived  from  a  calibrated  gravity 
model  of  the  area  [65]. 
The  hypothesis  proposed  by  this  research  is  that  the 
number  of  trips  generated  by  a  zone  are  a  function  of  the 
transportation  subsystem  that  connects  the  zone  under  con- 
sideration with  the  other  zones  of  the  study  area.   For 
this  purpose,  the  term  "relative  accessibility"  will  be 
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defined  conceptually  and  operationally. 

It  was  hypothesized,  earlier,  that  the  trips  produced 
by  or  attracted  to  a  zone  are  a  function  of  the  causal  or 
symptomatic  variables  modified  by  the  relative  ease  in  over- 
coming space  between  that  zone  and  all  other  zones.   Zones 
with  "relatively  more  accessible"  destinations  should,  in 
general,  produce  more  trips;  similarly,  zones  that  are 
"relatively  more  accessible"  to  origins  should,  in  general, 
attract  more  trips.   The  term  "relatively"  refers  to  the 
zone  under  consideration  as  compared  to  all  other  zones 
of  the  study  area.   This  implies  a  competitive  consideration 
among  zones  in  generating  trips.   Zones  of  similar  size 
activities  will  attract  trips  differently  according  to 
their  locational  and  accessibility  advantages. 

As  a  measure  of  the  ease  or  difficulty  of  overcoming 
space,  this  study  intends  to  use  the  set  of  "friction 
factors"  developed  from  the  calibrated  gravity  model  of 
trip  interchange  for  the  study  area.   This  has  the  advantage 
of  avoiding  the  use  of  a  constant  exponent  of  distance,  or 
time.   Friction  factors  as  developed  by  the  calibration  of 
the  gravity  model  are  a  function  of  travel  time  and  are 
classified  by  trip  purpose. 

This  study's  definition  of  "relative  accessibility"  is 
a  modification  of  Hansen's  [62].   In  notational  form, 
"relative  accessibility"  is  computed  as  follows: 
let 

i      =  zone  under  consideration,  (i  -    1,  2,  ...,  n); 
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jk 


^•j(£) 


A. 


i  •!<(£) 


Then  , 


A. 


i-k(£) 


any  zone  in  the  study  area,  including  zone 

i  ,  ( J  =  1  ,  2  ,  .  .  .  ,  i  ,  .  .  .  ,  n )  ; 

activity  under  consideration,  (k  =  1,  2,  ...,  m); 

trip  purpose,  {l    =  1,  2,  ...,  p); 

size  of  activity  k  in  zone  j ; 

friction  factor  corresponding  to  the  travel 

time  from  zone  i  to  zone  j  for  purpose  Z; 

and 

accessibility  of  zone  i  to  activity  k  for 

purpose  £. 


/.T  ^Jk  ^•j(£)'  ''' 


RA .  .  /  X  =  relative  accessibility  of  zone  i  to  activity 
kforpurpose£. 
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The  value  of  the  relative  accessibility  of  zone  i  to 
activity  k  for  purpose  £,  as  it  is  possible  to  infer  from 
its  formulation,  could  be  different  for  a  future  year  if  any 
or  all  of  the  following  study  area  parameters  change: 

1.   interzonal  travel  time,  consisting  of  interzonal 
driving  time,  terminal  time,  and  intrazonal  time, 
or 
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2.   size  of  activity  k,  in  any  or  all  zones  of  the 
study  area  . 

The  value  of  the  relative  accessibility  of  a  zone  to 
the  sane  activity  could  be  different  for  different  trip 
purposes.   The  reason  being  that  the  friction  factors, 
corresponding  to  a  certain  travel  time,  are  usually  differ- 
ent for  different  trip  purposes. 

Another  aspect  of  accessibility,  neither  directly  re- 
lated to  travel  time  nor  to  the  size  of  activity,  is  con- 
sidered.  The  relationship  of  location  to  trip  generation 
characteristics  is  investigated  by  stratifying  the  zones  of 
the  study  area  into  two  disjoint  sets:   central  and  non 
central.   The  conceptual  basis  being  that  "central  sites 
afford  maximum  accessibility..."  [27,  p.  108].   The  central 
area  is  also,  more  or  less,  equiaccessible  to  the  various 
zones, of  the  study  area,  because  of  the  convergence  of  the 
street  system  on  the  city  center.   This  stratification 
introduces  a  qualitative  factor  describing  the  general 
arrangement  of  the  land  uses  and  the  configuration  of  the 
street  system. 
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CHAPTER  IV.   DEVELOPMENT  OF  ACTIVITY-ACCESSIBILITY 
TRIP  GENERATION  MODELS 


The  following  methodology  was  used  in  this  investiga- 
tion.  Trip  generation  models  which  take  into  account  ac- 
cessibility variables  were  developed  from  data  from  an  opera- 
tional transportation  study.   These  models  were  then  compared 
with  the  conventional  models  developed  as  part  of  the  trans- 
portation study.   Both  sets  of  models  were  used  to  forecast 
1985  trip  generation.   The  two  sets  of  forecasts  were  com- 
pared by  testing  for  any  significant  differences,  on  a  zone- 
by-zone  basis,  between  the  two  forecasts. 

Models  that  take  into  consideration  the  stratification 
of  the  zones  of  the  study  area  into  central  and  non  central 
sets  were  also  developed.  This  stratification  was  investi- 
gated for  the  models  developed  by  the  transportation  study 
and  those  developed  by  this  investigation.  Also  investigated 
were  the  ranges  of  the  independent  variables  for  the  survey 
year  and  the  forecast  year. 

Since  the  main  purpose  of  this  investigation  was  to 
compare  the  sets  of  developed  models,  with  accessibility 
variables  and  zones  stratified  into  central  and  non  central 
sets,  to  models  developed  in  the  traditional  approach  by  an 
operational  transportation  study,  care  was  taken  to  keep  any 
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factors  that  might  disturb  the  comparison  out  of  the  developed 
models;  so  that  the  comparisons  would  be  most  valid. 

The  decision  to  develop  multiple  linear  regression 
models  of  trip  generation  using  data  summarized  by  zone  was 
mainly  in  the  interest  of  keeping  the  results  of  this  in- 
vestigation comparable  to  those  from  an  operational  trans- 
portation study. 

Data  Preparation 
The  data  used  in  this  investigation  werethe  results 
of  the  surveys  conducted  for  the  Indianapolis  Regional 
Transportation  and  Development  Study  (IRTADS).   This  deci- 
sion was  based  mainly  on  the  availability  of  data  and  the 
fact  that  the  IRTADS  survey  is  recent.   The  study  was  at 
a  stage  where  most,  if  not  all,  of  the  analyses  were  com- 
pleted, the  forecasts  established,  and  the  proposed  networks 
eval ua ted . 


Ready-Available  Data 
The  trip  generation  analysis  for  IRTADS  was  completed 
in  October,  1966;  and  was  reported  in  a  Technical  l\'ork 
Paper  [38].   All  the  models  developed  by  IRTADS  were  for 
total  person-trips  (except  for  truck  and  taxi  equations). 
The  dependent  variables  were  in  the  form  of  productions  and 
attractions  suitable  for  distribution  by  the  gravity  model. 
Nineteen  trip  generation  equations  were  developed  by  IRTADS. 
Two  of  which  were  for  truck  and  taxi  trip  ends,  another  four 
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were  for  control  totals,  and  the  rest  for  person-trips  pro- 
ductions and  attractions  by  purpose  [38]. 


The  Dependent  Variables .   This  investigation  was 
limited  to  six  trip  production  purposes,  five  trip  attrac- 
tion purposes,  and  two  control  totals  (one  for  all  produc- 
tions and  one  for  all  attractions). 

Trip  generation  equations  were  developed  for  the 
following  dependent  variables: 

1.  Home-based  work  person-trip  productions. 

2.  Ho me -based  shop  person-trip  productions. 

3.  Home-based  school  person-trip  productions. 

4.  Home-based  other  person-trip  productions. 

5.  Non  home-based  work-oriented  person-trip 
productions. 

6.  Non  home -based  non  work-oriented  person-trip 
productions. 

7.  Total  person-trip  productions. 

8.  Home-based  work  person-trip  attractions. 

9.  Home-based  shop  person-trip  attractions. 

10.  Home-based  other  person-trip  attractions. 

11.  Non  home-based  work-oriented  person-trip 
attractions. 

12.  Non  home-based  non  work-oriented  person-trip 
attractions . 

13.  Total  person-trip  attractions. 
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It  was  not  possible  to  develop  an  equation  for  home- 
based  school  person-trip  attractions  because  the  key  inde^ 
pendent  variable,  school  enrollment,  was  not  available. 


Socio-Economic  and  Land  Use  Variables.   A  total  of  29 
socio-economic  and  land  use  variables  were  originally  con- 
sidered by  IRTADS  in  their  trip  generation  analysis;  only 
15  were,  however,  eventually  retained  in  the  final  equa- 
tions [38].   Moreover,  only  10  of  these  were  available  for 
this  study  for  both  the  survey  year,  1964,  and  the  forecast 
year,  1985.   The  socio-economic  and  land  use  variables,  that 
were  available  and  used  by  this  study  as  independent  vari- 
ables in  developing  multiple  linear  regression  models  of 
trip  generation,  were: 

1 .  Total  empl oyment . 

2.  Retail  employment. 

3.  Service  employment. 

4.  Retail  floor  area. 

5.  Educational  floor  area. 

6 .  Dwelling  units. 

7.  Labor  force. 

8.  Popul ati  on . 

9.  Cars. 

10.   Single-famil}'  dwelling  units 
Appendix  A  includes  an  abbreviation  key  for  all  the  vari- 
ables used  in  this  investigation. 
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Generating  Accessibility  Variables 
Different  measures  of  relative  accessibility,  to  be 
used  as  independent  variables  in  trip  generation  regression 
models,  were  established.   The  operational  definition  of 
relative  accessibility,  presented  earlier  in  CHAPTER  III, 
isrestatedhere: 
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=    relative  accessibility  of  zone  i  to 

activity  k  for  trip  purpose  £; 
-    size  of  activity  k  in  zone  j; 
=  a  friction  factor  corresponding  to  the 
total  travel  time  from  zone  i  to  zone  j 
for  purpose  £ ; 
i  =  zone  under  consideration,  (i  =  1,  2,  .. 

j  =  any  zone  in  the  study  area,  including  zone 

i  ,  ( j  =  1  ,  2  ,  . . .  ,  i ,  . . . ,  n ) ;  and 
£  =  trip  purpose,  (£  =  1,  2,  ...,  p). 

Thus  the  required  inputs  are  the  size  of  various  ac- 
tivities in  each  zone,  and  skim  zone-to-zone  friction  factor 
trees  . 

Size  of  Activities.   The  definition  of  activity  was 
extended  for  this  purpose  to  include  all  the  ten  socio- 
economic independent  variables  that  were  available  from  the 
IRTADS  surveys.   If  the  information  had  not  been  available 


,  n)  ; 
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from  the  above  source,  the  size  of  activities  would  have 
been  collected  from  various  sources,  such  as:   employment 
data,  from  the  State  Employment  Securities  Division;  floor 
area  information,  from  the  land  use  survey;  and  others,  from 
the  home-interview  survey  or  census  records. 


Skim  Zone-to-Zone  Friction  Factor  Trees.   A  set  of 
friction  factors  for  each  of  six  trip  purposes  were  avail- 
able from  the  results  of  the  IRTADS  calibrated  gravity 
model  [66].   Table  Bl  shows  the  values  of  the  friction  fac- 
tors for  travel  times  of  one  to  forty  minutes  for  each  of 
the  six  trip  purposes.   The  procedure  of  calibrating  a  grav- 
ity model  is  well  documented,  and  computer  programs  are 
available  for  that  purpose  [67].  ' 

A  binary  zone-to-zone  trees  tape  was  made  available 
for  this  investigation  by  IRTADS.   This  binary  tape  was  up- 
dated (intrazonal  and  terminal  times  added)  and  skimmed  to 
give  a  skim  zone-to-zone  travel  time  binary  trees  using  the 
program  PR-130  [68,  pp.  A-59  -  A-64].   Six  skim  zone-to-zone 
friction  factor  trees  were  built,  one  for  each  of  six  trip 
purposes.   This  was  achieved  using  the  program  MATCH, 
written  for  the  IBM  7094.   Skim  zone-to-zone  travel  time 
binary  trees  tape  and  the  six  sets  of  friction  factors  were 
inputs  to  MATCH.   The  outputs  were  on  tape  in  BCD  (binary 
coded  decimal)  format,  which  is  FORTRAN  readable;  sample 
printouts  were  also  obtained  to  check  results.   The  program 
MATCH  is  listed  in  Appendix  B. 
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It  should  be  noted  that  only  the  highway  network  was 
considered  in  developing  relative  accessibility  variables. 
Although  the  trip  generation  models  were  for  person-trips, 
using  the  highway  network  only  would  not  introduce  any 
appreciable  bias  in  the  case  of  IRTADS;  mainly  because 
transit  passenger  trips  constituted  only  4.1  percent  of  all 
the  person-trips  [70,  p.  25].   Moreover,  the  transit  in 
IRTADS  area  was  entirely  bus  service  on  the  city  streets. 


The  Program  "ACCESS .  "   A  computer  program  "ACCESS"  was 
written  for  either  the  IBM  7094  or  the  CDC  6500.   The  pro- 
gram accepts  friction  factor  tapes  and  the  activity  size 
variables  as  inputs  and  generates  relative  accessibility 
measures  for  each  zone  in  the  study  area.   A  listing  of 
ACCESS  is  presented  in  Appendix  B.   The  sequence  of  computer 
programs  and  the  data  flow  for  calculating  relative  accessi- 
bilities are  shown  schematically  in  Figure  5. 

With  ten  activity  measures  (the  available  independent 
variables)  and  six  sets  of  friction  factors  (one  for  each 
trip  purpose),  sixty  measures  of  relative  accessibility 
could  be  generated.   However,  not  all  sixty  possible  combina- 
tions were  generated;  only  those  that  were  meaningful  to 
the  trip  generation  analysis  were  used.   For  example,  the 
relative  accessibility  of  a  zone  to  retail  floor  area  could 
be  meaningful  in  conjunction  with  home-based  shopping  person- 
trip  productions.   Also,  the  relative  accessibility  of  a 
zone  to  single  family  dwellings  would  be  meaningful  in 
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conjunction  with  home -based  shop  person-trip  attractions. 
Twenty  relative  accessibility  variables  were  generated  and 
considered  in  the  analysis;  these  are  presented  in  Figure 
6.   It  should  be  noted  that  the  same  measure  of  relative 
accessibility  could  be  meaningful  in  conjunction  with  both 
the  productions  and  attractions  of  some  of  the  trip  purposes 

A  notation  defining  an  accessibility  variable  should 
indicate  the  activity  under  consideration  and  the  trip 
purpose  whose  set  of  friction  factors  were  used  in  calcu- 
lating that  specific  accessibility.   A  relative  accessibil- 
ity variable  was  abbreviated  by  a  leading  letter  A  denoting 
accessibility,  followed  by  the  abbreviation  of  the  activity 
under  consideration,  and  ended  by  a  one  digit  number.   The 
number,  one  to  six,  indicated  the  trip  purpose  of  the 
friction  factor  used:   home-based  work,  home-based  shop, 
home-based  school,  home-based  other,  non  home-based  work- 
oriented,  or  non  home-based  non  work-oriented  respectively. 


Delimiting  the  Central  Area 
As  indicated  in  CHAPTER  III,  this  investigation  con- 
sidered stratifying  the  study  area  into  central  and  non 
central  areas.   It  was  assumed  that  the  central  and  non 
central  areas  might  reflect  two  different  trip  generation 
patterns  due  to  the  shape  of  the  study  area,  its  historical 
quasi -annular  urban  growth,  and  the  configuration  of  the 
transportation  system.   This  differentiation  of  the  central 
and  non  central  areas  was  categorical  instead  of  numerical, 
and  could  best  be  treated  through  stratification. 
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MEASURES  OF  ACTIVITY 

TRIP      PURPOSE 

HBWK 

1 

HBSHP 
2 

HBSCL 
3 

HBOTR 
4 

NHBWK 
5 

NHBNWK 
6 

Total     Employment 

P 

P  -  A 

P  -   A 

Retai  1     Employment 

Service     Employment 

Retail    Floor    Area 

P 

P   -   A 

P  -   A 

P  -  A 

Educational    Floor    Area 

P 

P    -   A 

P   -  A 

P  -  A 

Dwelling    Units 

A 

P   -   A 

Labor     Force 

A 

P   -    A 

Populat  ion 

A 

P  -   A 

Cars 

P  -   A 

Single      Family    Dwellings 

A 

P   -   A 

*   P  (and/or   A)    in   cells   indicates    that    the  corresponding 
relative     accessibility  variable    was   considered   in 
developing  models  of  trip   productions  (and/or  attractions) 
for    the    indicated     trip  purpose 


FIGURE  6- THE    GENERATED     RELATIVE     ACCESSIBILITY 
VARIABLES 
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Area!  Definitions  of  the  Study  Area.   The  IRTADS  Study 
Area  is  bounded  by  a  cordon  encompassing  all  of  Marion 
County  and  portions  of  Johnson  and  Hamilton  Counties.   The 
area  within  the  cordon  is  divided  into  nine  major  subareas 
termed  sectors  (sectors  zero  through  eight).   Sector  zero 
is  the  Central  Business  District  (CBD);  "...defined  as  the 
area  within  the  inner  loop."  [69,  p.  21].   The  remaining 
eight  sectors  radiate  outward  from  the  CBD,  dividing  the 
Study  Area  into  eight  pie-shaped  wedges. 

The  nine  sectors  are  further  divided  into  72  districts. 
The  districts  are  divided  into  395  zones  and  the  zones  are 
further  divided  into  1200  subzones.   Host  of  the  analysis 
by  IRTADS  and  by  this  investigation  were  carried  on  at  the 
zone  level.   The  delimitation  of  the  central  area,  however, 
was  carried  at  the  district  level;  because  the  suggested 
procedure  was  not  thought  to  be  refined  and  precise  enough 
to  warrant  an  analysis  on  the  zonal  level.   The  general 
layout  of  the  Study  Area,  and  its  sectors  and  districts  are 
shown  in  Figure  7. 


Criteria  for  Delimiting  the  Central  Area .   The  rationale 
behind  the  procedure  developed  to  delimit  the  central  area 
was  tied  to  the  expected  character  and  attributes  of  a 
central  area.   A  larger  oroportion  of  the  land  in  a  central 
area  is  expected  to  be  in  urban  use.   A  relatively  small 
proportion  of  the  land  in  a  central  area  is  expected  to  be 
devoted  for  residential  uses.   Conversely,  a  relatively  high 
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FIGURE   7  -    STUDY    AREA      SECTORS    AND     DISTRICTS 
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proportion  of  the  land  in  the  central  area  was  expected  to 
be  in  uses  which  are  known  to  seek  central  location. 

Land  was  considered  to  be  in  urban  use  if  it  did  not 
belong  to  any  of  the  following  classifications: 

1  .   quarry  and  mi  ni  ng  i 

2.  automobile  junkyard; 

3 .  water; 

4.  in  agricultural  use;  or 

5 .  vacant. 

The  percentage  of  urban  land  in  residential  use  was  the 
measure  of  the  intensity  of  residential  activity.   The  per- 
centage of  land  in  urban  use  was  the  measure  of  urbanity. 
Among  the  different  trade  and  service  uses  on  which  IRTADS 
had  floor  area  information,  the  following  were  chosen  as 
uses  that  seek  central  location: 

1.  wholesale    trade    (i'/ithout    warehousing); 

2.  general    retail    trade; 

3.  auto    retail; 

4.  apparel,  furniture,  and  appliance  retail; 

5.  retail  use  not  otherwise  classified; 

6.  finance,  business,  and  professional  services; 

7.  contract  construction  services; 

8.  governmental  services; 

9.  personal  services;  and 

10.   services  not  otherwise  classified  [71,  pp.  42-50]< 
Educational  services  were  excluded  because  schools  do  not 
necessarily  seek  central  locations.   The  floor  area  in 
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hundreds  of  square  feet  consumed  by  the  above  10  uses  in 
each  district,  per  acre  of  land  in  urban  use  in  the  same 
district  was  calculated.   This  ratio  measures  not  only  the 
amount  of  these  uses  which  are  represented,  but  also  re- 
flects the  intensity  of  use.   If  this  ratio  was  unitless, 
that  is,  the  numerator  and  denominator  were  of  the  same 
units,  it  would  be  indicative  of  the  average  number  of  floors 
per  structure  in  each  district. 

Using  the  above  three  measures,  the  following  condi- 
tions were  set  for  delimiting  the  central  area: 

1.  The  delimitation  was  to  be  performed  at  the 
district  level. 

2.  The  central  area  would,  probably,  include  all  of 
sector  zero  and  some  of  the  qualifying  surrounding 
districts. 

3.  The  districts  of  the  central  area  should  all  be 
contiguous  and  connected. 

4.  A  district  which  would  qualify  must  satisfy  at 
least  two  of  the  three  following  criteria: 

a.  In  the  lower  quart ile  of  all  the  districts  of 
the  study  area  in  percentage  of  urban  land  in 
residential  use. 

b.  In  the  upper  quartile  of  all  the  districts  of 
the  study  area  in  the  percentage  of  land  in 
urban  use. 
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c.   In  the  upper  quartile  of  all  districts  of  the 
study  area  in  the  ratio  of  hundreds  of  square 
feet  of  uses  usually  seeking  central  location 
to  acres  of  urban  land  in  each  district. 
For  each  district,  the  values  of  the  variables  defining 
the  three  criteria  were  calculated.   The  values  correspond- 
ing to  criterion  "a"  v/ere  sorted  in  ascending  order;  the 
values  for  the  other  two  criteria  were  sorted  des cendi ngly . 
Since  the  total  number  of  districts  is  72,  then  the  first 
18  districts  corresponding  to  the  sorted  values  would 
represent  the  lower  quartile  under  criterion  "a,"  and  the 
upper  quartiles  under  criterion  "b"  and  criterion  "c."   To 
avoid  any  illusive  precision,  the  values  of  each  criterion 
corresponding  to  the  18th  ordered  district  were  rounded  to 
the  nearest  integer.   Thus,  numerically,  the  cut-off  points 
for  the  three  criteria  were  determined  to  be: 

Criterion  a:   percent  of  urban  land  in  residential  use 

less  than,  or  equal  to  32%. 
Criterion  b:   percent  of  land  in  urban  use  greater 

than ,  or  equal  to  88% . 
Criterion  c:   hundreds  of  square  feet  of  selected 

trade  and  services  per  acre  of  urban  land 
greater  than,  or  equal  to  8. 
The  supporting  computations  are  presented  in  Table  CI  . 

All  the  nine  districts  of  sector  zero  and  districts: 
11,  21,  31,  41,  51,  52,  61,  71,  and  72  satisfied  the 
four   conditions.     District   81   satisfied   only 
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criterion  "a"  of  condition  4.   flowever,  it  was  decided  to 
include  it,  in  order  to  keep  the  central  area  completely 
connected . 

A  pictorial  representation  of  the  delimitation  process 
is  presented  in  Figure  CI  of  Appendix  C.   The  study  area 
as  stratified  into  central  and  non  central  areas  is  shown 
in  Figure  8.   The  districts  of  the  central  area  constituted 
105  zones  out  of  the  395  in  the  study  area. 

Model  Building 

Guidelines  for  Model  Building 
Multiple  linear  regression  models  of  trip  generation 
were  developed  using  the  computer  program  "BMD-2R,  stepwise 
regression"  [72].   In  addition  to  the  desired  statistical 
qualities  of  the  developed  models,  other  important  factors 
were,  also,  considered. 


Conceptual  Valid ity .   The  consideration  of  relative 
accessibility  was  mainly  to  achieve  a  sounder  conceptual 
basis  for  trip  generation.   In  addition,  only  independent 
variables  that  were  logically  related  to  the  specific  de- 
pendent variable  under  consideration  were  allowed  to  enter 
when  developing  regression  equations  for  that  dependent 
variable.   The  causal-logical  relationship  was  considered 
prior  to  the  mere  statistical  correlation  analysis.   As- 
sociation and  correlation  do  not  prove  causality;  causality 
should  only  be  hypothesized  on  theoretical  or  conceptual 
grounds  . 
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FIGURE   8  -   STUDY    AREA     STRATIFICATION    : 
CENTRAL    AND    NON   CENTRAL 
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Model  stability  is  one  of  the  desirable  products  of 
conceptual  validity.   Relationships  that  are  not  conceptually 
valid,  if  established  from  today's  data,  are  more  apt  not 
to  hold  in  the  future.   Predictive  equations  of  trip  gener- 
ation should  hold  for  the  future,  in  order  to  have  any 
forecasting  capability. 

Another  facet  of  conceptual  validity  is  the  sign  of 
the  regression  coefficient.   Because  of  collinearity  in  the 
variables,  the  coefficient  of  one  of  the  independent  vari- 
ables could  be  contrary  to  the  theoretical  relationship  and 
this  condition  might  be  statistically  acceptable.   In  spite 
of  this,~it  was  decided  to  delete  those  variables  whose 
coefficient  had  a  sign  contrary  to  conceptual  expectations. 
This  should  increase  the  statistical  validity  of  the  model 
as  it  tends  to  reduce  the  effects  of  collinearity. 


Siinpl  i  ci  ty .   To  keep  the  models  as  simple  as  possible, 
unnecessary  transformations  of,  and  interactions  among  the 
original  independent  variables  were  avoided.   Interactions 
beyond  the  product  of  two  independent  variables  were  con- 
sidered difficult  to  interpret;  thus  were  avoided  except  if 
the  third  variable  of  the  product  was  the  dummy- vari abl e 
defining  the  location  of  a  zone  in  the  central  or  non 
central  areas . 

Keeping  the  structure  of  the  model  as  simple  as  possible 
by  not  going  to  higher  order  interactions  curtails  the  propa- 
gation of  measurement  errors  [73].   Another  aspect  of 
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simplicity  is  parameter  parsimony.   While  it  is  valuable  to 
include  all  the  relevant  independent  variables,  thus  reduc- 
ing specification  errors,  it  is  doubtful  if  it  would  be 
advantageous  to  do  so  when,  as  is  the  case  for  transporta- 
tion studies,  the  input  data  are  inherently  plagued  by 
measurement  errors.   As  an  emphasis  of  this  research,  the 
number  of  independent  variables  in  the  model  was  kept  to  a 
minimum. 

Stabi 1 i  ty ■   In  order  that  the  developed  models  be 
stable  over  a  time  period,  the  prerequisite  for  allowing 
a  variable  to  enter  the  model  was  a  hypothesized  causal 
relationship  rather  than  a  mere  correlation. 

Stability  was  also  sought  over  the  range  of  the  values 
of  the  independent  variables.   This  could  be  quite  a  diffi- 
cult criterion  to  account  for  during  model  building.   A 
study  of  the  range  of  the  independent  variables  for  the  fore- 
cast year  was  undertaken,  and  possible  problem  zones  were 
identified.   Recommendations  will  be  made  to  ameliorate  this 
condition. 


Sens  i  ti  vi  ty ■   It  is  desirable  that  the  response  or  the 
dependent  variable  be  sensitive  to  changes  in  each  of  the 
independent  variables  in  the  model.   The  cost  of  adding  one 
more  independent  variable  would  not  be  justifiable  if  the 
dependent  variable  is  not  sensitive  to  changes  in  the  added 
independent  variable.   The  sensitivity  of  the  dependent  vari 
able  to  each  of  the  independent  variables  in  the  model  was 
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tested  by  calculating  the  standardized  regression  coefficients 
(the  regression  coefficients  multiplied  by  the  ratio  of  the 
standard  deviation  of  the  independent  variable  under  con- 
sideration to  the  standard  deviation  of  the  dependent  vari- 
able). 

Statistical  Considerations 


Stepwise  Regression.   The  computer  program  used  by  this 
research,  as  mentioned  earlier,  v;as  "Bf1D-2R,  stepwise  re- 
gression" [72].   Several  procedures  are  available  to  develop 
multiple  regression  models.   The  "tear-down"  or  "backward 
elimination"  method  starts  with  a  model  containing  all  the 
available  independent  variables,  and  subsequently  eli  m  i  n  a  t  e  s 
some  of  the  independent  variables  until  a  model  with  pre- 
described  statistical  features  is  reached.   The  "build-up" 
or  "forward  selection"  procedure  strives  for  a  similar  final 
outcome;  however,  working  in  the  opposite  direction  by  in- 
serting one  more  independent  variable  at  a  time.   Stepwise 
regression  is  an  improved  version  of  forward  selection  pro- 
cedure [14].   The  independent  variables  in  the  model  are 
re-examined  at  the  end  of  each  step.   The  variable,  that 
might  have  been  the  best  single  variable  to  enter  at  an 
earlier  step,  might  prove  to  be  unnecessary  at  a  later  stage 
because  of  the  relationship  between  it  and  other  variables 
now  in  the  equation.   Thus,  at  each  step,  the  partial  F-test 
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for  each  variable  in  the  equation  was  evaluated  and  compared 
to  a  preselected  percentage  point  of  the  appropriate  F  dis- 
tribution.  Stepwise  regression  evaluates  the  contribution 
of  each  independent  variable  in  the  model  at  the  end  of 
each  step,  regardless  if  the  independent  variable  has  entered 
at  the  last  step  or  at  any  earlier  step. 


Partial  F  or  Sequential  F-Test.   By  far,  the  most 
important  statistic  in  conjunction  with  regression  analysis 
is  the  multiple  coefficient  of  determination  (R  ).   It 

measures  the  proportion  of  total  variability  in  the  dependent 

2 
variable  explained  by  the  regression  model;  R   varies  between 

zero  and  one;  a  value  of  zero  indicates  a  complete  lack  of 

fit,  while  a  value  of  one  implies  a  perfect  correlation.  In 

stepwise  regression,  a  test  is  needed  at  each  step  to  check  if 

2 
the  increase  in  R   contributed  by  each  added  independent  vari- 
able in  the  equation  is  significantly  different  than  zero. 
The  following  F-statistic  tests  that  the  contribution  of  the 


k  independent  variables  is  significantly  greater  than  zero. 

F, 
where 


^'""•""^    (1  -  R^)/(n  -  k  -  1 


n    =  number  of  observations, 

k   =  number  of  independent  variables,  and 

2 

R,   =  coefficient  of  multiple  determination  of  a  model 

with  k  independent  variables. 
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The  calculated  F-statistic  is  compared  to  a  tabulated 

F,    ,       ,  v/here  a  is  the  probability  of  type  I  error,  or 
k>n-K-i,|  CL 

the  level  of  significance.   The  level  of  significance  chosen 
should  depend  on  the  consequences  of  rejecting  a  true 
hypothesis.   The  level  of  significance  for  including  a  vari- 
able was  set  at  0.010,  and  for  deleting  a  variable  at  0.005. 
The  selection  of  these  values  is  based  on  acceptance  of  a 
relatively  high  risk  of  including  a  variable  which  does  not 
belong.   Once  this  variable  has  been  accepted,  there  is  a 
lower  risk  acceptable  for  its  retention  in  the  equation 
based  on  the  entry  of  other  independent  variables. 

The  blind  use  of  the  F-test  may  result  in  developing 
a  regression  model  which  involves  more  independent  variables 
than  are  of  practical  significance.   In  transportation  stud- 
ies, the  number  of  observations  is  large  resulting  in  an 

F-statistic  which  is  statistically  significant  even  when  the 

2 
absolute  increase  is  R   is  very    small.   The  criterion  of  a 

2 

significant  increase  in  R   proved  to  be  superfluous  in  the 

majority  of  the  cases;  as  other  criteria  such  as  simplicity, 
parsimony,  and  reasonableness  controlled  the  number  of  vari- 
ables to  be  included  in  the  model. 


Standard  Error  of  Estimate .   Another  statistic  of 
interest  is  the  standard  error  of  the  estimate  (s).  It  is 
the  square  root  of  the  residual  mean  square.   The  smaller 
the  value  of  this  statistic  the  more  precise  the  predictions 
would  be.   The  criterion  of  reducing  s  must  be  used 
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cautiously;  since  s  can  be  made  small  by  including  enough 

2 
parameters  in  the  model,  just  as  R   can  be  increased 

As  more  independent  variables   are  included   in   the 

equation,  the  decrease  in  s  will  be  at  a  decreasing  rate. 

Reduction  of  s  is  desirable  if  many  degrees  of  freedom  for 

error  are  remaining  [14]. 

Another  way  of  looking  at  the  reduction  in  s  is  to 

consider  it  in  relation  to  the  dependent  variable  --  namely, 

as  a  percentage  of  the  mean  value  of  the  dependent  variable. 

Standard  error  of  estimate  as  a  percentage  of  the  mean  of 

the  dependent  variable  is  referred  to  as  the  coefficient  of 

variation  (C.V.). 

t-Test  on  Regression  Coefficients .   It  is  sometimes 
desirable  to  test  if  each  of  the  estimated  regression  parame- 
ters are  significantly  different  than  zero.   The  ratio  of 
each  regression  coefficient  to  its  standard  error  is  distrib- 
uted as  student-t.   If  the  regression  coefficient  of  one  of 
the  independent  variables  does  not  pass  the  t-test,  it  can 

be  deleted  from  the  equation. 

2 
The  three  criteria  of  R  ,  s,  and  significance  of  the 

regression  coefficient  are    not  independent.   Usually,  the 

2 
decision  can  be  made  on  the  basis  of  R   alone. 


Model  Identification 
The  first  set  of  models,  that  was  developed  by  this 
investigation,  was  a  rerun  for  each  of  the  13  dependent 
variables  using  the  same  independent  variables  established 
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by  IRTADS  for  their  equations.   Some  of  IRTADS  equations 
were  developed  using  389  zones  out  of  the  395  [38].   There, 
probably,  was  a  good  reason  for  excluding  some  of  the  zones 
from  the  analysis  by  IRTADS.   The  excluded  zones  and  the 
reason  for  their  exclusion  were  not  available  for  this  in- 
vestigation.  Therefore,  and  in  the  interest  of  compatibil- 
ity and  comparability,  data  from  all  the  395  zones  were 
used  to  re-estimate  the  parameters  of  the  models  developed 
by  IRTADS.   Those  models,  essentially  developed  by  IRTADS, 
were  used  as  a  basis  to  compare  with  other  developed  models. 

A  second  set  of  models  that  include  relative  accessi- 
bility variables  were  attempted  for  each  of  the  13  dependent 
variables.   It  was  not,  however,  possible  to  develop  satis- 
factory models  for  each  of  the  13  dependent  variables.   The 
first  two  sets  of  models  were  developed  with  data  from  the 
395  zones  with  no  distinction  relative  to  location  in  the 
central  or  non  central  areas.   Two  more  sets  of  models  were 
developed;  one  corresponding  to  the  set  developed  by  IRTADS, 
the  other  to  the  set  of  models  developed  by  this  investiga- 
tion.  The  models  of  these  two  sets  contained  a  dummy- 
variable  defining  the  location  of  a  zone  in  the  central  or 
non  central  areas. 

Thus,  basically,  four  sets  of  models  were  developed. 
Two  were  without  any  relative  accessibility  variables  among 
their  independent  variables;  one  of  those  two  was  the  set 
developed  by  the  traditional  procedures  for  IRTADS;  the 
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second  set  contained  a  dummy- var i abl e  which  defined  the  zone 
location  and/or  some  of  the  interaction  of  the  dummy- 
variable  with  the  other  independent  variables  in  the  equa- 
tion.  Of  the  remaining  two  sets,  each  had  relative  accessi- 
bility variables  and,  in  addition,  one  was  calibrated  with 
stratified  data. 

Figure  9  suggests  a  system  to  identify  the  developed 
models.   Models  for  each  dependent  variable  in  every    set 
were  not  possible  to  develop.   Cells  of  Figure  9,  where  an 
abbreviation  appears,  indicate  that  an  adequate  model  was 
developed  for  that  dependent  variable  in  the  corresponding 
set. 

Results  of  [lodel  Development 
The  results  of  developing  the  four  sets  of  trip  genera- 
tion models  are  presented  below.   The  four  sets  will  be  re- 
ferred to,  as  suggested  in  Figure  9,  as  follows: 

1.  Without  accessibility,  unstratified:   Set  W-U. 

2.  Without  accessibility,  stratified:   Set  W-S. 

3.  With  accessibility,  unstratified:   Set  A-U. 

4.  With  accessibility,  stratified:   Set  A-S. 
Reference  is  made  to  Appendix  A  for  the  definition  of  the 
symbols  that  appear  in  the  following  discussion. 


IRTADS  Trip  Generation  Models:   Set  W-U 
The  results  of  re-estimating  the  parameters  of  the 
IRTADS  models  are  presented  only  in  summary  form.   The  details 
of  model  development  are  not  reported  because  it  did  not 
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involve  decisions  by  this  investigation.   The  model  calibra- 
tion was  done  for  the  same  set  of  independent  variables 
using  the  complete  data  from  the  395  zones  of  the  study 
area.   Table  2  summarizes  the  set  W-U,  and  Table  3  gives 
the  summary  statistics. 

Activity -Accessibility  Models:   Set  A-U 
In  developing  trip  generation  models  with  accessibility 
variables,  the  independent  variables  which  proved  significant 
in  developing  set  W-U  were  allowed  to  enter  in  addition  to 
the  appropriate  relative  accessibility  variables;  meaning- 
ful cross  products  of  socio-economic  variables  and  accessi- 
bility variables  were  also  considered. 

It  is  restated  here  that  the  number  at  the  end  of  the 
abbreviation  of  an  accessibility  variable  represented  the 
trip  purpose  whose  friction  factor  was  used  in  generating 
the  specific  accessibility  variable.   The  numbers  one  to 
six  represented  the  folio vn'ng  trip  purposes:   home-based 
work,  home-based  shop,  home-based  school,  home-based  other, 
non  home-based  work-oriented,  or  non  home-based  non  work- 
oriented  respectively. 


Model  A-U-1  (HBWKP).   This  model  is  summarized  in 
Table  4.   The  only  accessibility  variable  considered  was 
the  accessibility  to  total  employment  (AEMPTOTl  )  .   The  cross 
products  of  AEMPTOTl  with  labor  force  (LF)  and  with  cars 
(CARS)  were  also  considered.   Only  the  cross  product  of 
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AEMPTOTl  by  CARS  entered  the  model.   Its  contribution,  as 

2 
measured  by  the  increase  in  R   and  the  reduction  in  the 

standard  error  of  the  estimate,  is  very    little.   Ho me -based 
work  productions  are  expected  to  be  inelastic  to  the  condi- 
tion of  the  transportation  system.   The  trip  from  home  to 
work  is  essential  and  regular.   The  frequency  of  this  trip 
is  least  affected  by  the  whims  of  the  trip  maker  and  the 
environmental  conditions. 

Model  A-U-2  (HBSHPP).   It  was  not  possible  to  develop 
this  model.   The  relative  accessibility  variable  considered 
was  that  to  retail  floor  area  (ARTLFLR2).   Neither  ARTLFLR2 
nor  any  of  its  cross  products  with  population  of  a  zone, 
dwelling  units  in  a  zone,  single  family  dwellings  in  a  zone, 
or  the  cars  garaged  in  a  zone  proved  significant.   This 
investigation  was  not  able  to  develop  a  model,  of  home-based 
shop  trips  from  a  zone,  that  had  accessibility  variables. 

Model  A-U-3  (HBSCLP).   Table  5  summarizes  this  model. 
The  only  accessibility  variable  considered  to  be  meaningful 
in  developing  a  model  of  home-based  school  trip  productions 
was  the  accessibility  to  educational  floor  area  (AEDFLR3). 
Cross  products  of  AEDFLR3  with  school  age  population,  d\'/ell- 
ing  units,  single  family  dwellings,  or  cars  were  also  con- 
sidered.  Only  SFD*AEDFLR3  appeared  in  the  chosen  model, 
second  to  CARS,  and  with  a  negative  coefficient.   The  nega- 
tive sign  of  the  coefficient  could  be  understood  if  it  is 
realized  that  a  high  accessibility  to  educational  floor 
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indicates  probably  spatial  proximity  to  schools;  thus  more 
school  trips  are  apt  to  be  walking  trips. 

Model  A-U-4  (HBOTRP).   It  was  not  possible  to  develop 
this  model.   The  accessibility  variables  considered  in  con- 
junction with  home-based  other  trip  productions  were  the 
accessibility  to  total  employment,  accessibility  to  retail 
floor  area,  and  the  accessibility  to  educational  floor  area. 
The  cross  products  of  each  of  the  accessibility  variables 
with  the  variables  population,  cars,  or  single  family 
dwellings  were  also  considered.   None  of  the  accessibility 

variables  or  its  cross  products  could  be  satisfactorily  in- 

2 
eluded.   Their  contribution  to  the  increase  in  R   or  the 

reduction  in  the  standard  error  of  estimate  would  have  been 

little.   The  signs  of  the  regression  coefficient  of  the 

accessibility  variables  would  have  been  negative  despite 

their  positive  first  order  correlation  with  the  dependent 

variable  because  the  relative  accessibility  variables  were 

correlated  to  the  independent  variables  population  and  cars. 

Model  A-U-5  (NHBWKP).   This  model  is  summarized  in 
Table  6.   Accessibilities  to  total  employment,  retail  floor 
area,  and  educational  floor  area  were  considered.   The 
cross  products  of  the  three  accessibility  variables  by 
total  employment,  retail  employment,  and  service  eniploynent 
respectively  were  also  made  available.   Other  cross  products 
could  have  been  considered,  but  because  of  the  inter- 
correlation  of  most  of  the  independent  variables  only  those 
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three  were  considered.   The  first  variable  to  enter  the 
equation  was  the  cross  product  of  total  eniployment  by 
accessibility  to  total  employment.   It  also  ranked  highest 
as  to  the  sensitivity  of  the  dependent  variable  to  the  vari- 
ous independent  variables  in  the  model. 

Model  A-U-6  (NHBNWA).   Table  7  presents  a  summary  of 
the  model.   Only  two  accessibility  variables  were  consider- 
ed:  the  accessibility  to  retail  floor  area  and  the  accessi- 
bility to  educational  floor  area.   Only  those  two  were  con- 
sidered logical  to  investigate  in  conjunction  with  non 
home-based  non  work-oriented  trip  productions.   The  cross 
products  of  each  of  the  accessibility  variables  by  total 
employment,  retail  floor  area,  and  educational  floor  area 
were  also  considered.   In  the  developed  model  only  the 
cross  product  of  total  employment  by  the  accessibility  to 
educational  floor  was  present.   It  entered  at  the  third  step 
and  ranked  fourth  as  to  the  sensitivity  of  the  dependent 
variable  to  the  independent  variables  in  the  model. 

Model  A-U-7  (TOTP).   This  model  is  summarized  in  Table 
8.   All  the  independent  variables  used  in  the  separate 
trip  purpose  production  models  were  allowed  to  enter.   Among 
the  accessibility  variables  only  the  cross  product  of  labor 
force  by  accessibility  to  total  employment  entered  the  final 
model  in  the  fourth  step.   Its  contribution  to  the  statisti- 
cal strength  of  the  model  is  quite  small. 
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Model  A-U-8  (HBWKA).   This  model  is  presented  in 
Table  9.   Eleven  independent  variables  were  allowed  to  enter 
this  model.   The  eleven  variables  were  constituted  of  one 
accessibility  variable:   the  accessibility  to  labor  force, 
five  socio-economic  and  demographic  variables:   total  em- 
ployment, retail  employment,  service  employment,  retail 
floor  area,  and  educational  floor  area,  and  the  five  cross 
products  of  the  accessibility  variable  with  each  of  the 
socio-economic  variables.   The  final  model  included  two 
independent  variables:   total  employment  and  the  cross 
product  of  retail  floor  by  the  accessibility  to  labor  force. 
The  contribution  of  the  latter  is  very    little  compared  to 
the  former. 


Model  A-U-9  (HBSHPA).   The  model  is  summarized  in 
Table  10.   The  independent  variables  considered  were  three 
accessibility  variables,  three  socio-economic  and  land  use 
variables,  and  their  nine  cross  products.   The  accessibili- 
ties to  population,  dwelling  units,  and  single  family  dwell- 
ings were  the  three  accessibility  variables.   Retail  employ- 
ment, service  employnient,  and  retail  floor  area  were  the 
socio-economic  and  land  use  variables.   The  chosen  model 
had  only  two  independent  variables.   The  first  variable  to 
enter  was  the  cross  product  of  accessibility  to  single 
family  dwellings  by  retail  employment  and  the  second  vari- 
able was  retail  floor  area.   The  variable  including  the 
accessibility  measure  entered  first  and  contributed  the  most 
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to  the  statistical  strength  of  the  model. 

Model  A-U-IQ  (HBOTRA).   Table  11  presents  a  summary  of 
this  model.   Four  accessibility  variables  were  considered: 
accessibility  to  population,  accessibility  to  dwelling 
units,  accessibility  to  single  family  dwellings,  and  access- 
ibility to  cars.   The  cross  products  of  retail  floor  area 
and  educational  floor  area  by  access ibilites  to  population, 
single  family  dwellings,  and  cars  were  also  considered.   The 
cross  product  of  retail  floor  area  by  accessibility  to 
single  family  dwellings  entered  at  the  sixth  step  and  was 
the  last  independent  variable  to  be  included. 

Model  A-U-11  (NHBWKA).   This  model  is  summfrized  in 
Table  12.   The  accessibilities  to  total  employment,  retail 
floor  area,  and  educational  floor  area  together  with  their 
cross  products  by  retail  floor  area  and  educational  floor 
area  were  considered.   The  selected  model  of  non  home-based 
work-oriented  trip  attractions  had  seven  independent  vari- 
ables.  The  cross  product  of  total  employment  by  accessi- 
bility to  total  employment  entered  in  the  fourth  step,  and 
the  cross  product  of  retail  floor  area  by  accessibility  to 
total  employment  entered  in  the  sixth  step.   Those  two 
variables  ranked  third  and  sixth  respectively  in  the  sensi- 
tivity of  the  dependent  variable. 


Model  A-U-12  (NHBNWA).   This  model  is  presented  in 
Table  13.   Among  the  independent  variables  allowed  to  enter 
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were  two  accessibilities:   to  retail  floor  area,  and  to 
educational  floor  area;  and  their  cross  products  by  total 
employment,  retail  floor  area,  and  educational  floor  area. 
The  cross  product  of  total  employment  by  accessibility 
to  educational  floor  area  entered  the  model  at  the  third 
step  but  ranked  first  as  to  the  sensitivity  of  the  dependent 
variable.   The  cross  product  of  educational  floor  area  by 
accessibility  to  educational  floor  area,  and  the  variable 
of  accessibility  to  educational  floor  area  entered  in  the 
fifth  and  sixth  steps  respectively,  and  ranked  fifth  and 
sixth  respectively. 


Model  A-U-13  (TQTA).   Table  14  summarizes  this  model. 
All  the  accessibility  variables  used  in  developing  models 
of  trip  attractions  together  with  all  the- other  independent 
variables  were  available.   Many  alternative  models  of  total 
trip  attractions  could  be  built  because  of  the  inter- 
correlations  in  the  independent  variables.   The  selected 
model  was  a  five-variable  model.   The  variable  entering  at 
the  fifth  step  was  a  cross  product  of  retail  employment  by 
the  accessibility  to  single  family  dwellings  calculated 
using  the  friction  factors  for  the  shopping  trips.   The 
same  variable  ranked  third  as  to  the  sensitivity  of  total 
trip  attractions  to  the  five  independent  variables  in  the 
model  . 
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The  Use  of  Dummy-Van"  abl  es 
Stratification  of  the  study  area  into  central  and  non 
central  areas  in  the  analysis  required  the  use  of  a  dummy- 
variable.   The  dummy- vari abl e  represents  a  qualitative 
variable  which  is  unmeasurable  on  a  continuous  numerical 
scale  [74].   This  variable  is  given  a  value  of  one  if  the 
zone  belongs  to  the  central  area  and  zero  if  the  zone  be- 
longs to  the  non  central  area.   Therefore,  the  coefficient 
of  the  dummy- vari abl e  represents  a  comparison  of  zones  in 
the  central  area  to  zones  in  the  non  central  area  [75]. 
Stated  differently,  the  coefficient  of  the  dummy-variable 
measures  the  comparison  in  the  impact  of  the  zone  location 
on  the  response  variable. 

A  Geometrical  Interpretation .   The  geometrical  inter- 
pretation of  the  use  of  dummy- vari  abl  es  v;i  1 1  be  illustrated 
by  an  example.   Consider  the  case  of  simple  linear  regres- 
sion: 


Y  -  a  +  b.,X^  +  b^Z 


where 


Y    =  the  dependent  variable, 

a   =  the  constant  term  of  the  equation, 

X,   =  a  continuously  measurable  independent  variable, 

b,   =  the  regression  coefficient  of  X,, 

Z    =  a  dummy- vari abl e ,  it  equals  one  for  zones  in 

the  central  area  and  zero  for  zones  in  the  non 

central  area  ,  and 
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b„   =  the  regression  coefficient  of  the  dummy- vari abl e 
Geometrically  the  preceding  equation  could  be  plotted  in 
cartesian  coordinates  as  illustrated  in  Figure  10. 

Consider,  now,  the  same  regression  equation  with  a 
cross  product  term.   Thus, 


Y  =  a  +  b^X^  +  b^Z  +  b^X^Z 


where 

Y  ,  a  ,  X  -,  ,  b  1  , 

Z,  and  bp     are  the  same  as  previously  defined:  and 

b^        =    the  regression  coefficient  of  the  cross 
product  term. 
The  above  relationship  could  also  be  plotted  in  cartesian 
coordinates  as  shown  in  Figure  11. 

As  can  be  seen  from  Figures  10  and  11,  a  dummy- vari abl e 
in  the  regression  model  connotes  a  shift  in  the  relationship 
depending  on  the  group  to  which  the  observation  belongs. 
A  du mmy- vari abl e  and  its  cross  products  indicate  both  a  shift 
in  the  intercept  and  a  change  in  the  slope.   Of  course,  a 
regression  model  might  contain  the  cross  product  term  only 
together  with  the  other  independent  variables,  indicating 
a  change  in  the  slope  without  any  shift  in  the  intercept. 

It  should  be  noted  that  interactions  of  the  dummy- 
variable  with  other  independent  variables  present^' a  problem; 
because  leaving  the-m  out  may  reduce  the  model  precision,  but 
if  included  the  model  becomes  unwieldy. 
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b|>0    ,  b2<0  ,b3=  0,Z  =1 
b|>0  ,  b2<0  ,b3<0,Z:^l 


CONTINUOUS    INDEPENDENT    VARIABLE 


FIGURE    II    _    LINEAR     REGRESSION     WITH    A     DUMMY. 
VARIABLE    AND    ITS   CROSS     PRODUCT 


no 


Statistical  Tests  of  Significance.   In  order  to  estab- 
lish a  criterion  for  the  statistical  validity  of  stratifying 
the  input  data,  statistical  tests  are  needed.   Stratified 
data  were  used  for  both  the  trip  generation  models  of  the 
IRTADS  type  and  those  developed  by  this  investigation  that 
included  accessibility  variables.   The  two  sets  with 
stratified  data  were  referred  to  as:   set  W-S  and  set  A-S 
respectively.   The  same  independent  variables  that  were 
chosen  for  sets  W-U  and  A-U  were  allowed  to  enter  plus  the 
previously  defined  dummy- vari abl e  ,  and  all  the  cross  pro- 
ducts of  the  dummy- va ri abl e  by  the  other  independent  vari- 
ables. 

In  order  to  establish  a  rule  for  limiting  the  size  of 
a  regression  model,  another  more  conservative  criterion 
was  utilized  in  addition  to  the  partial-F  test.   The  fol- 
lowing statistical  test  on  the  significance  of  the  net 
contribution  of  considering  the  factor  of  location  was  uti- 


lized.  Consider  the  following  F-ratio 


(R^  -  R2)(n  -  k., 


k^  -  1 


1  -  RA)(k^) 


where 


=  number  of  independent  variables  representing 

the  factor  of  zone  location, 
=  number  of  the  independent  variables  other  than 

those  representing  the  factor  of  location, 
=  number  of  observations. 
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R«    =  coefficient  of  multiple  determination  for  equa- 
tion with  (k,  +  kp)  variables,  and 

2 

Rp,    =  coefficient  of  multiple  determination  for  equa- 
tion with  kp  variables  [75]. 
"Again,  the  significance  of  the  net  relationship  depends 
heavily  on  the  number  of  observations  and  on  the  number  of 
variables."  [75,  p.  381].   The  level  of  significance  was 
conservatively  set  at  a  =    0.0005. 

The  t-test  used  for  testing  if  an  individual  regression 
coefficient  is  significantly  different  from  zero,  is  of  low 
utility  in  the  case  of  dummy- vari abl es  [75].   This  should 
be  clear  because  the  choice  of  the  values  of  one  and  zero 
to  indicate  if  a  zone  belongs  to  the  central  or  the  non 
central  areas  respectively,  is  arbitrary.   Testing  whether 
the  coefficient  of  a  d ummy- vari abl e  is  significantly  dif- 
ferent from  zero  was  in  effect  testing  for  the  significance 
of  the  difference  between  the  two  classes. 


IRTADS  Models  with  Dummy-Variables:   Set  H-S 
For  each  of  the  13  models,  the  dummy- vari abl es  and  each 
of  its  cross  products  by  the  independent  variables  were 
allowed  to  enter.   Thus,  for  a  model  with  k  independent 
variables  in  set  W-U,  the  corresponding  model  in  set  W-S 
could  have  at  most  2k  +  1  independent  variables.   In  most 
of  the  models,  ho v/ ever,  the  number  of  independent  variables 
did  not  reach  2k  +  1.   The  possibility  of  almost  doubling 
the  number  of  independent  variables  is  not  a  detriment  to 
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parsimony,  because  one  component  of  each  of  the  added  cross 
products  was  already  in  the  equation.   The  only  added  cost 
is  due  to  the  collection  of  the  information  necessary  to 
classify  zones  into  central  and  non  central.   Thus  strati- 
fication is  believed  to  enable  a  more  efficient  use  of  the 
other  available  independent  variables. 

Model  W-S-1  (HBWKP).   This  model  is  summarized  in  Table 
15.   Only  the  cross  product  term  was  added.   The  negative 
sign  of  the  cross  product  term  could  indicate  that  the  rate 
of  home-based  work  productions  for  a  zone  in  the  central 
area  per  labor  force  is  lower  than  the  corresponding  rate  for 
a  non  central  zone . 

Model  W-S-2  (HBSHPP).   This  model  was  not  developed  be- 
cause neither  the  dummy- vari abl e  nor  any  of  its  cross  pro- 

2 
ducts  increased  R   significantly.   A  possible  interpretation 

is  that  the  rates  of  home-based  shop  productions  per  car  or 

per  single  family  dwelling  were  not  significantly  different 

in  the  central  and  non  central  areas. 


Model  W-S-3  (HBSCLP).   The  model  which  is  summarized 
in  Table  16  has  only  the  cross  oroduct  term  added.   The 
negative  sign  of  the  regression  coefficient  shows  that  home- 
based  school  productions  per  car  in  the  non  central  area 
are,  almost,  twice  that  of  the  central  area. 

Model  W-S-4  (HBOTRP) .   The  nodel  is  summarized  in  Table 
17.   It  does  not  correspond  exactly  to  model  W-U-4,  which 
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used  the  variable  population  (POP).   Instead,  the  variable 
single  family  dwelling  (SFD)  v/as  used  in  W-S-4.   The  cross 
product  of  SFD  and  the  dummy-vari abl e  was  the  only  addi- 
tional variable. 

Model  W-S-5  (NHBWKP).  The  model  is  presented  in  Table 
18.  The  cross  product  of  dwelling  units  by  the  dummy- 
variable  and  the  dummy-vari abl e  itself  were  not  included  in 
the  model.  The  other  three  cross  products  entered  consecu- 
tively in  the  fifth,  sixth  and  seventh  steps.  No  interpre- 
tation of  the  signs  of  the  regression  coefficients  v/ill  be 
attempted  because  the  relationships  are  probably  confounded 
due  to  the  inter-correlation  of  the  independent  variables. 

Model  W-S-6  (NHBNWP ) .   The  model  is  summarized  in 
Table  19.   It  includes  cross  products  of  the  dummy- vari abl e 
by  retail  floor  area  and  by  retail  employment.   Both  of 
these  had  negative  regression  coefficients  indicating  a 
greater  rate  in  the  increase  of  non  home-based  non  v/ork- 
oriented  productions  due  to  an  increase  in  retail  floor  area 
or  retail  enployment  in  the  central  area  than  in  the  non 
central  area.   The  dummy- vari abl e  itself  entered  at  the 
last  reported  step  with  a  positive  coefficient,  indicating 
a  shift  in  the  intercept  in  addition  to  the  different  slopes 
for  zones  in  the  central  and  non  central  areas.   The  posi- 
tive sign  of  the  regression  coefficient  of  the  dummy-vari abl e 
could  be  interpreted  that  the  location  of  a  zone  in  the 
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central  area  increases  the  potential  of  producing  non  home- 
based  non  work-oriented  trips.   Ho v^ ever,  this  interpretation 
is  very    naive  because  of  the  complex  relationship  among  the 
other  seven  independent  variables.   The  dummy-vari abl e  in 
this  model,  actually,  has  a  damping  effect  on  the  response  of 
the  model  to  changes  in  retail  floor  area  and  retail 
empl oyment . 

Model  W-S-7  (TOTP )  .   This  model  is  summarized  in  Table 
20.   The  last  two  steps  reported  include  the  cross  products 
of  retail  floor  area  and  total  employment  by  the  dummy- 
variable.   Again,  any  attempt  to  interpret  the  sign  of  the 
regression  coefficients  could  be  invalidated  by  the  complex 
inter-correlations  of  the  independent  variables. 


Model 


HBWKA).   The  model  is  summarized  in  Table 


21.   It  shows  that  the  rate  of  home-based  work  attractions 
per  total  employment  is  higher  in  the  non  central  area 
than  in  the  central  area. 


Model  W-S-9  (HBSHPA)  .   Table  22  is  a  summary  of  this 
model.   Two  cross  product  terms  are  in  the  model.   Both  are 
with  negative  regression  coefficients  indicating  the  dif- 
ferent drav/ing  powers  of  activities  in  the  central  area 
compared  to  activities  in  the  non  central  area. 

Model  W-S-10  (HBOTRA)  .   This  model  is  summarized  in 
Table  23.   The  only  cross  products  included  in  the  model  were 
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those  of  the  dummy- vari abl e  with  each  of  retail  floor  area 
and  educational  floor  area.   The  regression  coefficients  in 
both  cases  were  negative.   This  indicates  that  attraction 
rates  of  different  land  uses  are  different,  and  that  the 
rate  for  the  same  land  use  differs  with  location.   In  the 
presented  model  of  home-based  other  attractions,  the  rates 
of  trip  attractions  per  unit  of  educational  floor  area  or 
retail  floor  area  were  higher  for  non  central  locations 
than  for  central  locations  for  constant  values  of  other 
independent  variables  in  the  model. 

Model  W-S-1 1  (NHBWKA ) .   The  model  is  summarized  in 
Table  24.   The  cross  products  of  each  of  the  four  independent 
variables  from  set  W-U  by  the  dummy-vari abl es  appear  in 
this  model. 

Model  W-S-1 2  (NHBNWA )  .   The  model  is  summarized  in 
Table  25.   It  included  in  addition  to  two  cross  product 
terms  the  dummy- vari  abl  e  itself.   Model  W-S-12  is  sii'iilar 
to  model  W-S-6  in  this  respect.   Only  those  two  models  of 
set  W-S  included  the  pure  dummy- va ri abl e  in  addition  to 
cross  product  terms  indicating  a  shift  in  the  intercept  from 
central  to  non  central  locations  in  the  non  home-based  non 
work-oriented  trip  purpose. 


Model  W-S-1 3  (TOTA )  .   This  model  is  summarized  in 
Table  26.   Three  cross  product  terms  appear  in  the  model.   It 
would  be  unsafe  to  interpret  the  signs  of  the  regression 
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coefficients  due  to  the  inter-correlation  of  the  independent 
vari  abl es . 

Ac  ti  vi  ty- Access  i  bi  1  i  ty  Models  vn'th  Dummy-Variables: 

Set  A-S 

The  models  of  set  A-U  that  were  developed  with  accessi- 
bility variables  in  them  were  recalibrated  by  allowing  in 
addition  to  the  independent  variables  already  in  the  equa- 
tion a  dummy- vari abl e  defining  location  and  its  cross  pro- 
ducts by  the  other  independent  variables.   Additional 
accessibility  or  other  independent  variables  were  only  al- 
lowed as  cross  products  involving  the  dummy- vari abl es  . 

Model  A  - S -  1  ( H  B  W  K  P ) .   It  was  not  possible  to  develop 
this  model.   The  model  of  home-based  work  productions  was 
least  possible  to  improve  by  introducing  accessibility  vari- 
ables or  by  stratification. 

Model  A-S-2  (HB5HPP).  This  model  was  not  developed 
because  it  v.' as  not  possible  to  develop  the  corresponding 
model  A-lJ-2. 


Model  A-S-3  (HBSCLI 


The  model  as  summarized  in  lable 


27  indicates  that  the  rate  of  home-based  school  productions 
per  car  in  a  zone  is  lower  for  the  central  area  than  for  the 

non  central  area.   The  addition  of  the  cross  product  of 

2 
cars  by  the  dummy- vari abl e  increased  the  R   of  this  model 

substantially. 
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Model  A-S-4  (HBOTRP).  This  nodel  was  not  developed 
because  it  was  not  possible  to  develop  the  corresponding 
model  A-U-4. 

Model  A-S-5  (NHBWKP) .   This  nodel  is  summarized  in 
Table  28.   The  cross  products  of  the  dummy- vari abl e  by  re- 
tail employment  and  service  employment  appear  in  the  model 
and  both  with  a  negative  regression  coefficient.   A  three- 
term  product  of  total  employment  by  the  accessibility  to 
total  employment  by  the  dummy-variable  appears  in  the  model 
with  a  positive  regression  coefficient.   After  accounting 
for  the  other  variables  in  the  model,  the  rate  of  non  home- 
based  work  productions  per  total  employment  by  the  accessi- 
bility to  total  employment  in  the  central  area  is  approxi- 
mately three  times  that  in  the  non  central  area. 

Model  A-S-6  (NHBNWP).   A  summary  of  the  model  appears 
in  Table  29.   The  cross  products  of  the  dummy-variable  by 
retail  floor  area,  retail  employment,  service  employment, 
and  total  employment  are  present  in  the  model.   It  is  almost 
futile  to  try  to  interpret  the  signs  of  the  regression  coef- 
ficients of  the  dummy- vari abl e  cross  products  because  of 
the  inter -correlation  of  the  independent  variables. 

Model  A-S-7  (TOTP )  .   The  model  is  summarized  in  Table 
30.   As  can  be  seen,  the  contribution  of  accessibility  vari- 
ables and  of  the  dummy- vari abl e  to  the  improvement  of  the 
model  of  total  productions  is  quite  small. 
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Model  A-S-8  (HBWKA).   This  nodel  is  summarized  in 
Table  31.   The  cross  product  of  total  enipl  oyi'ient  by  the 
dummy-vari abl e  was  added.   It  indicates  that  the  rate  of 
ho me -based  work  attractions  per  total  employment  is  higher 
for  zones  of  the  non  central  area  than  for  those  of  the 
central  area.   This  is  valid  only  after  accounting  for  the 
third  term  in  the  model  which  is  a  cross  product  of  retail 
floor  area  by  the  accessibility  to  labor  force. 


Model  A-S-9  (HBSHPA).   As  summarized  in  Table  32,  the 
model  includes  the  triple  product  of  the  dummy-variable  by 
retail  employment  by  the  accessibility  to  single  family 
dwellings.   From  the  sign  and  value  of  the  regression  coef- 
ficients, it  could  be  said  that  after  accounting  for  retail 
floor  area  in  a  zone  the  attractive  power  of  the  zone  to 
home-based  shopping  trips  is  increased  if  the  zone  lies 
in  the  non  central  area  and  is  decreased  if  the  zone  lies 
in  the  central  area  as  a  function  of  the  product  of  the 
accessibility  of  the  zone  to  single  family  dwellings  by  the 
retail  employment.   To  understand  this  relationship  the 
results  of  model  W-S-9  are  recalled.   The  latter  model  showed 
that  the  attraction  of  shopping  trips  per  unit  of  retail 
floor  area  or  per  retail  employment  was  higher  in  the  non 
central  area  than  in  the  central  area.   Since  the  first 
order  correlation  of  retail  employment  to  retail  floor  area 
was  relatively  high  (r  =  0.810),  this  same  relationship  was, 
probably,  reflected  also  in  model  A-S-9  in  the  sign  and 
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value  of  the  regression  coefficient  of  the  variable  that 
included  the  product  of  retail  employment  by  the  dummy- 
variable  namely  the  triple  product  of  relative  accessibility 
to  single  family  dwellings  by  retail  employment  by  the 
dummy- vari  able. 

Model  A-S-10  (HBOTRA).   As  summarized  in  Table  33, 
this  model  included  three  cross  products  of  the  dummy- 
variable.   The  sign  and  magnitude  of  the  regression  coef- 
ficients suggest  that  service  employment  and  educational 
floor  area  are  weak  contributors  to  the  attraction  of  home- 
based  other  trips  in  the  zones  of  the  central  area. 

Model  A-S-11  (MHBWKA).   As  shown  in  Table  34,  the  model 
contained  cross  products  of  the  dummy- vari abl e  with  three 
out  of  the  ten  independent  variables  in  the  model.   In 
spite  of  the  complex  relationships  among  the  indeoendent 
variables  in  the  model,  the  si  an  and  magnitude  of  the  re- 
gression coefficients  suggest  that  retail  employment  and 
retail  floor  area  have  very  little  effect  on  the  attractive 
power  of  a  zone  in  the  central  area  for  non  home-based  work 
trips. 

Model  A-S-1 2  (NHBNMA)  .   This  model  is  presented  in 
Table  35.   The  same  observation  could  be  I'lade  about  the 
effect  of  retail  employment  and  retail  floor  area  as  noted 
above  on  model  A-S-11.   It  is  further  observed  that  the 
interaction  of  total  employment  by  the  accessibility  to 
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educational  floor  area  will  increase  the  attractive  power 
of  zones  in  the  central  area  for  non  home-based  non  work- 
oriented  trips. 


Model  A-$-13  (TOTA).   The  node!  is  summarized  in  Table 
36.   It  would  be  difficult  to  draw  conclusions  on  the  ef- 
fect of  different  independent  variables  in  this  model.   How- 
ever, it  can  be  observed  that  the  stratification  of  the  data 
into  central  and  non  central  contributes,  probably,  more 
to  the  soundness  of  the  relationship  than  the  inclusion  of 
the  accessibility  variables. 
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CHAPTER  V.   DISCUSSION  OF  RESULTS 


Comparison  of  the  Developed  Models 
Four  sets  of  trip  generation  models  were  developed  by 
this  investigation.   One  of  which  was  a  recalibration  of 
the  set  of  models  developed  by  IRTADs.   The  i  mproveinents 
introduced  by  each  of  the  other  three  sets  over  the  set 
developed  by  IRTADS,  along  the  traditional  approach, were 
investigated.   Because  the  desirability  of  a  model  is 
determined  by  other  than  statistical  criteria  alone,  the 
scope  of  the  comparison  was  not  restricted  to  statistical 
considerations.   These  criteria  were  discussed  in  Chapter 
IV. 

Since  siiiiilar  criteria  governed  the  building  of  the 
four  sets  of  models,  the  statistical  comparisons  were 
limited  to  the  two  criteria  of  coefficients  of  multiple 
detenii nat i on  (R  )  and  the  coefficient  of  variation  (C.V.). 
It  was  not  always  possible  to  conduct  statistical  tests  of 
the  significance  of  the  increase  in  R"^  of  the  m odels  of  each 
set  compared  to  the  models  of  all  others.   Comparisons  of 
two  models  were  made  where  the  independent  variables  of 

one  were  a  subset  of  the  independent  variables  of  the  other. 

2 
The  significance  of  the  increase  in  R   achieved  in  con- 
sidering the  factor  of  location  of  a  zone  could  be  tested 
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using  an  F-statistic  as  described  in  Chapter  IV,  however, 
this  was  not  always  possible.   When  the  independent  vari- 
abiles of  one  model  are  not  a  subset  of  the  independent 
variables  of  the  other  model  or  when  the  two  models  do  not 

have  the  same  number  of  independent  variables,  the  compari- 

2 
son,  based  upon  changes  in  R  ,  must  be  made  entirely  on 

personal  judgment. 

The  scope  of  the  comparisons  included: 

1.  improvements  achieved  by  introducing  relative 
accessibility  variables  to  the  basic  IRTADS 
models,  that  is,  models  of  set  A-U  versus  models 
of  set  W-Ui 

2.  improvements  achieved  by  calibrating  the  models 
with  data  stratified  according  to  the  zone  loca- 
tion over  the  basic  IRTADS  models,  that  is,  models 
of  set  W-S  versus  models  of  set  W-U; 

3.  improveiiients  achieved  by  calibrating  the  models 
containing  relative  accessibility  variables  with 
data  stratified  by  location  over  similar  models 
calibrated  with  unstratified  data,  that  is,  models 
of  set  A-S  versus  models  of  set  A-U;  and 

4.  improvements  achieved  by  introducing  both  relative 
accessibility  variables  and  calibrating  the  model 
with  stratified  data  over  the  basic  IRTADs  models, 
that  is,  models  of  set  A-S  verus  models  of  set 
W-U. 
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Coraparison  of  Models  with  Accessibility  Variables 
Versus  Basic  IRTADS  Models  (Set  A-U  Versus  Set  W-U) 

Introducing  relative  accessibility  variables  to  models 
of  trip  generation  was  proposed  by  this  investigation  to 
increase  the  conceptual  validity  of  the  models.   It  was 
hypothesized  that  trips  generated  by  a  zone  are  a  function 
of  the  relative  accessibility  of  the  zone.   In  addition  to 
the  conceptual  validity  considerations,  the  statistical 
strength  of  the  relationship  was  tested.   Whenever  the  in- 
dependent variables  in  a  model  of  set  W-U  were  a  subset  of 
the  independent  variables  of  the  corresponding  model  of  set 

A-U,  the  statistical  test  of  the    significance  of  the  increase 

2 
in  R  was  applicable.   Table  37  summarizes  the  statistical 

2 

test  of  the  significance  of  the  increase  in  R  .   The  sta- 
tistical test  was  not  applicable  for  models  where  the  trip 
purpose  is  the  only  entry  in  Table  37. 

HBWKP:  Models  with  and  without  Accessibility.   Model 
A-U-1  included  two  more  independent  variables  than  model 

W-U-1,  only  one  of  the  added  independent  variables  was  a 

2 

relative  accessibility  measure.   Thus  the  increase  in  R 

atti^ibuted  to  it  was  considered  in  the  test.   The  increase 

2 
in  R   was  not  significant  at  a  =  0.0005. 


HBSCLP:   Models  with  and  without  Accessibility.   The 

2 
increase  in  R   was  significant  at  a  =  0.0005. 
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NHBWKP:   Models  with  and  without  Accessibility.   The 

addition  of  a  relative  accessibility  variable  did  not  in- 

2 
crease  the  R  ,  even  with  more  independent  variables  in  the 

equation.   However,  it  is  important  to  note  that  the  rela- 
tive accessibility  independent  variable  was  the  first  to 
enter,  since  it  had  the  highest  si  in  pie  correlation  with  the 
dependent  variable  among  all  the  independent  variables 
which  were  considered  for  this  model. 

NHBNWP:   Models  with  and  without  Accessibility .   The 

model  with  relative  accessibility  variables  did  not  achieve 

2 

a  higher  R  . 

TOTP:   Models  with  and  without  Accessibility.   No  in- 

2 

crease  in  R  was  achieved  by  including  relative  accessibil- 
ity variables.   The  model  W-U-7  is  probably  a  more  desirable 
model . 


HBWKA:   Models  with* and  without  Accessibility.   The 

2 
increase  in  R   attributed  to  including  a  relative  accessi- 
bility variable  was  not  significant  at  a  =  0.0005. 

HBSHPA:   Models  with  and  without  Accessibility.   I n- 

2 
eluding  a  relative  accessibility  variable  increased  the  R 

for  the  same  number  of  Independent  variables,  and  the 

standard  error  of  the  estimate  was  reduced.   A  statistical 

test  was  not  possible  because  the  independent  variables  of 

one  of  the  models  were  not  a  subset  of  the  Independent 

variables  of  the  other. 
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HBOTRA:   Models  with  and  without  Accessibility.   The 

2 
increase  in  R   due  to  introducing  a  relative  accessibility 

variable  was  significant  at  a  =  0.0005. 

NHBWKA:   Models  with  and  without  Accessibility.   Intro- 
ducing relative  accessibility  variables  did  not  Increase 

2 
R  .   The  variable,  dwelling  units,  in  model  W-U-11  was  not 

allowed  in  model  A-U-11  because  a  meaningful  relationship 

with  non  home-based  trips  was  not  justifiable. 

NHBNWA:   Models  with  and  without  Accessibility.   The 

2 

increase  in  R   due  to  introducing  relative  accessibility 

variables  was  small.   A  statistical  test  on  the  significance 
of  the  increase  was  not  possible,  the  Independent  variables 
of  model  W-U-12  were  not  a  subset  of  the  Independent  vari- 
ables of  model  A-U-12. 


TOTA:   Models  with  and  withou^  Accessibility.   The 
two  models  di'ffer  only  in  the  last  Independent  variable  to 
enter.   At  the  fifth  step,  model  U-U-13  accepted  the  inde- 
pendent variable:   retail  employment  (EMPRTL).   At  the  same 

step,  model  A-U-13  accepted  the  variable:   EMPRTL  *  ASFD2 

2 
which  increased  the  R   by  0.0383  as  compared  to  0.0089 

achieved  by  adding  EMPRTL.   The  model  A-U-13  is  probably 
more  desirable.   The  product  of  retail  employment  and  the 
relative  accessibility  to  single  family  dwellings  (calcu- 
lated using  friction  factors  for  the  shopping  trips)  con- 

2 
t  r  1  b  u  t  e  d  in  ore  to  the  Increase  in  R   than  the  variable  retail 

employment  alone. 
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Comparison  of  Stratified  Models  without  Accessibility 

Variables  Versus  Basic  IRTADS  Models 

(Set  W-S  Versus  Set  W-U) 

Table  38  summarizes  the  statistical  test  of  significance 

described  in  CHAPTER  IV.   Stratification  of  the  zones  into 

2 

central  and  non  central  areas  increase  the  R   significantly 

at  an  a  level  of  0.0005  for  all  the  models,  except  for  the 
trip  purpose  of  HBWKA. 


Comparison  of  Stratified  Models  with  Accessibility 

Variables  Versus  Unstratified  Models 
with  Accessibility  Variables  (Set  A-S  Versus  Set  A-U) 

The  purpose  of  the  comparison  was  to  assess  the  increase 

2 
in  R   achieved  by  stratifying  the  zones  by  location  for  trip 

generation  models  which  included  relative  accessibility  vari- 
ables.  The  tests  of  significance  are    suinmarized  in  Table  39. 

Model  A-S-1  was  not  possible  to  develop  as  were  models 

2 
A-U-2  and  A-U-4.   The  increase  in  R   achieved  by  the  nodels 

of  set  A-S  over  the  models  of  set  A-U  were  all  significant 

at  a  =  0.0005  except  for  HBWKA. 

Comparison  of  Stratified  Models  with  Accessibility 

Variables  Versus  Basic  IRTADS  Models 

(Set  A-S  Versus  Set  W-U) 

The  purpose  of  this  comparison  was  to  investigate  the 

overall  improvement  attributed  both  to  including  relative 

accessibility  variables  and  stratifying  zones  by  location 

in  set  A-S  over  the  basic  IRTADS  models  of  set  W-U.   Again, 

statistical  tests  were  only  possible  if  the  independent 

variables  of  the  model  of  set  W-U  were  a  subset  of  those  of 
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the  corresponding  model  of  set  A-S.   Table  40  sumn arizes  the 
statistical  test  of  significance. 

Models  A-S-1,  A-S-2,  and  A-S-4  were  not  possible  to 
develop.   Models  where  a  statistical  test  was  applicable 

are  shown  with  complete  entries  in  Table  40.   The  increase 

2 
in  R  was  significant  at  the  o.    level  of  0.0005  for  the  fol- 
lowing trip  purposes:   HBSCLP,  HBSHPA,  HBOTRA,  and  TOTA. 

2 

No  significant  increase  in  R   was  achieved  for  HBWKA.   Those 

models  for  which  a  statistical  test  did  not  apply  d.re    dis- 
cussed below. 


Models  of  NHBWKP .   Model  A-S-5  included  nine  independent 
variables,  five  of  which  were  socio-economic  variables. 
Model  W-U-5  included  four  socio-economic  variables  which  were 

not  a  subset  of  the  independent  variables  of  model  A-S-5. 

2 
The  increase  in  R   due  to  including  four  independent  vari- 
ables, measuring  relative  accessibility  and  specifying  zone 
location,  in  addition  to  the  five  socio-economic  variables 
in  model  A-S-5  was  significant  at  a  =  0.0005. 

Models  of  NHBNWP.   Model  A-S-6  had  a  higher  R^  and 

lower  S.E.  than  those  of  model  A-U-6.   However,  a  statistical 

2 
test  of  the  significance  of  the  increase  of  R   was  not 

possible  to  perform. 

Models  of  TOTP.  The  two  models  are  identical  in  the 
first  three  steps.  At  the  fourth  step,  model  W-U-7  added 
the  variable  population  while  model  A-S-7  added  3  more 
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variables  at  the  fourth,  fifth,  and  sixth  steps.   The  com- 

2 
bined  effect  of  the  increase  in  R   due  to  those  3  added 

variables  was  approximately  equivalent  to  including  the 
variable  population.   Thus  model  A-S-7  could  not  be  consid- 
ered as  an  improvement  over  model  W-U-7. 

Models  of  NHBWKA.   In  developing  model  A-S-11,  the  in- 
dependent variable  dwelling  units,  which  entered  second  in 
model  W-U-11,  was  excluded  because  the  trip  purpose  being 
modeled  was  non  home-based  work,  and  neither  ends  of  the 
trip  would  be  expected  to  start  or  terminate  at  a  dwelling 

unit.   The  contribution  of  the  variable  dwelling  units  to 

2 

the  R   of  model  W-U-11  was  large.   Even  though  model  A-S-11 

2 
with  10  independent  variables  had  a  higher  R   than  model 

W-U-11  with  4  independent  variables,  no  statistical  test  of 

significance  could  be  used. 

Models  of  NHBNWA.   Model  A-S-12  with  7  independent  vari^ 

2 
ables  achieved  a  higher  R   than  model  W-U-12  with  5.   The 

increase  could  not  be  statistically  tested. 

Overall  Comparisons 

2 
Values  of  R  ,  C.V.,  and  the  number  of  independent  vari- 
ables for  each  model  of  the  four  developed  sets  are  summar- 
ized in  Table  41.   A  discussion  of  the  improvements,  if  any, 
introduced  to  the  models  of  the  thirteen  trip  purposes 
f ol 1 ows  . 
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Model s  of  HBWKP .   Introducing  relative  accessibility 
variables,  and/or  stratifying  the  zones  by  location  did  not 
make  any  practical  improvement  over  the  model  of  set  W-U 
with  a  single  independent  variable. 

Models  of  HBSHPP.   No  acceptable  alternatives  to  model 
W-U-2  were  possible  to  develop. 

2 
Models  of  HBSCLP .   The  increase  in  R   was  significant 

at  a  =  0.0005  for  models  W-S-3,  A-U-3,  and  A-S-3  over  W-U-3. 

Including  relative  accessibility  variables,  stratification 

by  zone  location,  and  a  combination  of  those  two  actions 

improved  the  modeling  of  this  trip  purpose. 

Models  of  HBCTRP .   Mo  acceptable  alternatives  to  model 
W-U-4  were  possible  to  develop. 

Models  of  N  H  B  W  K  P  .   They>e  was  a  significant  increase  in 

9 

R   attributed  to  stratification  of  zones  by  location  at  ct  = 
0.0005  for  both  models  W-S-5  and  A-S-5  over  W-U-5  and  A-U-5 
respectively.   Adding  only  relative  accessibility  variables 
did  not  improve  the  model.   The  improvement  due  to  strati- 
fication was  more  evident. 


Models  of  NHBNWP.   Stratification  increased  R   signifi- 
cantly in  models  W-S-6  and  A-S-6  over  models  W-U-6  and 
A-U-6  respectively.   Introducing  relative  accessibility  vari 
ables  did  not  contribute  as  much  as  stratification  did. 
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Models  of  TOTP.   Both  stratification  and  relative 
accessibility  variables  did  not  drastically  improve  other 
models  over  W - U -  7 . 


Models  of  HBWKA.   Both  stratification  and  relative 
accessibility  variables  did  not  drastically  improve  other 
models  over  W-U-8. 

Models  of  HBSHPA ,   This  trip  purpose  had  the  weakest 
model  in  the  IRTADS  set.   Stratification  strengthened  the 
model  but  relative  accessibility  variables  did  not  contribute 
as  much . 

2 
Models  of  HBOTRA .   It  was  possible  to  increase  R   sig- 
nificantly by  either  stratifying  with  respect  to  location, 
or  including  relative  accessibility  variables,  or  both. 

Models  of  NHBWKA.   Stratification  of  both  n.odels  W-U-11 

2 
and  A-U-11  increased  R   significantly.   Including  relative 

accessibility  variables  only,  without  stratification,  did 

not  improve  the  model. 

Models  of  NHBNWA.   The  increase  in  R^  of  model  W-S-12 

over  model  W-U-12  was  significant.   Adding  relative  accessi- 

2 

bility  variables  without  stratification  did  not  increase  R 

as  much  as  stratification  alone  achieved.   However,  a  seven 
independent  variable  model  with  both  relative  accessibility 
variables  and  stratification  was  slightly  better  than  the 
eight  independent  variable  model  with  stratification  only. 
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Models  of  TOTA .   Stratification  increased  R   signifi- 
cantly for  both  models  with  and  without  relative  accessi- 
bility variables.   Model  A-S-13  with  both  relative  accessi- 

2 
bility  variables  and  stratification  increased  R   significantly 

over  model  W-U-13. 
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CHAPTER  VI.   AN  APPLICATION  OF  THE  DEVELOPED  MODELS 


The  four  sets  of  developed  models  were  solved  with 
1985  forecast  of  the  independent  variables.   The  productions 
and  attractions  forecasted  for  each  zone  by  purpose,  as  pre- 
dicted by  each  of  the  four  sets  of  models,  were  compared  to 
detect  trends  and  establish  if  they  were  significantly  dif- 
ferent from  each  other.   The  forecasts  of  the  socio-economic 
variables,  and  the  minimum  time  skim  trees  of  the  proposed 
future  network  were  obtained  from  IRTADS.   The  procedure 
described  in  CHAPTER  IV  v/as  used  to  generate  the  1985  rela- 
tive accessibility  variables  that  entered  the  regression 
equations. 

A  regression  i^iodel  applies  in  the  region  of  values 
covered  by  the  range  of  each  of  the  independent  variables 
in  the  model.   In  the  case  of  multiple  regression  the  in- 
ference space  could  be,  indeed.  Quite  restricted;  especially, 
when  considering  the  joint  space  defined  simultaneously  by 
the  range  of  all  the  independent  variables  in  the  model. 
Defining  this  inference  space  could  be  an  unwieldy  combina- 
torial problem. 

In  the  case  of  trip  generation  analysis  with  data 
aggregated  at  the  zone  level,  the  applicability  of  the  re- 
gression models  for  forecasting  might  be  reduced  if  some  of 
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the  zones  expect   growth  that  results  in  forecast  of  inde- 
pendent variables  outside  the  range  of  the  base  year    values, 
In  addition  to  range,  the  distribution  of  the  values  of  an 
independent  variable  in  the  base  year  and  the  forecast  year 
could  be  of  importance  to  the  applicability  of  the  regres- 
stion  model.   The  strength  of  a  model  lies  in  the  region 
where  most  of  the  values  of  the  independent  variables  fall. 
Theoretically,  independent  variables  with  uniform  distribu- 
tions provide  regression  models  of  equal  strength  over  the 
whole  range  of  the  values  of  the  independent  variables. 


1985  Forecast  of  the  Socio-Economic  Variables 
The  following  observations  about  the  forecast  values 
of  the  socio-economic  variables  were  by  no  means  extensive 
or  complete.   The  values  of  each  of  the  ten  socio-economic 
variables  for  each  zone  were  tabulated  in  frequency  distri- 
butions for  the  base  year  and  the  forecast  year.   It  v/as 
observed  that  the  shape  of  the  distributions  was  similar  for 
the  two  years;  however,  the  distribution  of  the  1985  values 
was  shifted  to  the  right,  that  is,  less  frequent  in  the  low 
range  of  values  and  more  frequent  in  the  high  range  of  val- 
ues.  This  is  to  be  expected,  it  is  because  of  the  urban 
growth  of  the  area.   Appendix  D  oresents  the  tabulated  fre- 
quency distributions  and  the  corresponding  histograms. 

The  values  of  the  1985  forecasted  socio-economic  vari- 
ables in  each  zone  were  checked  against  the  ranges  of  1964 
values.   The  summary  of  this  check  is  presented  in  Table  42. 
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It  shows  the  zones  where  values  of  the  1985  forecast  socio- 
economic variables  are  outside  the  range  of  the  1964  values. 
Twenty  three  zones  of  the  395  zones  of  the  study  area  had 
between  one  and  five  of  the  socio-economic  variables  that 
had  forecasted  values  outside  the  range  of  1964  values. 
The  total  employment  forecasted  value  for  1985  exceeded  the 
1964  range  for  only  zone  170.   The  educational  floor  area 
forecasted  values  for  1985  exceeded  the  1964  range  for  each 
of  zones  72,  82,  326,  327,  328,  and  361. 

It  should  be  noted  that  only  three  zones  out  of  23 
belong  to  the  central  area.   Those  were  zones  326,  327, 
and  328,  which  are  adjacent  and  constitute  district  number 
71.   The  forecasted  educational  floor  area  in  each  of 
them  was  outside  the  range  of  1964  values. 

A  third  group  of  zones  constituted  of  zones  100,  102, 
108,  110,  175,  177,  209,  232,  240,  314,  324,  364,  382,  390, 
and  391  had  forecasted  va-lues  outside  the  range  of  the 
values  of  demographic  variables  that  characterize  residential 
areas . 

The  latter  group  of  zones  were  observed  to  be  at  the 
borders  of  the  study  area,  and  were  of  very    large  land  area 
compared  to  the  other  zones  of  the  study  area.   The  percent- 
age of  land  in  urban  use  for  all  the  zones  of  this  group 
was  below  50?^  which  indicates  underdevelopment  in  1964,  but 
high  potential  for  growth  by  1985.   The  high  growth  potential 
and  the  large  area  of  those  zones  resulted  in  a  forecast  of 
socio-economic  variables  outside  the  range  of  1964  data. 
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Other  1985  Independent  Variables 
The  developed  models  included  seven  relative  accessi- 
bility variables.   They  were: 

1.  ALFl; 

2.  AEMPTOTl; 

3.  ASFD2; 

4.  AEDFLR3; 

5.  ASFD4; 

6.  AEMPT0T5i  and 

7.  AEDFLR6. 

Those  seven  variables  were  generated  for  1985  using  the 
procedure  described  in  CHAPTER  IV.   The  values  of  the  fric- 
tion factors  were  assumed  to  stay  the  same  for  the  forecast 
year.   Preliminary  research  has  indicated  that  such  an 
assumption  was  warranted  [76]. 

The  definition  of  the  central  and  non  central  areas 
was  assumed  the  same  for  the  forecast  year  as  for  the  base 
year.   A  check  was  not  possible  because  forecast  values  for 
all  the  variables  that  were  used  to  reach  a  stratification 
were  not  available.   The  procedure  described  in  CHAPTER  IV 
depended  essentially  on  the  ranking  of  the  study  area  dis- 
tricts with  respect  to  three  measures  of  land  use  and  land 
use  intensity.   The  underlying  assumption  in  using  the  same 
stratification  for  1985  was  that  the  ranking  of  the  districts 
with  respect  to  the  three  decision  variables  would  be, 
essentially,  the  same  in  1985.   The  exact  ranking,  of  course. 
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was  not  the  criterion;  rather,  the  belonging  to  the  upper 
or  lower  quart ile.   Thus  the  implied  assumption  was  not 
over  demanding. 


1985  Productions  and  Attractions  Forecast 
Computer  programs  were  written  to  solve  the  developed 
■regression  equations  of  trip  generation  for  the  forecast 
values  of  the  independent  variables.   Four  sets  of  fore- 
casts for  each  of  13  trip  purposes  for  every    zone  in  the 
study  area  were  obtained.   In  the  event  that  the  forecast 
value  happened  to  be  negative,  due  to  a  negative  constant 
term  in  the  equation,  the  forecast  value  was  set  to  zero. 
The  programs  supplied  outputs  in  the  form  of  tabulated 
printouts  and  punched  cards. 

The  results  of  the  forecasts,  by  trip  purpose  for  each 
of  the  zones  of  the  study  area,  were  compared  against  each 
other.   The  objective  was  to  test  if  the  four  sets  of  models 
resulted  in  significantly  different  forecasts.   Five  fore- 
cast comparisons  were  considered: 

1.  models  without  accessibility  and  no  stratification 
(basic  IRTADS  models)  versus  models  with  accessi- 
lity  and  no  stratification  (set  W-U  versus  set  A-U); 

2.  models  without  accessibility  and  no  stratification 
(basic  IRTADS  models)  versus  models  without  accessi- 
bility but  with  stratification  (set  W-U  versus 

set  W-S); 
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3.  models  with  accessibility  and  no  stratification 
versus  models  v/ith  accessibility  and  stratifica- 
tion (set  A-U  versus  set  A-S); 

4.  models  without  accessibility  but  with  stratifica- 
tion versus  models  with  accessibility  and  strati- 
fication (set  W-S  versus  set  (A-S);  and 

5.  models  without  accessibility  and  no  stratification 
(basic  IRTADS  models)  versus  models  with  accessi- 
bility and  stratification  (set  W-U  versus  set  A-S) 

Each  of  the  above  comparisons  were  made  for  all  the  395 
zones  of  the  study  area  as  one  group,  for  the  105  zones  of 
the  central  area  as  a  second  group,  and  for  the  290  zones 
of  the  non  central  area  as  a  third  group.   In  total,  twelve 
comparisons  were  conducted. 

The  Pai  red  t-Tes t 

The  statistical  test  used  for  the  above  comparisons 
was  the  paired  t-test.   In  any  one  comparison  two  sets  of 
forecasts,  of  equal  size  for  each  trip  purpose,  were  under 
consideration.   The  two  sets  to  be  compared  were  denoted  by 
Y  and  Y'.   Each  component  of  the  two  sets  can  be  logically 
paired  with  respect  to  zone  number. 

The  hypotheses  to  be  tested  were: 


H 1  :   y  Y  7*  u  Y 
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The  set  of  the  differences  of  the  paired  observations  v/as 


defined  as  D  =  Y  -  Y',  or  otherwise  stated  as  d. 


y  ■  -  y  • 


where  i  =  1  to  n .   The  above  hypotheses  could  be  stated: 


H^:   Ud  f    0 


Let  d  and  s-r  be  the  estimate  of  w^      and  the  standard  devia- 
tion of  the  estimate  of  pp  respectively,  then 

t  =  3_ 

s_ 

d 
would  be  the  test  statistic.   The  above  hypotheses  were  re- 
j e c t e d  if 


t|  ->_    t(i-a/2)(n-l 


It  could  be  shown  that  a  paired  t-test  is  equivalent  to  a 
randomized  complete  block  design  with  two  blocks.   Equiva- 
lent conclusions  could  thus  be  reached  using  analysis  of 
variance  techniques. 

A  two-tailed  test  was  used  because  there  were  no  £ 
priori  reasons  to  suspect  that  the  forecasts  by  the  models 
of  one  of  the  sets  would  yield  a  higher  or  lower  forecast 
than  that  by  the  models  of  another  specific  set. 

Result  of  the  Tests 
The  mean  value  of  the  differences  between  the  pair  of 
zonal  forecast,  for  each  trip  purpose,  by  models  of  the  two 
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compared  sets  and  the  associated  t-statistic  are  tabulated 
in  Appendix  E.   The  test  results  are  discussed  below. 

All  Zones  of  the  Study  Area 
A  summary  of  the  results  is  presented  in  Figure  12.   The 
five  comparisons  are  discussed  below. 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Models  with  Accessibility  Variables  (Set  W-U  Versus  Set  A-U). 
Among  the  trip  purposes  compared  three  models  did  not  fore- 
cast significantly  different  zonal  productions  and  attrac- 
tions.  These  three  trip  purposes  were:   HBSHPA,  NHBNWA, 
and  TOTA.   It  was  observed  that  the  forecast  for  all  trip 
productions  were  significantly  different  at  a  =  0.05,  and 
on  the  average,  models  with  accessibility  variables  con- 
sistently forecasted  zonal  productions  larger  than  those  by 
basic  IRTADS  models. 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  without  Accessibility  Variables  (Set  W-U 
Versus  Set  W-S).   For  only  two  of  the  trip  purposes,  namely 
HBWKA  and  HBSHPA,  the  differences  in  the  zonal  forecasts, 
on  the  average,  were  significant  at  a  =  0.05.   It  is  observed 
again  that  the  forecasts  by  the  stratified  models  were,  on 
the  average,  larger  than  the  forecasts  by  the  basic  IRTADS 
models,  whenever  the  difference  v/as  significant. 
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TRIP 
PURPOSE 

COMPARISON  OF  FORECASTS  BY  VARIOUS  MODEL  SETS 
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Hq  :   fjiy=  fj.y'  vjQS  rejected  at    a  =  0.05 
Hq  :  fi^  =  fiy'  was  not  rejected  at  a  =0.05 
d    for  indicated  purpose  was  positive 
d      II  II  II  M     negative 

indicates  that  either  or  both   models  do  not  exist 


FIGURE     12-    PAIRED    t_TEST:  SUMMARY  OF   RESULTS 
ALL     ZONES 
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Forecasts  by  Unstratified  Models  with  Accessibility 
Variables  Versus  Forecasts  by  Stratified  Models  with 
Accessibility  Variables  (Set  A-U  Versus  Set  A-S).   Forecasts 
of  only  two  of  the  trip  purposes  were  significantly  dif- 
ferent at  a  =  0.05.   The  two  trip  purposes  were  HBOTRA  and 
NHBNWA. 

Forecast  by  Stratified  Models  without  Accessibility 
Variables  Versus  Forecasts  by  Stratified  Models  with  Accessi- 
bility Variables  (Set  W-S  Versus  Set  A-S)  .   All  trip  purooses 
except  TOTA  had  significantly  (a  =  0.05)  different  fore- 
casts.  No  generalization  was  possible  regarding  the  sign 
of  the  mean  of  the  difference. 


Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  with  Accessibility  Variables  (Set  W-U 
Versus  Set  A-S).   The  mean  differences  between  the  forecasts 
were  significant  at  a  =  0.05  for  all  the  trip  purposes 
tested.   No  generalizations,  however,  were  possible  to  nake 
regarding  the  sign  of  the  mean  of  the  difference. 

Zones  of  the  Central  Area 
The  same  tests  were  conducted  for  the  zones  of  the 
central  area  only.   A  summary  of  the  test  results  is  presented 
in  Figure  13.   The  five  comparisons  are  discussed  below. 
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TRIP 
PURPOSE 

COMPARISON  OF  FORECASTS  BY  VARIOUS    MODEL  SETS 
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was  rejected  at  a  =  0.05 


Hq  •  /iy  =  /^y'  was  not   rejected   at  a  =  0.05 
d     for  Indicated  purpose  was  positive 
d       11  "  "  11     negative 

indicates   that  either  or  both  nnodels  do  not  exist 


FIGURE      13-    PAIRED    t  _  TEST  :    SUMMARY    OF    RESULTS 
ZONES    OF    THE    CENTRAL    AREA 
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Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Models  with  Accessibility  Variables  (Set  W-U  Versus  Set  A-U). 
Forecasts  of  total  productions  were  significantly  different 
at  a  =  0.05,  while  differences  between  forecasts  of  all 
other  trip  purposes  were  not  significant.   Models  with  rela- 
tive accessibility  variables  resulted  in  forecasts  close, 
on  the  average,  to  those  by  IRTADS  models  in  the  zones  of 
the  centra  1  area . 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  without  Accessibility  Variables  (Set  W-U 
Versus  Set  W-S)  .   All  the  compared  trip  purooses  had  signif- 
icantly different  forecasts,  except  for  the  NHBWKP  trip 
purpose.   It  was  observed  that  for  the  trip  purposes  where 
the  differences  were  significant,  the  means  of  the  differ- 
ences for  zones  of  the  central  area  were  positive.   There- 
fore, the  basic  IRTADS  models  would  yield,  on  the  average, 
forecasts  in  the  central  area  higher  than  the  forecasts  by 
the  correspondinq  stratified  models. 


Forecasts  by  Unstratified  Models  with  Accessibility 
Variables  Versus  Forecasting  by  Stratified  Models  with 
Accessibility  Variables  (Set  A-U  Versus  Set  A-.S).   All  the 
differences  were  significant  except  for  the  HBOTRA  trip 
purpose.   Again,  the  mean  differences  for  all  the  trip  pur- 
poses that  had  a  significant  difference  were  positive.   Thus, 
stratification,  probably,  has  the  same  effect  on  models  with 
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relative  accessibility  variables  as  on  models  without  them. 
Stratification  tends  to  result  in  models  that  predict  on 
the  average  lower  trip  productions  and  attractions  for  the 
zones  of  the  central  area. 

Forecasts  by  Stratified  Models  without  Accessibility 
Variables  Versus  Forecasts  by  Stratified  Models  with  Ac- 
cessibility Variables  (Set  W-S  Versus  Set  A-S).   Forecasts 
by  models  of  the  two  compared  sets  were  significantly  dif- 
ferent only  for  the  two  trip  purposes:   TOTP  and  HBOTRA. 
The  sign  of  the  mean  difference  for  the  latter  two  purooses 
was  negative.   Because  only  two  trip  purposes  are  involved, 
it  is  difficult  to  draw  meaningful  inferences. 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  with  Accessibility  Variables  (Set  W-U 
Versus  Set  A-S )  .   The  zonal  forecasts  were  significantly 
different  on  the  average.   It  also  appeared  that  forecasts 
by  the  traditional  models  of  IRTADS  would  result  in  fore- 
casts higher  than  those  predicted  by  the  models  with  rela- 
tive accessibility  and  stratified  data  for  all  but  one  of 
the  trip  purposes,  for  zones  of  the  central  area.   One  would 
expect  that  future  growth  of  the  urban  area  would  cause  a 
decrease  in  the  relative  importance  of  the  central  area, 
and  as  such  it  would  attract  a  lower  proportion  of  the 
tri  ps . 
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Zones  of  the  Non  Central  Area 
The  comparison  of  the  forecasts  for  the  zones  of  the 
non  central  area  are  discussed  below.   The  summary  of  the 
results  of  the  difference  tests  are  presented  in  Figure  14. 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Models  with  Accessibility  Variables  (Set  H-U  Versus  Set  A-U) 
The  differences  between  the  forecasts  by  the  models  of  those 
two  sets  were  significant  at  a    =  0.05  for  all  but  one  of 
the  trip  purposes,  namely  HBSHPA.   No  trends  were  detected 
in  the  sign  of  the  mean  of  the  differences. 

Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  without  Accessibility  Variables  (Set  W-U 
Versus  W-S).   Forecasts  of  HBWKP,  HBSCLP,  NHBWKP,  HBWKA, 
HBSHPA  and  HBOTRA  by  the  models  of  the  two  sets  were  siqnifi' 
cantly  different  at  a  =  0.05.   It  v/as  also  observed  that 
the  mean  differences  for  all  the  above  trip  purposes  were 
negative.   This  indicated  a  probable  trend  that  models  cali- 
brated with  stratified  data  predicted,  on  the  average,  zonal 
productions  and  attractions  higher  than  the  predictions  by 
the  basic  IRTADS  models,  for  the  above  trip  purposes,  in 
the  zones  of  the  non  central  area. 
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d    for  indicated  purpose  was  positive 
d     II  M  II  II      negative 

indicates    that    either    or    both    models    do  not    exist 


FIGURE      14-    PAIRED    t_TEST   :    SUMMARY    OF    RESULTS 
ZONES   OF  THE    NON    CENTRAL    AREA 
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Forecasts  by  Unstratified  Models  with  Accessibility 
Variables  Versus  Forecasts  by  Stratified  flodels  with  Ac- 
cessibility Variables  (Set  A-U  Versus  Set  A-S).   The  models 
of  the  two  sets  resulted  in  forecasts  that  were  signifi- 
cantly different  at  a  =  0.05  for  five  trio  purposes,  namely, 
HBSCLP,  HBWKA,  HBSHPA,  HBOTRA,  and  TOTA .   The  mean  of  the 
difference  was  negative  for  the  five  above  purposes.   This 
substantiated  the  observation  that  models  calibrated  with 
stratified  data  tend  to  predict  higher  trips,  on  the  aver- 
age, for  zones  of  the  non  central  area. 

Forecasts  by  Stratified  Models  without  Accessibility 
Variables  Versus  Forecasts  by  Stratified  Models  with  Ac- 
cessibility Variables  (Set  W-S  Versus  Set  A-S).   Forecasts 
by  the  models  of  the  two  sets  were  significantly  different 
at  a  =  0.05  for  all  but  the  HBOTRA  trip  purpose. 


Forecasts  by  Basic  IRTADS  Models  Versus  Forecasts  by 
Stratified  Models  with  Accessibility  Variables  (Set  W-U 
Versus  Set  A-S)  .   Forecasts  for  all  the  trip  purposes  as 
predicted  by  the  models  of  the  two  sets  were  significantly 
different  at  a  =  0.05. 

It  can  be  noted  from  Figure  15  that,  in  the  non  central 
area,  the  basic  IRTADS  models  (set  W-U)  overf orecas ted  com- 
pared to  stratified  models  with  accessibility  (set  A-S),  on 
the  average,  trips  for  the  following  five  purposes:   NHBWKP, 
HBWKA,  HBSHPA,  NHBWKA,  and  NHBNWA.   For  these  trip  purposes 
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only,  those  zones  which  did  not  individually  have  overfore- 
casted  trips  are  shov/n  in  Figure  15.   These  zones  tend  to 
fall  in  the  vicinity  and  along  the  corridors  defined  by 
major  thoroughfares.   In  other  words  the  proposed  accessi- 
bility stratified  models  (set  A-S)  tend  to  forecast  greater 
numbers  of  trips  in  zones  in  the  non  central  area  for  most 
trip  purposes;  and  for  the  five  purposes  noted  above,  the 
models  still  forecast  greater  numbers  where  the  zones  have 
greater  accessibility  as  provided  by  the  corridors. 


The  Proposed  Trip  Generation  Process 
The  trip  generation  models  proposed  by  this  research 
are  functions  of  the  status  of  the  transportation  system. 
In  an  operational  transportation  study  future  forecasts  of 
trip  generation  would  then  be  affected  by  the  nature  of  the 
proposed  transportation  network.   And  since  the  proposed 
network  should  be  designed  to  serve  future  trip  generation; 
therefore,  an  iterative  process  should  be  followed.   It 
would  be  terminated  when  an  equilibrium  between  the  future 
supply  of  transportation  (proposed  plan)  and  the  demand  for 
transportation  (travel  forecast)  is  reached.   This  iterative 
process  is  illustrated  in  Figure  15. 
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FIGURE     15 


NON    CENTRAL    AREA    ZONES     WHERE     BASIC 
IRTADS     MODELS    DID    NOT    OVERFORECAST 
TRIPS    FOR   :    NHBWKP,    HBWKA,   HBSHPA, 
AND /OR     NHBWKA    TRIP    PURPOSES 
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CHAPTER  VII .   CONCLUSIONS 


Based  on  the  preceding  results  and  analyses  the  follow- 
ing conclusions  can  be  drawn. 

1.  Among  all  of  the  relative  accessibility  variables 
considered,  the  following  variables  were  included 
in  the  trip  generation  models  which  were  developed: 

a.  accessibility  to  employment  in  conjunction  with 
hone-based  work  person-trip  productions  and 

non  home-based  work-oriented  person-trip 
attractions  ; 

b.  accessibility  to  labor  force  in  conjunction 
with  home-based  work  person-trip  attractions; 

c.  accessibility  to  single  family  dwellings  in 
conjunction  with  home-based  shop  person-trip 
attractions  and  home-based  other  person-trip 
attractions;  and 

d.  accessibility  to  educational  floor  area  in  con- 
junction with  home-based  school  productions, 
and  non  home-based  non  work-oriented  oerson- 
trip  attractions. 

The  preceding  accessibilities  were  each  calculated 
with  the  friction  factor  corresponding  to  the  same 
trip  purpose  as  the  i-^iodel  under  consideration. 
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2.  Relative  accessibility  variables  in  trip  neneration 
models  improved  the  statistical  strength  of  models 
of  person-trip  attractions  more  than  that  of  models 
of  person-trip  productions.   Competition  is  a  more 
important  locational  consideration  for  high- 
attraction  zones  which  indicates  their  need  for 
greater  accessibility. 

3.  Calibrating  trip  generation  models  with  data  strati- 
fied according  to  the  location  of  the  zone  in  the 
central  or  non  central  areas  always  improved  the 
statistical  strength  of  the  models  whether  the 
models  had  accessibility  variables  or  not.   Models 
by  hone-based  person-trip  attractions  were  least 
improved  by  stratification,  indicating  substantially 
similar  attracting  characteristics  for  work  trips 

by  zones  in  the  central  and  non  central  areas. 

4.  In  general,  the  statistical  strength  of  the  models 
was  better  achieved  by  stratification  alone  than 

by  including  relative  accessibility  variables  only. 

5.  Models  of  home-based  work  person-trip  productions 
or  attractions  were  improved  least  by  including 
relative  accessibility  variables  and/or  stratifica- 
tion.  This  is  expected  because  work  trips  are 
inelastic  to  trip  length,  due  to  their  regularity 
and  essentiality. 

The  four  sets  of  developed  models  were  solved  with  the 
1985  forecasted  values  of  the  independent  variables.   The 
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forecasts  were  analyzed  to  identify  comparative  forecasting 
trends  of  the  different  models.   The  following  conclusions 
were  drawn  . 

6.  It  was  observed  that  stratified  models  consistently 
forecasted  more  trip  productions  and  attractions 
for  zones  of  the  non  central  area  and  less  for 
zones  of  the  central  area  than  models  without 
stratification.   Stratified  models  are  thus  sensi- 
tive to  the  situation  of  equilibrium  and  saturation 
being  reached  in  the  central  area,  and  also,  the 
faster  rate  of  traffic  growth  in  the  non  central 
area . 

7.  1985  forecasts  of  person-trip  productions  and 
attractions  by  models  with  relative  accessibility 
variables  and  that  were  calibrated  with  stratified 
data  were  significantly  different  than  forecasts 

by  basic  IRTADS  models.   There  was  not  a  detectable 
trend  as  to  the  sign  of  the  mean  difference  between 
zones  of  the  central  and  non  central  areas.   Further 
analysis  indicated  that  stratified  models  with 
relative  accessibility  variables  forecasted  more 
productions  and  attractions  than  forecasts  by  basic 
IRTADS  models,  in  general,  for  zones  located  in  the 
vicinity  and  along  corridors  defined  by  the  major 
thoroughfares  of  the  study  area.   This  reflects  a 
possible  locational  aspect  of  trip  generation  in 
addition  to  the  central-non  central  stratification. 
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From  the  study  of  the  range  of  the  values  of  the 
forecasted  socio-economic  variables  as  comoared 
to  their  range  in  the  base  year,  it  was  concluded 
that  the  anticipated  growth  in  the  outer  parts  of 
the  study  area  should  be  taken  into  consideration 
at  the  outset  of  the  study  namely,  in  defining 
the  zones.   Care  should  be  taken  to  limit  the  size 
of  these  outlying  zones  to  keep  from  having  fore- 
casted values  of  the  socio-economic  variables  that 
are  outside  the  range  of  the  base  year  values. 
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CHAPTER  VIII.   RECOMMENDATIONS  FOR  EXTENSIONS  AND 
FURTHER  RESEARCH 


Based  on  the  preceding  results,  analyses,  and  conclu- 
sions, the  following  extensions  and  further  research  are 
recommended . 

1.  A  study  of  the  effects  of  changes  in  the  transpor- 
tation network  level  of  service  on  the  trip  gener- 
ation characteristics  of  an  urban  area  should  be 
conducted  as  soon  as  data  for  tv/o  points  in  time 
are  available  for  niedium  and  large  size  urban  areas. 

2.  It  should  be  tested  if  the  same  relative  accessi- 
bility variables  identified  by  this  i n ves ti qati on 
would  be  relevant  to  trip  generation  in  other 
urban  areas.   The  research  should  include  urban 
areas  with  a  range  of  sizes,  larger  and  smaller  than 
the  IRTADS  area. 

3.  The  sensitivity  of  the  relationship  of  trip  genera- 
tion to  relative  accessibility  variables  should  be 
investigated  for  different  size  urban  areas,  which 
would  test  the  hypothesis  that  relative  accessibil- 
ity is  a  more  important  factor  in  larger  urban 
areas . 
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7. 


For  larger  urban  areas,  the  importance  of  non 
highway  transit  accessibility  as  a  factor  in  the 
trip  generation  models  should  be  investigated. 
Further  research  is  needed  to  establish  the  cri- 
teria for  differentiating  the  districts  (or  zones) 
between  the  central  or  non  central  strata. 
Consideration  should  be  given  to  the  possibility  of 
a  three  level  stratification:   the  first  level 
being  the  central  or  CBD  area;  the  second  level 
taking  in  those  districts  (or  zones)  which  are  non 
central  and  along  major  transportation  corridors; 
and  lastly,  the  remainder  of  the  non  central  area. 
This  research  should  provide  knowledge  on  the  ef- 
fects on  trip  generation  of  different  urban  growth 
patterns  and  transportation  systems. 
When  formulating  relative  accessibility  variables 
to  be  included  in  nodels  of  auto-trip  attractions, 
the  consideration  of  the  availability  of  parking 
at  the  destination  end  of  the  trip  should  be  in- 
vestigated. 
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APPENDIX  B:   GENERATING  ACCESSIBILITY  VARIABLES 


The  Program  MATCH 


$  IBFTC    MA  TCH 

INTEGER       TREES(513) 

DIMENSION    FRIC(3,'fO),FF(3,513) 

REWIND    1 

REWIND    i^ 

REWIND       8 

READ{5,300) (FRIC(2,J»,J=1,^0) 


CALL   SLITE(O) 

CALL   SLITEll) 

CALL   KDBIN(TREES) 

DO  3  IZON=l,^27 

5 

CALL   RDBIN 

DO  10    J=2,^28 

K=l 

I 

IF(TREES( J) .EQ.K)  GO  TO  20 

IF( TREES! J) .GE.^0)  GO  TO  30 

K=K+1 

GO  TO  1 

20 

FF  (2,J)  =  FRIC(2,K) 

GO  TO  10 

30 

FF(2»J)=FRICl2,'fO> 

10 

CONTINUE 

200 
300 
400 


V.RI  TE(4,200) 

XRI TE(6,400) 

CONTINUE 

CALL    SLITE(2) 

CALL    KDBIN 

ENDFILE    4 

REWIND    4 

FORMAT (I6,19F6 

F0RMAT{12X,F6.2) 

F0RMAT(I6,6F20.2) 

STOP 

END 


IZON, !FF(2, J ), J=2,428) 
IZQN»(FF{2,JI,J=2f402tB0) 


2/(20F6.2) ) 
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$  IBMAP 

BI  NTR 

RCBIN 

SAVE 

SLT 

2 

TRA 

*+2 

TRA 

END 

SLT 

1 

TRA 

RDON 

CLA 

3,^ 

STA 

PLACE 

TSX 

.QPEN.'t 

PZE 

TREE 

TRA 

BACK 

RCON 

TSX 

.READ,«^ 

PZE 

TREE,, ERR 

PZE 

END,, ERR 

PLACE 

lORT 

♦*,,♦♦ 

BACK 

RETURN 

ROBIN 

END 

TSX 

.CLOSE, 4 

PZE 

TREE 

TRA 

BACK 

ERR 

TRA 

SYSDMP 

TREE 

FI  LE 

END 

,CK2,BLK  =  5I3,LOW,BIN,HCLD,  IMPL)  T,  >^OJNT 

$CATA 


Note:   RDBIN  is  a  modified  version  of  a  similar 

SUBROUTINE  that  appeared  in  Reference  [77] 
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"able  Bl  . 

Final 

IRTADS  Fri 

ction  Factors. 
PURPOSE 

1 

2 

3 

'l 

5 

6 

T«AVEL 

T  IME 

HB 

HB 

HB 

HB 

MHB 

MHB 

(  MINUTES) 

WORK 

SHOP 

SCHOOL 

OThFR 

WORK 

MON-WORK 

L 

12.00 

35.00 

85.00 

27.00 

11.50 

25.00 

2 

12.00 

32.00 

66.00 

24.00 

9.80 

15.00 

3 

8.60 

25.00 

46.50 

19.30 

7.80 

9.90 

4 

6.40 

19.00 

32.50 

14.50 

5.90 

6.70 

5 

4.80 

14.00 

23.00 

10.70 

4.50 

4.60 

6 

3.65 

9.70 

16.00 

7.50 

3.20 

3.15 

7 

2.82 

6.05 

1  1.00 

5.10 

2.15 

2.30 

8 

2.25 

3.50 

7.10 

3.35 

1.45 

1.60 

9 

1.85 

2.  10 

4.40 

2.85 

1.05 

1.17 

10 

1.55 

1.35 

2.60 

1.65 

0.76 

0.82 

11 

1.3^ 

0.90 

1.75 

1.30 

0.60 

0.64 

12 

1.15 

0.63 

1.25 

1.00 

0.47 

0.48 

13 

1.02 

0.45 

0.85 

0.79 

0.38 

0.37 

14 

0.89 

0.33 

n.61 

0.62 

0.30 

0.27 

15 

0.80 

0.24 

0.43 

0.49 

0.24 

0.20 

16 

0.72 

0.18 

0.33 

0.41 

0.21 

0.  16 

17 

0.65 

0.15 

0.26 

0.35 

0.17 

0.13 

18 

0.58 

0.  12 

0.20 

0.29 

0.15 

0.10 

19 

0.54 

0.09 

0.15 

0.25 

0.13 

0.08 

20 

0.48 

0.07 

0.12 

0.21 

0.11 

0.07 

21 

0.44 

0.06 

0.10 

0.18 

0.09 

0.06 

22 

0.41 

0.05 

0.08 

0.16 

0.08 

0.05 

23 

0.37 

0.0^ 

0.06 

0.14 

0.07 

0.04 

24 

0.35 

0.03 

0.05 

0.12 

0.06 

0.03 

25 

0.37 

0.03 

0.04 

0.11 

0.06 

0.03 

26 

0.29 

0.02 

0.04 

0.10 

0.05 

0.02 

27 

0.27 

0.02 

0.03 

0.09 

0.05 

0.02 

28 

0.25 

0.02 

0.03 

0.08 

0.04 

0.02 

29 

0.23 

0.02 

0.02 

0.07 

0.0^ 

0.02 

30 

0.21 

0.01 

0.02 

0.06 

0.03 

0.01 

31 

0.18 

0.01 

0.02 

0.06 

0.03 

0.01 

32 

0.17 

0.01 

0.02 

0.05 

0.03 

0.01 

33 

0.16 

0.01 

0.01 

0.05 

0.03 

O.Ol 

34 

0.15 

0.01 

0.01 

0.04 

0.02 

0.01 

35 

0.14 

O.Ol 

0.01 

0.04 

0.02 

O.Ol 

36 

0.  13 

O.Ol 

0.0  1 

0.04 

0.02 

0.01 

37 

0.12 

0.01 

0.01 

0.03 

0.02 

0.01 

38 

0.1 1 

0.01 

0.01 

0.03 

0.02 

0.01 

39 

0.10 

0.01 

0.01 

0.03 

0.02 

0.01 

40 

0.10 

O.Ol 

0.01 

0.03 

O.Ol 

0.01 

Source : 

Reference  [66] . 
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The  Program  ACCESS 


10 


20 


30 


PROGRAM 

ITAPE5=IN 

DIMENSIO 

READ(5,2i 

RE&Dl5,6i 

REWIND    1 

READd  ,li 

WRI TE(6, 

00     10    K= 

XCI ZON.K 

DO    10    J= 

XI  ZON.K 

IF(  IZON.I 

GO     TO    1 

CCNTINUE 

REWIND    I 

REWIND    2 

READ(2,li 

WRI TE(6, 

DO    20    K= 

r'  =  K-^ 

Xd  ZON,K 

DO    20    J= 

X(I ZON.K 

IF  (  IZOiN. 

GO    TO    3 

CCNTI.MUE 

REWIND    2 

REWIND    3 

READ(3,1 

WRI TE(6, 

DO    30    K= 

Xd  ZDN,K 

N=K-7 

DO    30    J  = 

xn  ZQN.K 

IF(  IZON. 

GO    TO    5 

CCNTINUE 

ACCESS! INPUT,UUTPUT,PUNCH,TAPE1,TAPE2, TAPE  3, 

PUT,TAPE6=0UTPUT,TAPE7=PUNCH) 

N  FF«^27),V(39  5,9),X(39  5tl6J,SX(16»,ACC(395,l6) 

00)      (  {\/(  I,  J),J=1,5  ),  1  =  1,  395) 

00)      (  (V(  I  ,J  ),J  =  6,9),  1=1, 395) 

00) IZON,  (FF( J)  , J=l,^27) 

5  00)IZ0N,(FF(J),J=l,40l,80) 

I  ,9 

)=0. 

1,395 

)=X(IZON,K)+V( J,K)«FF( J ) 

EQ.395)  GO  TO  2 


00)  IZON,  (FF( J)  ,  J=  1,427) 
500) IZON, (FF(J),J=1, 40 1,80) 
10,13 

)=0. 

1  ,395 

)=X(IZON,K)+V( J,M)*FF( J ) 

EQ.395)  GO  TO  4 


00)  IZON,  (FF( J)  ,J=1, A27) 

500)iZ0N,(FF(J),J=l,40l,80) 

14,16 

)=0. 

1  ,395 

)=X(IZON,K)+VJ  J,N)«FF( J ) 

E0.395)GO  TO   6 
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60 


70 


80 
100 
200 
300 
^00 


REWIND 

DO    60 

SX( K)= 

DO    60 

SX(K)= 

DC    70 

URI TE( 

DO    70 

ACC  (IZ 

DO    80 

PUNCH 

WRI TEC 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

2  ,*ZONE 

3,*  2 

510F8.2 
62/8X,* 
7*90  * 
8F8.2/7 
9.2/7X, 
07X, *19 
17X,*22 
37X,*25 
37X,*28 
47X,*31 
57X,«34 
67X,*37 
500  FORMAT 
FORMAT 
FORMAT 
STOP 
END 


3 

K=l, 

0. 

IZON 

SX(K 

K=l, 

6,70 

IZON 

ON,K 

K=l, 

300, 

6,A0 

(16, 

(«tX, 

(16/ 

(IHl 

S  TO 

9/ 

/8X, 

60 

,7X, 

X,*l 

*160 


600 
700 


0 

0 

0 

0 

0 

0 

0 

(16, 

CfX, 

(16) 


16 

=  1,3 

)*^X( 

16 

0) 

=  1.3 

)=X( 

16 

K,( 
0) 

19F6 
2F6. 
(  12F 
,/// 

ACT 
3 
/8X, 
*30 
*,10 
*100 
30 

*, 
,10F 
,10F 
,  lOF 
,10F 
,10F 
,10F 
,10F 
18F6 
1F6. 


95 

IZONJ,K) 

K 

95 

IZOM,K)/SX(K )*100. 


ACC(  I 
K,  (AC 
,11  KZ 
0,6X, 
6.2)  ) 
/,  19X 
IVITY 

*00*, 
*,10 
F8.2/ 
*,1 
*,10F 
10F8. 
8.2/7 
8.2/7 
3.2/7 
8.2/7 
8.2/7 
8.2/7 
8.2/7 
.2) 
0,42X 


ZON,K) , IZ0N=1, 395) 
C( IZUN,K ), IZ0N=l,395 
0F6.2)  ) 
2F6.0,  18X,  1F6.0,22X) 


,*IR 
*,  14 

lOX, 
F8.2 
3X,* 
0F8. 
8.2/ 
2/7X 
X,*2 
X,*2 
X,*2 
X,*2 
X,*3 
X,*3 
X,«3 


TADS 
//6X 

9F8. 

/8X, 

70 

2/7X 

7X,* 

f*17 

00 

30 

60 

90 

20 

50 

80 


-PUR 
*ZON 

5 
2/8X 

♦  40 
*,  10 
.♦11 
140 
0   * 
*,  10 

♦  ,  10 

♦  ,  10 
*,  10 

♦  ,  10 

♦  ,  10 
*,  10 


CUE  ACCE 
E 

6 
.*10   *, 

*,  10F8. 
F8.2/8X, 
0   *,10F 

*, 10F8. 
,10F8.2/ 
F8.2/7X, 
F8.2/7X, 
F8.2/7X, 
F8.2/7X, 
F8.2/7X, 
F8.2/7X, 
F8.2/7X, 


SSIBILI 

0 

7 
10F8.2/ 
2/8X,*5 

♦  80  (  , 
8.2/7X, 
2/7X,*l 
7X,*180 

♦  210   ♦ 

♦  240 

♦  270 

♦  300 

♦  330 

♦  360 

♦  390 


TY  DF 
1 

8X,^2 
0  ♦, 
10F8. 
♦  120 
50  ♦ 
♦  10 
,  10F8 
,  10F8 
,  10F8 
,  10F8 
,  10F8 
,  10F8 
,6F8. 


ALL^ 

8    ♦ 
0      ♦, 
10F8. 
2/8X, 
♦  ,10 
,10F8 
F8.2/ 
.21 
.2/ 
.2/ 
.2/ 
.21 
.If 
2) 


,3F6.0, lOX ) 


APPENDIX  C:   DELIMITING  THE  CENTRAL  AREA 
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APPENDIX  D:   FREQUENCY  DISTRIBUTION  OF 
INDEPENDENT  VARIABLES 


Tables  and  histograms  comparing  the  frequency  distri- 
bution of  the  socio-economic  independent  variables  for  the 
survey  year,  1964,  and  the  forecast  year,  1985,  are  presented 
in  the  following  pages. 
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APPENDIX  E:   COMPARISON  OF  TRIP  FORECASTS 

Summary  tables  of  the  comparisons  of  the  trip  fore- 
casts by  the  four  sets  of  developed  models  for  each  of  the 
central,  non  central,  and  the  v/hole  study  area  are  presented 
on  the  following  pages. 

In  the  following  tables  d  represents  the  riean  of  the 
difference  in  the  forecast  by  the  two  models  being  compared, 
and  the  t-statistic  was  calculated  as  follows: 


t  =  ^ 

s 


The  degrees  of  freedom  were  394,  289,  or  104  when  testing 
zones  of  the  whole  study  area,  the  non  central  area,  or 
the  central  area  respectively. 


224 


k.     -> 


r-l 

(\j 

".0 

*o 

in 

<>• 

<T 

o 

fSI 

ir\ 

•f 

00 

h3 


^        ^ 


o 

in 

CO 

o 

ro 

r«- 

o 

in 

si- 

o 

^^ 

(M 

sO 

sO 

CT 

o 

cr 

a 

• 

• 

• 

• 

t 

• 

<NI 

rf\ 

O 

o 

in 

r^ 

>J- 

O 

f—t 

^-^ 

^^ 

fO 

--<       I 


^ 

o 

00 

r»- 

(N 

00 

CsJ 

.^ 

o 

r>- 

^ 

(M 

* 

00 

o 

o 

(NJ 

C 

(M 

00 

^^ 

o- 

in 

r- 

ro 

^-4 

^-4 

>4- 

vO 

• 

¥ 

• 

■Jf 

• 

• 

• 

• 

« 

t 

• 

o 

1 

1 

1 

o 

1 

in 
1 

1 

^H 

O 

^-4 

in 

f\/ 

O 

sO 

•o 

r^ 

«— 1 

>o 

ro 

00 

in 

(T 

>J- 

o 

vO 

f~- 

>J- 

GO 

r^ 

fO 

(M 

cr 

O 

+ 

m 

^ 

o 

-J- 

o 

00 

ro 

vO 

h- 

00 

CO 

• 

•x 

• 

< 

• 

« 

• 

• 

• 

• 

• 

• 

^-4 

CD 

^ 

-r 

00 

vO 

m 

r>- 

f>- 

■O 

pn 

cr 

1 

1 

1 

1 

(T 

1 

1 

1 

1 

1 

LU 

o 


<J- 

(T 

s}- 

CO 

cr 

* 

>r 

* 

t 

00 

1 

*. 

t 
1 

* 

-*• 

ro 

O 

O 

PO 

CT> 

t 

• 

1 

o 
1 

in       "^       "^ 


* 

:¥ 


LU 

a 

Q.  Q. 

>—  Q. 


_c 

oo 

a 

CD 

OD 

CD 

X 

I 

I 

1— '         -J 


225 


'■a 


UJ 

a 

Q.  Q- 


r-<  — 4  00 


r~ 

"r- 

<o 

<o 

00 

CsJ 

r- 

00 

r-< 

•—4 

t 

^ 

o 

• 

;k 

• 

• 

0* 

• 

>o 

^ 

•-^ 

^ 

fSJ 

o 

o 

,-^ 

pri 

1 

1 

1 

u^ 

*,^ 

(M 

ro 

a> 

o 

.—1 

00 

•t 

< 
* 

•J- 

^ 

>* 

■o 

• 

^ 

• 
O 

1 

• 

^^ 

O 

>i- 

<t 

ir\ 

*-4 

(T 

f\j 

o 

(NJ 

cr 

CO 

• 

t 

• 

• 

rri 

1 

o 

1 

o 

rO 

«fl 

f-^ 

(T 

o 

in 

r» 

cr 

cr 

in 

CO 

o 

r^ 

CNJ 

* 

(NJ 

m 

•o 

<r 

<o 

o 

o 

• 

t 

• 

• 

• 

• 

• 

X. 

cr 

ir> 

m 

o 

o 

r>- 

■c 

1 

•—4 

1 

1 

1 

1 

1 

vO 

—1 

r- 

r- 

fO 

•t 

r»- 

•c 

00 

in 

in 

rvi 

o 

■j- 

00 

o^ 

X 

00 

>t 

O 

p~ 

^i- 

• 

^ 

t 

• 

• 

• 

• 

• 

^^ 

^ 

m 

%. 

>r 

#-H 

r\j 

r^ 

in 

ro 

in 

1 

^-4 

-o 

in 

CT 

>— ) 

^^ 

f 

in 

Q. 

a. 

Q. 

a. 

Q. 

Q. 

Q. 

^ 

a: 

i^ 

:x 

^ 

I 

u 

h- 

^ 

z 

i 

(/) 

00 

o 

CD 

CO 

CD 

CO 

CD 

CD 

X 

X 

X 

X 

X 

I 

^ 

z 

< 

<t 

<r 

< 

< 

Q. 

cc 

iC 

:e 

5C 

X 

t— 

i 

CD 

i 

00 

a 

(C 

m 

CD 

CD 

CD 

I 

X 

X 

X 

X 

z 

z 

-o 

in 

in 

(M 

sT 

o 

00 

O 

S- 

•—4 

N- 

in 

in 

ro 

o- 

f<^ 

o 

in 

t 

t 

rvj 

cr 

IM 

in 

.— t 

o 

in 

>o 

17< 

ffl 

o 

".0 

^ 

m 

^H 

cr 

-o 

sO 

~4 

m 

rvi 

.£1 

o 

r^ 

o 

CO 

fV) 

o 

>*• 

y 

in 

*. 

00 

rg 

■t 

-J- 

in 

ro 

00 

f\j 

l«- 

• 

• 

• 

• 

• 

t 

f 

• 

• 

t 

( 

^ 

^^ 

^ 

1 

^^ 

O 

INI 
1 

(M 
1 

«— 4 

o 

1 

—t 

00 

>J- 

po 

>t 

fC 

rvj 

o 

O 

-* 

—1 

fo 

>*• 

in 

t-i 

(NJ 

o 

r-t 

•— ( 

r~- 

00 

1— 1 

CD 

s*- 

^ 

r>- 

^ 

.^ 

rr\ 

r- 

•& 

<t- 

OD 

rf- 

in 

J- 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

1 

"it 

■^ 

00 

1 

00 

•T) 

1 

o 

(\J 

rg 

1 

p»^ 
(NJ 

226 


•■o 


UJ 

o 


'■D 


I/) 

a 
a  a 
—       zc 

a:        Z> 

I—  u. 


r-  - 

o       r 


'<^        — I       — I 


-4  I 

I 


in 

s*- 

O 

(T 

m 

m 

<r 

m 

0^ 

XI 

o 

;K 

r\j 

O 

vO 

r-4 

in 

00 

c> 

■o 

in 

(^ 

(M 

in 

o 

>0 

<o 

(M 

X 

B 

« 

^ 

• 

• 

t 

• 

• 

f 

• 

• 

00 

f\J 

r- 

m 

in 

r^ 

^-4 

c^ 

'^ 

o 

1 

1 

1 

1 

1 

1 

1 

1 

fO 

>^_ 

ro 

*0> 

V 

r- 

o 

oo 

f^ 

o 

•o 

c 

^0 

rsj 

o 

-o 

(Nl 

•c 

00 

r- 

r- 

o- 

in 

^ 

fSJ 

f— 4 

■J- 

o 

■ 

-)f- 

• 

in 

o 

1 

*. 

«— 4 

■5fc 

o 

O- 

Q. 

Q. 

^ 

I 

3. 

vo 

00 

CD 

r 

X 

(T  -<  -^ 


00 

•o 

.—1 

1^ 

(*^ 

r-4 

(M 

^^ 

>o 

CO 

<M 

•o 

sD 

t\) 

(M 

CO 

>t 

-o 

r~ 

ro 

* 

>< 

(M 

>k 

<t 

— 

pn 

in 

(^ 

r- 

h- 

fM 

tNJ 

-^ 

^ 

r*- 

%. 

oo 

m 

o 

>^ 

rn 

o 

1^ 

in 

00 

in 

r^ 

v}- 

r- 

CO 

•o 

ro 

^J• 

^M 

r^ 

•c 

in 

>t 

00 

4- 

o 

in 

00 

•J- 

c> 

vO 

r^ 

o 

o 

in 

(M 

in 

^-4 

(T 

(M 

in 

O 

-*. 

r^ 

r 

#— « 

in 

CO 

vO 

(NJ 

tNJ 

CO 

(M 

ro 

• 

1 

• 

• 

• 

fl 

• 

• 

• 

• 

f 

o 

1 

X. 

<-4 

^ 

o 

<»^ 

o 

1 

—4 

1 

1—1 

rj 

O 

1 

00 

pa 

in 

o 

-J3 

-H 

>J- 

•O 

O 

(T> 

<NJ 

i\i 

• 

• 

• 

• 

in 

O 

r>- 

in 

1 

p^ 

m 

1 

Q. 

Q. 

CL 

<I 

<l 

< 

< 

a 

ii: 

.1 

<r 

a 

QC 

^ 

i 

h— 

^ 

^ 

Q. 

^ 

X 

H- 

J 

^ 

<r 

o 

CD 

CD 

►— 

i 

go 

o 

CD 

CD 

t— 

CO 

X 

X 

O 

CD 

CD 

CO 

X 

X 

o 

X 

z 

2 

t— 

X 

X 

X 

z 

TT 

t— 

227 


»T3 


l>0 

a 

Q.  CL 

a:       3 

I—  Q. 


i    2 


^ 
% 


O  <M 


— <       o       — <       o       o 

fM  >*•  J-  J-  O 

r~        o       r-        >t        ^t 


in 


i/> 


in        -J 


in 


-       X 


t\l 

m 

fM 

r- 

—4 

o 

iT 

CO 

in 

ro 

o 

ir 

fO 

vJ- 

C 

o 

O 

1 

o 

O 
1 

I 


Q. 

a 

Q- 

a. 

Q. 

<r 

<3 

<I 

< 

a 

Q. 

-J 

Qt 

^ 

JE 

<I 

Q. 

a. 

:«: 

Ji 

^ 

X 

o 

►— 

JC 

^ 

a 

^ 

I 

t— 

:s 

^ 

JC 

1/1 

uo 

O 

m 

CE 

^— 

J. 

1/1 

o 

GD 

CD 

OD 

CD 

CC 

CD 

I 

X 

O 

CD 

CD 

CD 

X 

X 

r 

X 

T 

X 

^ 

Z 

>— 

X 

X 

X 

?• 

z 

CvJ 


o 

<o 

m 

J- 

r\j 

J- 

o 

^ 

o 

(NJ 

■o 

CO 

o 

m 

^- 

4- 

O 

O 

o 

r~ 

* 

1^ 

m 

J|^ 

X 

in 

cr 

in 

O" 

rr, 

h- 

vO 

in 

* 

% 

•-4 

"k 

r>- 

o 

>0 

in 

o 

■XI 

O 

^i^ 

<f 

in 
1 

rsj 

>J- 

flO 

1 

1 

rvj 

1 

1 

ri^ 

(M 

• 
O 

n 

• 

00 

t 
o 

• 
O 

t 
(J- 

cr 

• 
o 

• 

1 

— < 
1 

1 

in 

rr. 
1 

1 

(NJ 

Ji< 

V 

•*r 

K 

•k 

11. 

V 

m 

^O 

r- 

INJ 

f\j 

0- 

00 

^co 

o 

^^ 

<M 

r- 

rr, 

>3- 

o 

in 

OD 

o 

• 

in 

o 

• 

• 

in 

• 

>3- 

• 

IM 

t 

• 

o 

^^ 

r- 

(NJ 

O 

.—1 

in 

(^ 

r- 

tNI 

(M 

r- 1 

cr 

f^ 

vj- 

O 

^^ 

in 

(NJ 

f\J 

^ 

^ 

r~- 

:k 

m 

in 

po 

«!■ 

fw 

(N) 

o 

o 

CD 

-¥- 

% 

^4 

^ 

ro 

o 

in 

-^ 

^^ 

O- 

o 

ro 

fo 

J- 

fM 

ro 

(T 

~M 

^ 

o 

(NJ 

(NJ 

<J- 

1 

1 

1 

1 

- 

1 

1 

I 

228 


00 

>t 

s. 

cr 

OO 

tr 

fSI 

--J 

00 

in 

o 

r^ 

«t 

t-i 

1 

^4 

^4 

*  *  *  *  + 

'O  .*       rvj  iri  00 

>}■  ITt         ,^  _<  .-I 

00  .— I       00  o  >*■ 


ro 


r~         n        CO 
>^         (\i        .^ 


>-D 


•T3 


<\* 

^-4 

in 

r-4 

o 

ro 

o 

00 

CO 

CT 

>o 

in 

^0 

rn 

>^ 

(T 

o 

>»• 

^t 

r>- 

■^ 

% 

*-4 

^ 

^^ 

o 

^J■ 

in 

r^ 

o 

cr 

•o 

:nj 

^ 

^ 

-4 

^ 

o 

st 

in 

1—4 

>!■ 

o 

r^ 

^ 

r- 

IT 

in 

o 

in 

cr 

CD 

CO 

o 

f^ 

•o 

1 

ro 

ro 

vO 

r— 1 

(N/ 

m 

•t 

in 

in 

vO 

1 

r^ 

h- 

1— t 
1 

t--4 

>!• 

>o 

* 

J^ 

^ 

i(. 

i' 

* 

> 

*- 

V 

+ 

^CT" 

h- 

f«^ 

^— t 

t- 

CO 

>o 

o 

CO 

00 

ca 

00 

h- 

in 

in 

r- 

>o 

o 

a 

f<> 

% 

^ 

in 

i 

^^ 

O 

vO 

>t 

(NJ 

CO 

o 

CT- 

-• 

% 

ic 

•o 

M 

in 

vD 

1 

m 

CO 

1 

CO 

0- 

cr 

1 

r\j 

in 

-4- 

o 

CT> 

rsj 

>t 

(\J 

00 

in 

>t 

CO 

o 

cr 

^-t 

in 

^^ 

in 

fSI 

CO 

■x- 

i^ 

r»- 

^ 

vO 

r- 

cr 

00 

1— 1 

cr 

r-» 

.t 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

% 

■X- 

00 

^ 

f\i 

in 

ro 

cr 

m 

rsj 

o 

rg 

in 

ro 

f- 

in 

>o 

sO 

•-4 

(VJ 

(T- 

<«t 

^•^ 

^-t 

t-l 

f— 1 

1 

(M 

00 

p-4 

1 

f^ 

1 

•* 

k 

it 

i 

;k 

fc 

* 

^ 

;*t 

> 

00 

o 

O 

O 

^-4 

^^ 

00 

f\J 

CO 

(N 

o 

vO 

rc 

(M 

CD 

J- 

<£) 

t>- 

OO 

•O 

■K- 

% 

^^ 

it 

f^ 

r^ 

O 

ro 

a> 

^^ 

O 

^4 

>J- 

% 

% 

rsj 

t 

in 

>o 

00 

^ 

o 

CO 

cr 

<M 

O 

1 

1— 1 

1 

t-l 

1 

1— < 

1 

^-1 

C\J 

1 

O 

CO 

>r 

cr 

in 

cr 

>J- 

^^ 

PO 

-^ 

sJ- 

r- 

CO 

<o 

<r 

r- 

>t 

^d 

.£) 

00 

^-4 

rr\ 

o 

r- 

r- 

r^ 

rri 

1^ 

^^ 

in 

r- 

p^ 

00 

r- 

tM 

>t 

in 

CO 

in 

^M 

<r 

r* 

00 

o 

PO 

o 

r-l 

ro 

■c 

vD 

•—1 

(\j 

in 

»-^ 

in 

1 

1 

in 

>-> 
1 

r^ 

•J- 

1 

UJ 

1/1 

o 

a.       Q- 

—        oi 
a:        13 


i 

I/) 

00 

o 

CO 

CD 

(C 

CD 

CO 

CO 

X 

X 

X 

I 

r 

X 

2 

Z 

< 

<r 

<i 

< 

< 

Ci. 

Of 

:^ 

3 

^ 

X 

>— 

i 

2: 

3 

1/1 

a 

CD 

CD 

00 

CD 

00 

X 

X 

X 

X 

X 

2 

z 

VITA 


229 


VITA 


Taminam  Zaki  Nakkash  was  born  May  4,  1942;  in  Beirut, 
Lebanon.   He  attended  elementary  school  there  and  was  gradu- 
ated from  Ali  Ibn  Abi  Taleb  Secondary  School  in  June,  1959; 
and  from  the  International  College  in  June,  1960. 

He  attended  the  American  University  of  Beirut,  Beirut, 
Lebanon,  in  September,  1960  majoring  in  civil  engineering. 
He  was  awarded  the  Bachelor  of  Engineering  degree,  with 
distinction,  in  June,  1964. 

From  July  to  December,  1964  he  was  employed  as  a  Con- 
struction Engineer  by  the  Trans -Arabi an  Pipe  Line  Company, 
in  Rafha,  Saudi  Arabia. 

He  was  employed  by  the  University  of  South  Carolina, 
Columbia,  South  Carolina,  as  a  Graduate  Research  Assistant, 
while  attending  graduate  school  there  from  January  1965  un- 
til he  was  av/arded  the  Master  of  Engineering  degree  in 
June  ,  1 956 . 

He  has  been  employed  as  a  Graduate  Instructor  by  Purdue 
University  since  August,  1966;  and  has  been  working  on  the 
Ph.D.  since  that  time. 

He  is  a  Junior  Member  of  the  Institute  of  Traffic 
Engineers,  an  Associate  Member  of  the  American  Society  of 
Civil  Engineers,  a  Student  Associate  Member  of  the  Operations 


230 


Research  Society  of  America,  a  Supporting  Member  of  the 
Highway  Research  Board,  and  a  Member  of  the  Order  of  Engi- 
neers and  Architects,  Beirut,  Lebanon.   He  is  a  member  of 
the  Tau  Beta  Pi  and  an  Associate  Member  of  the  Society  of 
Sigma  X  i . 

He  is  a  citizen  of  Lebanon. 


